AD/A-006  i»13 


i.o  at 


,,  Z£SJUZ2Z2£k  fLJJ* 

DOCUMENT  CONTROL  DATA  .RAD 

Security  elmt tttic mtt&n  **><*>•  of  •batrmet  *r>d  tndrutng  onnotstiem  wmf  b*  mntir+d  wh*.,  tS-  **ormtl  r+p*+1  t*  f Imairttod^ 

I 0«iCiN*rif»(i  ACTiviT*  (Corporate  eutbof)  J RC  POR  T SfC^NlTV  CL* SSlVi  C *T?ON 

„ . . UNCLASSIFIED 

Cambridge  Collaborative  — 

Cambridge,  Massachusetts 

I NE^OBt'tiTlC 

STATISTICAL  ENERGY  ANALYSIS  FOR  DESIGNERS,  PART  I.  Basic  Theory 


« Ot*C*M»T(ve  NOTH  (Typo  c t roport  mod  inclut ivo  dot**) 

Final  Report  - April  1972  to  December  1973 

'•T' WthoniJ)  (Ftrai  rtomo,  middle  intttof,  lm»t  nom*)  ~~  " 

Richard  H.  Lyjn 

I REPORT  DATE  " 

September  1974 

M CONTUAC  T 0«  GRANT  NO 

F33615-72-C-11S6 

b.  PROj C C T NO 

1370 

r Task  No.  137002 


|7«.  TOTAL  NO  OF  PACO 


7*.  no  or  7IFI 


M ORIGINATOR'S  REPORT  NuMBERlSl 


Cambridge  Collaborative  73-1 

•b.  OTHER  REPORT  now  (Any  otht  numbdrd  thst  m»yl 
thlm  report) 

AFFDL-TR- 74-56,  Part  I 


I"  Dir«l»UTIO«  lT»TtMtHT 


Approved  for  public  release;  distribution  unlimited. 


Ml  SoPPltMFNTART  NOTES 


112  SPONSORING  MILITARY  ACTlVITV 


Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  AFB,  Ohio 


This  report  describes  the  theoretical  development  of  Statistical 
Energy  Analysis  (SEA' . The  historical  development  of  the  subject  and 
an  extensive  annotated  bibliography  are  included  to  add  perspective 
to  the  technical  discussions.  Those  aspects  of  random  vibration  theory 
are  reviewed  first.  Then,  the  coupled  vibration  of  randomly  excited 
resonators  is  discussed.  Using  these  results  as  a basis,  the  SEA 
model  of  coupled  systems  is  introduced,  and  the  consequences  of  the 
model  are  derived.  From  these  basic  relationships,  a scheme  for  the 
prediction  of  vibration  levels  in  coupled,  randomly  excited, 
resonantly  vibrating  systems  is  developed. 


NATIONAI  1 fcLHNIC  AL 
INFORMATION  SERVICE 


DD  ,'“".“..1473 


BICES  SWiRT  TO  CHANS! 


UNCLASSIFIED 


Security  Classification 


Response  Prediction,  Statistical  Energy 
Anal /sis.  Preliminary  Design, 
Environmental  Estimation 


STATISTICAL  ENERGY  ANALYSIS  FOR  DESIGNERS 
PART  I.  BASIC  THEORY 


Richard  H.  Lyon 
Cambridge  Collaborative 


Approved  for  public  release;  distribution  unlimited. 

I'C 


FOREWORD 


The  research  described  in  this  report  was  performed  by 
Cambridge  Collaborative,  Cambridge,  Massachusetts  for  the 
Aerospace  Dynamics  Branch,  Vehicle  Dynamics  Division,  Air 
Force  Flight  Dynamics  Laboratory,  Wright-Patterson  Air  Force 
Base,  Ohio  under  Contract  F32615-72-C-1196,  Project  No.  1370, 
"Dynamic  Problems  in  Flight  Vehicles,"  and  Task  No.  137002 
"Flight  Vehicle  Vibration  Control."  Mr.  A.  R.  Basso  and 
later  Mr.  J.  D.  Willenborg  were  the  Project  Engineers  at  the 
Air  Force  Flight  Dynamics  Laboratory. 

Part  I of  this  report  covers  the  basic  theory  of 
statistical  energy  analysis  methods  while  Part  II  covers  the 
application  of  statistical  energy  analysis. 

The  report  covers  work  conducted  from  April  1972  to 
December  1973.  Cambridge  Collaborative  report  number  is 
73-1. 

This  report  was  released  by  the  author  in  December  1973. 


This  technical  report  has  been  reviewed  and  approved. 

/ ' ft-  / ' ' 

WALTER. J.  MYRYTOW 

Asst,  for  Research  & Technology 

Vehicle  Dynamics  Division 


ii 


ABSTRACT 


This  report  describes  the  theoretical  development  of 
Statistical  Energy  Analysis  (SEA) . The  historical  develop- 
ment of  the  subject  and  an  extensive  annotated  bibliography 
are  included  to  add  perspective  to  the  technical  discussions. 
Those  aspects  of  random  vibration  theory  are  reviewed  first. 
Then,  the  coupled  vibration  of  randomly  excited  resonators 
is  discussed.  Using  these  results  as  a basis,  the  SEA  model 
of  coupled  systems  is  introduced,  and  the  consequences  of 
the  model  are  derived.  From  these  basic  relationships,  a 
scheme  for  the  prediction  of  vibration  levels  in  coupled, 
randomly  excited,  resonantly  vibrating  systems  is  developed. 
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CHAPTER  1.  THE  DEVELOPMENT  OF  STATISTICAL  ENERGY  ANALYSIS 


1.0  Introduction 


This  report  is  a presentation  of  the  basic  theory  and  pro- 
cedure for  application  of  a branch  of  study  of  dynamical 
systems  called  Statistical  Energy  Analysis,  which  we  shall 
refer  to  as  "SEA".  The  name  SEA  was  coined  in  the  early 
1960 's  to  emphasize  certain  aspects  of  this  new  field  of 
study.  Statistical  emphasizes  that  the  systems  being  studied 
are  presumed  to  be drawn  from  statistical  populations  having 
known  distributions  of  their  dynamical  parameters.  Energy 
denotes  the  primary  variable  of  interest.  Other  dynamical 
variables  such  as  displacement,  pressure,  etc.,  are  found  from 
the  energy  of  vibration.  The  term  Analysis  is  used  to  empha- 
size that  SEA  is  a framework  of  study,  rather  than  a par- 
ticular technique. 

Statistical  Vibration  Analysis.  Statistical  approaches 
in  dynamical  analysis  have  a longbi story . In  mechanics,  we 
are  most  familiar  with  their  application  to  the  vibration  that 
is  random  in  time  of  a deterministic  system.  We  shall  use 
this  analysis  in  parts  of  Chapters  2 and  3.  It  is  useful  to 
emphasize  here  that  the  important  feature  of  SEA  is  the 
description  of  the  vibrating  system  as  a member  of  a statistical 
population  or  ensemble,  not  whether  or  not  the  temporal  be- 
havior is  random. 

Traditional  analyses  of  mechanical  system  vibration  of 
machines  and  structures  have  been  directed  at  the  lower  few 
resonant  modes  because  these  modes  tend  to  have  the  greatest 
displacement  response  in  many  instances,  and  the  frequencies 
of  excitation  were  fairly  low.  Of  course,  the  vibration  of 
walls  and  floors  and  their  high  frequency  sound  radiation  have 
been  of  interest  for  a long  time,  but  mechanical  and  structural 
engineers  have  been  generally  unaware  of  or  unconcerned  with 
this  work.  The  advent  of  fairly  large  and  lightweight  aero- 
space structures,  and  their  attendant  high  frequency  broad- 
band loads,  has  meant  much  more  attention  to  higher  order 
modal  analysis  for  purposes  of  predicting  structural  fatigue, 
equipment  failure  and  noise  production. 

A characteristic  of  higher  order  mode  analysis,  however, 
is  basic  uncertainty  in  modal  parameters.  The  resonance 
frequencies  and  mode  shapes  of  these  modes  show  great  sensi- 
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tivity  to  small  details  of  geometry  and  construction.  In 
addition , the  computer  programs  used  to  evaluate  the  mode  shapes 
and  frequencies  are  known  to  be  rather  inaccurate  for  the 
higher  order  modes,  even  for  rather  idealized  systems.  In 
light  of  these  uncertainties,  a statistical  model  of  the 
modal  parameters  seems  quite  natural  and  appropriate. 

If  there  is  cause  for  statistical  approaches  from  the 
nature  of  the  dynamical  problem,  there  is  equal  motivation 
from  the  viewpoint  of  application.  Mechanical  and  structural 
designers  are  often  faced  with  making  environmental  and 
response  estimates  at  a stage  in  a project  where  structural 
retail  is  not  known.  These  estimates  are  made  for  the 
qualification  of  equipments  and  the  design  of  isolation, 
damping,  or  structural  configurations  to  protect  equipment 
and  protect  the  integrity  of  the  structure.  Highly  detailed 
analyses  requiring  specific  knowledge  of  shape,  construction, 
loading  functions,  etc.,  are  not  appropriate.  Simpler 
statistical  analytical  estimates  of  response  to  environment 
that  preserve  parameter  dependence  (such  as  damping,  average 
panel  thickness,  etc.)  are  appropriate  to  the  designer’s 
need  at  this  stage. 

Inspirations  for  the  SEA  Approach.  There  is  experience 
in  dealing  with  dynamical  systems  described  by  random  para- 
meters.  Two  notable  examples  that  have  served  as  "touch- 
stones" in  early  developments  of  SEA  are  the  theory  of  room 
acoustics,  and  statistical  mechanics.  Room  acoustics  deals 
with  the  excitation  of  systems  of  very  many  degrees  of  free- 
dom (there  may  be  over  a million  modes  of  oscillation  of  a 
good  sized  room  in  the  audible  frequency  range)  and  the  inter- 
actions between  such  systems  (sound  transmission  through  a 
wall  is  an  example) . The  analyses  are  carried  out  using 
both  modal  and  wave  models.  The  very  large  number  of  degrees 
of  freedom  is  an  advantage  from  a statistical  viewpoint  — 
it  tends  to  diminish  the  fluctuations  in  prediction  of 
response.  We  shall  show  why  this  is  so  in  Chapter  4. 

Statistical  mechanics  deals  wit!  the  random  motion  of 
systems  with  either  a few  or  very  many  degrees  of  freedom. 
However,  it  is  random  motion  of  a very  special  type,  which 
we  may  call  "maximally  disordered."  In  this  state  of 
vibration,  all  modes,  whether  they  resonate  at  frequencies 
near  each  other  or  are  far  apart,  tend  to  have  equal  energy 
of  vibration  and  to  have  incoherent  motion.  We  must  add  the 
proviso,  "ignoring  quantum  effects."  The  energy  of  the 
modes  is,  aside  from  a universal  constant,  the  system  tem- 
perature. The  state  of  equal  modal  energy  is  spoken  of  as 
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"equipartition  of  energy".  In  SEA,  we  sometimes  make  the 
squipartition  assumption  for  modes  that  resonate  in  the 
same  frequency  band,  but  not  for  all  modes. 

Statistical  mechanics,  and  its  related  science,  heat 
transfer,  also  teach  us  that  thermal  (random  vibration) 
energy  flows  from  hotter  to  cooler  systems,  and  that  the  rate 
of  flow  is  proportional  to  temperature  (modal  energy)  dif- 
ference. In  Chapter  3,  we  show  that  this  result  also  applies 
to  dynamical  systems  excited  by  broad  band  noise  sources. 

Not  only  that,  but  since  we  also  show  in  Chapters  3 and  4, 
that  narrow  band  sources  are  equivalent  to  broad  band  sources 
when  system  averages  are  taken,  the  result  can  be  generalized, 
with  proper  care,  to  pure  tones. 

The  study  of  the  statistical  mechanics  of  electrical 
circuits  shows  that  a resistor  at  known  temperature  is 
equivalent  to  a thermal  reservoir  (or  temperature  bath) . The 
interaction  of  the  dynamical  system  with  this  reservoir  is 
represented  as  a white  noise  generator  in  circuit  theory.  We 
can  turn  the  argument  around  and  say  that  a damper  (mechanical 
resistance  element),  in  conjunction  with  a noise  generator, 
represents  a thermal  reservoir,  and  we  should  not  be  surprised 
when  certain  "thermal"  results  develop  from  its  analysis. 

The  Advantages  and  Limitations  of  Statistical  Analysis. 

One  advantage  of  statistical  analysis  of  systems  may  be  seen 
from  the  practical  aspects  of  room  acoustics.  If  one  truly 
has  10 6 modes  to  deal  with,  even  the  m.s.  pressure  associated 
with  ^ach,  changing  with  time  as  the  flute  gives  way  to  the 
tympani,  would  be  a hopeless  mass  of  information  to  assimilate. 
What  one  does  instead  is  to  describe  the  field  by  a few 
coherent  features  of  the  modal  pattern  (direct  field  and  a 
few  early  reflections)  and  the  incoherent  energy  (reverberant 
fielc'  totalled  into  a few  frequency  bands.  Thus,  instead  of 
10 6 measures  of  the  sound  field  (which  would  be  incomplete  in 
any  case  without  the  coherence  data)  , we  are  able  to  describe 
the  sound  field  ef fective  y by  10-20  measures. 

The  statistical  analysis  also  allows  for  much  simpler 
description  of  the  system,  whether  one  describes  the  field 
by  modes  or  waves.  In  the  former  case,  modal  density, 
average  mod.'l  damping,  and  certain  averages  of  modal  impedance 
to  sound  sources  are  required.  In  a wave  description,  such 
parameters  as  mean  free  path  for  waves,  surface  and  volume 
absorption,  and  general  geometric  configuration  are  required. 
The  number  of  input  parameters  is  generally  in  balance  with, 
the  number  of  measures  (10-20)  to  be  taken. 
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The  most  obvious  disadvantage  of  statistical  approaches 
is  that  they  give  statistical  answers , which  are  always  subject 
to  some  uncertainty.  In  very  high  order  systems,  this  is  not 
a great  problem.  Many  of  the  systems  we  may  wish  to  apply 
SEA  to,  however,  may  not  have  enough  modes  in  certain  frequency 
bands  to  allow  predictions  with  an  acceptable  degree  of  cer- 
tainty. To  keep  track  of  this,  we  may  calculate  the  variance 
as  well  as  the  mean,  and  also  calculate  the  confidence  for 
prediction  intervals.  This  is  discussed  in  Chapter  4. 

In  addition  to  hard  and  fast  computational  problems, 
there  are  certain  difficulties  in  the  psychology  of  statistical 
methods.  A designer  is  not  dealing  with  a gas  of  complicated 
molecules  in  random  collision  — he  is  concerned  about  pre- 
dicting the  structural  response  of  a wing,  for  which  he  has 
engineering  drawings,  to  a loading  environment,  for  which  he 
has  flight  data.  Instead  of  following  a deterministic  cal- 
culation (probably  computer  based)  it  is  suggested  that  he 
will  get  a "better"  estimate  if  he  represents  the  wing  as  a 
flat  plate  of  a certain  average  thickness  and  total  area! 

His  incredulity  may  be  imagined.  But  he  must  remember  that 
his  knowledge  of  the  wing  at  the  50th  mode  of  vibration  may 
be  just  as  well  represented  by  the  flat  plate  as  it  is  by 
his  drawings.  Also,  the  answers  he  gets  by  SEA  will  be  in  a 
form  that  is  usable  to  him,  generally  retaining  parameter 
dependence  that  will  allow  him  to  interpret  the  effect  of 
certain  simple  design  changes  on  response  level. 

To  close  this  introductory  section,  it  may  be  revealing 
to  quote  from  M.L.  Mehta,  a theoretical  physicist  concerned 
with  applying  statistical  method.*,  to  large  nuclei. 

"Statistical  theory  does  not  predict  the  detailed  level  se- 
quence (resonance  frequencies]  of  any  one  nucleus,  but  it 
does  describe  the... degree  of  irregularity .. .that  is  expected 
to  occur  in  any  nucleus. .. .Here  we  shall  renounce  knowledge 
of  the  system  itself .. .There  is  a reasonable  expectation, 
though  no  rigorous  mathematical  proof,  that  a system  under 
observation  will  be  described  correctly  by  an  ensemble 
average....  If  this  particular  [system]  turns  out  to  be  far 
removed  from  the  ensemble  average,  it  will  show  that... [the 
system] .. .possesses  specific  properties  of  which  we  are  not 
aware.  This,  then.,  will  prompt  us  to  try  to  discover  the 
nature  and  origin  c>f  these  properties."  The  problems  we 
face  seem  to  be  universal. 
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A Brief  Historical  Survey 


Beginnings.  The  earliest  work  in  the  development  of 
SEA  as  an  identifiable  entity  were  two  independent  calculations 
in  1959  by  R.  H.  Lyon  and  P.  W.  Smith,  Jr.  In  England,  on 
an  NSF  postdoctoral  fellowship,  Lyon  calculated  the  power 
flow  between  two  lightly  coupled,  linear  resonators  excited 
by  independent  white  noise  sources.  He  found  that  the  power 
flow  was  proportional  to  the  difference  in  uncoupled  energies 
of  the  systems  and  that  it  always  flowed  from  the  system  of 
higher  to  lower  modal  energy. 

The  other  calculation  wa3  by  Smith  at  BBN  working  under 
U.S.  Air  Force  support.  Smith  calculated  the  response  of  a 
resonator  excited  by  a diffuse,  broad  band  sound  field,  and 
found  that  the  response  of  the  system  reached  a limit  when 
the  radiation  damping  of  the  resonator  exceeded  its  internal 
damping,  but  that  this  limit  did  not  depend  on  the  precise 
value  of  the  radiation  damping. 

This  result  of  Smith's  was  somewhat  surprising  since 
many  workers  regarded  an  acoustic  noise  field  simply  as  a 
source  of  broad  band  random  excitation.  When  a resonator, 
excited  by  broad  band  noise,  has  its  internal  damping  reduced 
to  zero,  the  response  diverges,  i.e.,  goes  to  infinity.  The 
limit  involved  in  Smith's  result  was  of  course  due  to  the 
reaction  of  the  sound  field  itself  on  the  resonator. 

After  Lyon  joined  BBN  in  the  fall  of  1960,  it  developed 
that  Smith's  limiting  vibration  amounted  to  an  equality  of 
energy  between  the  resonator  and  the  average  modal  energy  of 
the  sound  field.  The  two  calculations  were  consistent,  and 
power  would  flow  between  resonators  until  equilibrium  would 
develop.  If  the  coupling  were  strong  enough  compared  to 
internal  damping,  equipartition  would  result.  But  how 
specifically  did  the  wave-field  result  of  Smith  relate  to 
the  two-mode  interaction  that  Lyon  had  studied? 

To  answer  this  question,  Lyon  and  Maidanik  wrote  the 
first  paper  that  may  be  said  to  be  an  SEA  publication,  although 
the  name  SEA  had  not  then  been  coined.  This  paper  combined 
Lyon's  work  in  England  with  extensions  to  make  it  able  to  deal 
with  the  kind  of  problem  Smith  had  analyzed.  Formulas  for 
the  interaction  of  a single  mode  of  one  system  with  many 
modes  of  another  were  developed , and  experimental  studies  of 
a beam  (few  mode  system)  with  a sound  field  (multi  modal 
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system)  were  reported.  This  work  also  showed  the  importance 
of  the  basic  SEA  parameters  for  response  prediction:  modal 
density,  damping,  and  coupling  loss  factor. 
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Early  Extensions  and  Improvements  of  the  Theory. 

Almost  immediately,  the  activity  in  SEA  spilt  along two 

lines.  One  line  was  the  clarification  of  basic  assumptions 

and  improvement  in  the  range  of  approximation  to  real 

system  performance.  The  second  line  was  the  application  cf 

SEA  to  other  systems.  The  earliest  application  was  to  sound- 

structure  interaction,  largely  a result  of  Smith’s  work,  but 

also  because  it  seemed  "obvious"  that  SEA  would  work  best  < 

when  a sound  field,  with  all  of  its  many  degrees  of  freedom 

was  involved.  Very  soon,  however,  applications  were  also  made 

to  structure-structure  interactions. 

One  basic  assumption  in  SEA  had  to  deal  with  light  i 

coupling.  How  much  of  a restriction  did  this  represent?  1 

Also,  what  were  the  uncoupled  systems?  Two  independent 
studies,  by  Ungar  and  Scharton,  showed  that  the  light 
coupling  assumption  was  unnecessary  if  the  uncoupled  systems 
were  defined  as  the  blocked  system,  meaning  that  the  other 
system  was  held  fixed  while  the  system  being  considered  was 
allowed  to  vibrate.  Also,  the  energy  flow  relations  were 
valid  whether  one  was  using  the  "blocked"  energies  of  the 
system  as  the  driving  force  or  the  actual  energy  of  each 
subsystem  with  the  coupling  Intact.  Of  course,  the  con- 
stant of  proportionality  would  be  different  for  these  two 
calculations. 

The  question  of  quantifying  the  uncertainty  in  the 
prediction  of  energy  flow  was  examined  by  Lyon,  who  developed 
a theory  of  response  variance  and  prediction  intervals  for 
SEA  calculations.  The  calculation  of  variance  for  structure- 
structure  interactions  in  which  relatively  few  modes 
participate  in  the  energy  sharing  process  was  first  included 
in  a paper  by  Lyon  and  Eichler. 

An  important  extension  of  the  two  system  theory  was 
made  by  Eichler,  who  developed  predictions  from  energy 
distribution  for  three  systems  connected  in  tandem.  The 
practical  application  of  SEA  often  involves  the  flow  of 
vibratory  energy  through  several  intervening  "substructures" 
before  it  gets  to  a vibration  sensitive  area.  It  is 
important,  therefore,  to  be  able  to  predict  the  energy  dis- 
tribution for  such  cases. 
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Improvements  in  Range  of  Application.  The  earliest 
work  on  structure- structure  vibration  transmission  was  under- 
taken with  Air  Force  sponsorship  and  was  concerned  with 
electronic  package  vibration.  An  early  paper  on  this  topic 
by  Lyon  and  Eichler  dealt  with  plate  and  beam  interaction, 
very  similar  to  an  example  discussed  in  Chapter  4,  and  two 
plates  connected  together.  A subsequent  paper  on  a three 
element,  plate-beam-plate  system  by  Lyon  and  Scharton  made 
use  of  the  earlier  formulation  of  three  element  systems  by 
Eichler  and  some  plate  edge-admittance  calculations,  also 
by  Eichler. 

The  basic  SEA  theory  was  pretty  much  directly  applicable 
to  these  new  systems.  The  major  problem  was  evaluation  of 
modal  densities  and  coupling  loss  factors  for  various  inter- 
acting junctions  between  systems.  For  example,  the 
radiation  of  sound  by  reinforced  plates  was  evaluated  by 
Maidanik,  and  a similar  study  of  the  radiation  of  sound  by 
cylinders  was  undertaken  by  Manning  and  Maidanik.  We  have 
already  mentioned  the  plate-edge  admittance  calculations  of 
Eichler.  These,  along  with  earlier  (pre-SEA)  calculations 
of  force  and  moment  impedance*  of  beams  and  plates  have 
allowed  a wide  variety  of  structural  coupling  loss  factors 
to  be  evaluated.  Quite  *ec-*mtly,  a series  of  soil-foundation 
impedance  evaluations  by  Kurzweil  allows  one  to  apply  SEA  to 
certain  structure-ground  vibration  problems. 

Modal  densities  of  acoustical  spaces  have  been  studied 
for  quite  a long  time.  Also,  the  modal  densities  of  some 
flat  and  curved  panel  structures  pre-date  SEA.  However,  the 
activity  in  SEA  has  motivated  work  in  modal  density  evaluation. 
For  example,  the  modal  density  of  cylinders  has  been  studied 
by  Heckl,  Manning,  Chandiramani,  Miller,  and  Szechenyi.  The 
modal  density  of  cones  has  been  calculated  by  Chandiramani, 
and  by  Wilkinson  for  curved  sandwich  panels.  Generally, 
modal  density  prediction  has  not  been  as  difficult  as  the 
calculation  of  coupling  loss  factors. 

Understanding  and  prediction  of  system  damping  has 
improved  very  little  over  the  years  since  SEA  began.  In  most 
part,  the  improvements  that  have  taken  place  were  not 
particularly  related  to  SEA  work,  although  the  work  by  Heckl 
on  plate  boundary  absorption  and  by  Ungar  and  Maidanik  on 
air  pumping  along  riveted  beams  were  generally  related  to 
SEA.  Despite  this  work,  our  ability  to  predict  damping  in 
built-up  structures  is  still  based  largely  on  empiricism. 


7 


Other  Developments . We  should  also  note  certain 
developments  that  might  be  termed  "sword  into  plowshares" 
activity.  Most  early  applications  of  SEA  have  been  aerospace 
related  because  that  is  where  the  problems  have  been  (and 
still  are)  and  DOD  and  NASA  were  paying  the  bills.  There 
have  been  recent  applications  of  SEA  to  certain  architectural 
and  building  problems  that  are  worthy  of  note.  In  England, 
Crocker  and  his  associates  have  carried  out  studies  of  sound 
transmission  through  double  septum  walls,  following  up  on 
some  earlier  work  by  White  and  Powell.  More  recently,  Ver  has 
applied  SEA  to  the  prediction  of  sound  transmission  through 
floors  consisting  of  a main  slab  with  another  lighter  slab 
floated  on  many  small  point  springs.  Rinsky  has  also  studied 
sound  transmission  through  double  and  single  stud-reinforced 
walls  using  SEA. 

Finally,  the  attempts  to  better  understand  the  theoretical 
basis  for  SEA  and  the  limiting  effect  of  its  assumptions  on 
the  range  of  applications  has  continued.  The  most  notable 
effort  along  these  lines  has  been  by  Bogdanoff  and  Zeman. 

Other  work  includes  a recent  thesis  by  Lotz  and  calculations 
comparing  ideal  deterministic  with  averaged  systems  by 
Scharton,  Manning  and  Remington.  So  far,  the  attempts  at 
further  reducing  the  number  of  assumptions  that  one  must  make 
in  SEA  have  not  been  very  effective,  but  we  must  not  be  dis- 
couraged by  this.  Every  step,  even  the  small  ones  taken  to 
this  point,  have  been  very  productive  in  extending  the  reach 
and  usefulness  of  SEA. 


1 . 2 The  General  Procedures  of  SEA 


In  the  following  chapters,  we  will  derive  the  basic 
equations  of  SEA  and  give  motivation  for  the  modeling  and 
computational  procedures.  To  provide  an  overview  of  this 
process,  however,  in  this  section  we  describe  the  way  that 
SEA  calculations  are  made,  and  the  steps  that  are  necessary 
to  arrive  at  a prediction  of  response.  Hopefully,  this 
will  provide  a framework  for  a better  understanding  of  the 
purpose  of  the  later  chapters  and  how  the  various  elements 
of  SEA  fit  together. 

Model  Development.  In  its  simplest  elements,  SEA 
results  in  a procedure  for  calculating  the  flow  and  storage 
of  vibrational  energy  in  a complex  system.  The  energy 
storage  elements  are  groups  of  "similar  modes".  Energy 


input  to  each  of  the  storage  elements  comes  from  a set  of 
external  (usually  random)  sources.  Energy  is  dissipated  by 
mechanical  damping  in  the  system,  and  transferred  between 
the  storage  elements.  The  analysis  is  essentially  that  of 
linear  R-C  circuit  theory  with  energy  playing  a role 
analogous  to  electric  charge  and  modal  energy  taking  the 
role  of  electrical  potential.  A typical  SEA  model  is  drawn 
in  Fig.  1.1  and  the  analogous  electrical  circuit  that  might 
be  used  to  represent  it  is  shown  in  Fig.  1.2.  We  shall  not 
make  direct  use  of  this  electrical  analogy  in  this  report, 
but  rather  directly  work  with  the  simultaneous  equations. 

The  fundamental  element  of  the  SEA  model  is  a group  of 
"similar"  energy  storage  modes.  These  modes  are  usually 
modes  of  the  same  type  (flexural,  torsional,  etc.)  that 
exist  in  some  section  of  the  system  which  we  may  call  a "sub- 
system” (an  acoustic  volume,  a beam,  a bulkhead,  etc.).  In 
selecting  the  modal  group,  we  are  concerned  that  it  meets 
the  criteria  of  similarity  and  significance.  Similarity 
means  that  we  expect  the  modes  of  this  group  to  have  nearly 
equal  excitation  by  the  sources,  coupling  to  modes  of  other 
subsystems,  and  damping.  If  these  criteria  are  met,  they 
will  alro  have  nearly  equal  energy  of  vibration.  Significance 
means  that  they  play  an  important  role  in  the  transmission, 
dissipation,  or  storage  of  energy.  Inclusion  of  an 
"insignificant"  modal  group  will  not  cause  errors  in  the  cal- 
culations, but  may  needlessly  complicate  the  analysis. 

Input  power  from  the  environment,  labelled  II.  , may 
result  from  a turbulent  boundary  layer,  acoustical  noise,  or 
mechanical  excitation.  It  is  usually  computed  for  some  fre- 
quency band,  possibly  a one-third  or  full  octave  band.  In 
order  to  evaluate  U.  , we  need  to  know  certain  input 
impedances  for  the  subsystem.  The  important  requirement  is 
that  this  input  power  not  be  sensitive  to  the  state  of 
coupling  between  subsystems.  If  it  is,  then  the  system  pro- 
viding the  power  (a  connecting  structure  for  example)  has 
important  internal  dynamics  and  should  be  modeled  as  another 
subsystem. 

The  dissipation  of  power  for  each  subsystem  IK. 
represents  energy  truly  lost  to  the  mechanical  vibration  and 
will  depend  only  on  the  amount  of  energy  stored  in  that  sub- 
system. It  may  be  truly  dissipated  by  friction  or  viscosity, 
or  it  may  be  merely  radiated  away  into  the  air  or  surrounding 
structure.  The  important  proviso  is  that  this  power  cannot 
be  returned  to  the  system.  If  it  can,  then  it  is  part  of  the 
power  flow  through  coupling  elements  and  will  require  the 
addition  of  another  subsystem  or  coupling  path  to  the  diagram. 

9 


a«<  4.*— V—1 * X* 


» pi, 1 • 


The  transmitted  power  II,  2 , represents  the  rate  of  energy 
exchange  between  subsystems  1 and  2.  All  energy  quantities 
that  we  deal  with  here  arr  time  averaged  values.  There  may 
be  very  large  temporal  variations  in  power  flow  between  the 
subsystems , even  larger  than  the  average  flow,  but  these  are 
ignored.  The  transmitted  power  depends  on  the  difference  in 
modal  energy  of  the  two  subsystems  and  the  strength  of  the 
coupling  between  them. 

Parameter  Evaluation.  Evaluation  of  the  quantities 
that  appear  In  fig.  1.1  and  that  are  discussed  in  the  pre- 
ceding paragraphs  require  the  evaluation  of  certain  para- 
meters , which  we  may  call  SEA  parameters,  that  mostly  pre- 
date SEA.  We  group  them  as  "energy  storage"  and  "energy 
transfer"  parameters.  Energy  storage  is  determined  by  the 
number  of  available  modes  N^,  N2...for  each  subsystem  in  the 
chosen  frequency  band  Aw.  The  ratio  of  N to  Aw  i3  called 
the  modal  density  n of  the  subsystem  and  is  frequently  used 
in  SEA  calculations  instead  of  the  mode  count  N. 

Energy  transfer  parameters  include  the  input  impedance 
to  the  system  for  the  determination  of  input  power,  the  loss 
factor,  which  relates  subsystem  energy  to  dissipated  power, 
and  the  coupling  loss  factor  relating  transmitted  power  to 
subsystem  modal  energy.  In  the  following  paragraphs  we  pro- 
vide a brief  indication  of  how  these  parameters  are  usually 
evaluated . 

Modal  density  may  be  measured  by  exciting  the  subsystem 
with  a pure  tone  and  gradually  varying  the  frequency.  As  a 
resonance  is  encountered  a maximum  in  response  will  occur. 

If  a chart  of  response  amplitude  as  a function  of  frequency 
is  drawn,  these  peaks  may  be  counted.  This  technique  may 
miss  some  modes  and  is,  therefore,  not  perfectly  reliable. 
Experimental  methods  have  been  devised  to  reduce  the  number 
of  modes  missed,  but  the  error  cannot  be  completely 
eliminated. 

The  most  commonly  used  way  of  evaluating  modal  density  is 
simply  to  calculate  it  from  theoretical  formulas.  Most 
systems  have  modal  densities  that  may  be  calculated  in  terms 
of  relatively  simple  gross  parameters , such  as  overall 
dimensions,  and  the  average  speed  o.!*  waves  in  the  system. 

A few  examples  of  the  calculation  of  modal  density  are  given 
in  Chapter  2. 
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The  dissipation  of  energy  is  measured  by  the  loss 
factor  of  the  system  (defined  in  Chapter  2).  Unless  the  dis- 
sipation mechanism  is  of  a very  particular  type,  it  is  more 
reliable  to  measure  the  loss  factor  than  it  is  to  calculate 
it.  This  is  done  by  measuring  the  rate  of  energy  decay  in 
the  system  when  the  excitation  has  been  removed,  or  by 
measuring  the  response  bandwidth  of  individual  resonance 
peaks.  The  relation  of  both  of  these  measures  to  the  loss 
factor  is  developed  in  Chapter  2. 

The  input  power  from  the  environment  is  sometimes 
measured,  but  more  often  calculated.  It  may  be  measured  by 
isolating  the  subsystem  of  known  damping  (loss  factor)  and 
observing  its  response  to  the  environment.  By  equating  input 
to  dissipated  power,  IIin  is  known.  The  input  power  may  be 
calculated  by  evaluating  the  load  exerted  on  the  subsystem 
by  the  environment  and  the  impedance  of  the  subsystem  to  this 
load.  Such  calculations  may  be  quite  involved,  but  a simple 
example  is  given  in  Chapter  2 for  flat  plates  excited  by  a 
point  force. 

The  coupling  loss  factor  is  the  parameter  governing 
transmitted  power.  It  is  defined  in  Chapter  3 and  has  been 
measured  for  some  systems,  although  it  is  often  calculated. 

It  is  also  related  to  quantities  that  may  have  been  cal- 
culated or  measured  for  other  reasons  - the  junction 
impedance j of  mechanical  systems,  the  transmission  loss  of 
walls,  the  acoustic  radiation  resistance  of  a piston,  and 
other  similar  measures.  Since  it  in  general  depends  on  both 
subsystems,  and  the  variety  of  systems  of  application  for 
SEA  is  very  large,  the  number  of  coupling  loss  factors  that 
we  may  be  concerned  with  increases  as  the  square  of  the 
number  of  subsystems  on  the  list.  The  better  strategy  would 
appear,  therefore,  to  express  the  coupling  loss  factor  where 
possible  in  terms  of  subsystem  impedances,  as  we  have  done  in 
Chapter  3,  and  then  list  the  impedances. 

Calculation  of  Response.  After  the  model  has  been 
decided  upon  and  the  parameters  evaluated,  the  final  step 
is  to  calculate  the  response.  The  first  part  of  this  cal- 
culation is  the  evaluation  of  the  vibrational  energy  in  the 
various  energy  storage  elements  or  mode  groups.  As  noted 
above,  this  amounts  to  solving  a set  of  linear  algebraic 
equations,  one  equation  to  each  subsystem.  A simple 
example  is  worked  out  in  Chapter  3.  These  average  energies 
will  depend  on  the  various  input  power  values,  the  modal 
densities,  and  the  coupling  and  dissipative  loss  factors. 
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On  the  basis  of  the  formalism  adopted  in  Chapters  3 and  ] 

4 , the  most  direct  relation  between  vibrational  energy  and 
another  dynamical  variable  is  with  the  velocity  of  motion.  ’3 

By  formulas  developed  in  the  analysis , the  velocity  may  in 
turn  be  related  to  still  other  variables  — displacement,  i 

strain,  stress,  pressure,  etc.  These  variables  are  still 
forms  of  a spatially  averaged  response.  However,  on  the 
basis  of  system  geometry  and  mode  shapes , estimates  earn  also  j i 

be  made  of  the  spatial  distribution  of  response.  j ! 


I The  variability  just  referred  to  occurs  even  with  the 

SEA  assumption  of  equal  energy  for  every  mode  in  the  sub- 
system. In  addition,  however,  there  is  variability  because 
we  deal  with  a particular  system  in  the  laboratory,  not  an 
ensemble  of  systems.  As  noted  in  the  introduction  to  this 
| chapter,  the  ensemble  average  will  not  fit  each  member 

| exactly,  and  we  may  expect  some  variation  in  the  various  para- 

I meters  that  we  have  been  discussing.  The  analysis  of  these 

* variations  is  an  important  part  of  the  response  prediction 

| process  and  tells  us  how  much  reliance  we  may  place  on  a 

• response  estimate  based  on  the  average  behavior.  The  dis- 

| cussion  of  this  topic  is  developed  in  Chapter  4. 


1 . 3 Future  Developments  of  SEA 


The  future  development  of  SEA  is,  in  large  part,  likely 
to  be  a continuation  of  certain  features  of  its  development 
to  date,  i.e.,  its  application  to  a larger  group  of  sub- 
systems. This  group  will  surely  include  new  structures  and 
seismic  (ground  vibration)  systems,  and  possibly  water 
waves  and  ship  or  offshore  structures.  Such  applications 
will  increase  the  glossary  of  SEA  parameters  and  make  SEA 
even  more  useful  than  it  has  been  until  now. 


The  wider  use  of  SEA  is  also  likely  to  involve  a larger 
group  of  professionals  within  each  subject  area.  This  means 
that  research  workers  and  analysts,  test  engineers,  and 
designers  are  all  likely  to  use  it  as  one  tool  among  many 
that  are  available  for  dynamical  analysis.  We  may  expect  to 
see  SEA  become  more  computerized,  not  only  for  the  purpose 
of  solving  the  simultaneous  equations  governing  response 
energy,  but  also  for  evaluating  coupling  loss  factors  for 
very  complex  structures  using  finite  element  methods. 
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Finally,  we  may  expect  to  see  developments  in  the  basic 
theory  of  the  statistical  analysis  of  systems  that  will 
greatly  expand  the  conceptual  framework  of  SEA.  To  be  more 
specific,  let  us  see  how  SEA  has  expanded  upon  the  framework 
of  statistical  mechanics,  and  what  may  develop  in  this  regard 
in  the  future. 

Statistical  mechanics  represents  the  most  disordered 
state  of  motion  possible  for  a system.  Every  mode  through- 
out the  system  has  the  same  average  energy.  This  includes 
modes  in  different  subsystems  and  those  that  resonate  at 
different  frequencies.  Thermodynamics  is  the  macroscopic 
counterpart  to  this  totally  equilibrium  system  of  analysis. 

The  possibility  for  different  temperatures  in  different 
parts  of  a system  is  allowed  in  the  kinetic  theory  version 
of  statistical  mechanics  and  in  "non-equilibrium"  thermo- 
dynamics. The  modes  within  each  subsystem  all  have  the  same 
energy,  but  the  average  modal  energies  of  the  subsystems 
may  be  different. 

This  latter  situation  is  the  one  we  deal  with  in  SEA, 
except  that  we  proceed  even  farther  from  equilibrium  by 
allowing  modes  in  different  frequency  bands  to  have  differing 
energies.  Basically,  we  can  get  away  from  this  because  the 
systems  that  we  deal  with  are  presumed  linear,  and  modes  that 
resonate  in  different  frequency  bands  may  be  considered  to 
be  uncoupled  from  each  other.  Thus,  SEA  represents  still 
another  step  away  from  the  "maximally  disordered"  state 
described  by  statistical  mechanics. 

At  this  point,  it  is  logical  to  ask  whether  there  is 
another  step  to  be  made  in  this  sequence  toward  a model 
that  is  less  disordered  than  the  current  SEA  model  that  would 
provide  use  with  useful  answers  in  situations  in  which  SEA 
has  limitations.  One  answer  might  be  — forget  about  this 
chain  of  logic,  go  all  tl  j way  back  to  the  deterministic 
system.  That  is  a possible  answer,  but  it  may  not  be  the 
most  useful  one. 

The  most  glaring  deficiency  of  SEA  is  its  inability  to 
deal  with  modal  coherence  effects.  In  simplest  terms,  modal 
coherence  may  lead  to  such  phenomena  as  "direct  waves," 
discussed  in  Chapter  2.  It  is  not  clear  at  present  whether 
this  represents  the  "next  logical  step"  in  the  chain  that 
we  spoke  of,  but  it  bears  examination. 
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A second  deficiency  of  SEA  is  its  assumption  that  all 
inodes  that  resonate  in  Aw  are  equally  probable  over  that 
interval,  and  that  their  resonance  frequencies  are  unaffected 
by  the  values  of  the  resonance  frequencies  of  the  other  modes.  .ij 

There  have  already  been  some  developments  on  improving  upon  \ 

this  hypothesis  and  we  may  expect  more.  We  are  on  more  solid  \ 

ground  in  predicting  developments  here  than  we  are  in  our  j 

concern  for  modal  coherence  effects.  j 


1.4  Organization  of  the  Report 


This  report  is  concerned  with  presenting  the  basic 
theoretical  elements  of  SEA.  We  may  cite  these  as: 

(1)  The  theory  of  energy  exchange  between  two 
resonators.  This  is  a problem  in  random 
vibration  theory. 

(2)  The  extension  of  the  two  resonator  inter- 
action to  the  interaction  of  two  systems 
having  modal  behavior.  This  extension 
requires  that  we  introduce  the  basic  idea 
of  statistical  populations  of  systems. 

(3)  The  representation  of  system  interactions 

by  the  SEA  model  with  the  attendant  necessity 
to  identify  and  evaluate  SEA  parameters. 

(4)  The  calculation  of  average  energy.  This  is 
exemplified  by  applying  the  formulas  to  some 
fairly  simple  cases. 

(5)  The  calculation  of  variance  in  response  and 
the  use  of  variance  to  generate  estimation 
intervals  and  their  associated  confidence 
coefficients. 

Since  the  emphasis  in  the  report  is  the  explication  of 
SEA,  we  have  included  generous  discussion  of  the  ideas  of 
energy  storage  and  transfer,  statistical  ensembles  of  systems, 
isolated  vs.  coupled  systems,  etc.  Our  intention  is  to  pro- 
vide the  most  clear  cut  statement  on  these  matters  that  can 
be  made  at  present.  The  reader  will  have  to  decide  how  well 
we  have  succeeded. 
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Since  our  attention  is  on  SEA,  however,  we  do  not  spend 
time  on  other  subjects  that  we  employ  in  the  various 
examples.  Thus,  for  example,  we  use  the  equations  of  bending 
motion  of  plates  and  beams,  but  we  do  not  derive  them.  We 
also  use  certain  input  impedances  to  beams,  plates  and  sound 
fields,  but  these  are  not  derived.  For  the  most  part,  such 
derivations  are  readily  available  to  the  student,  and  in 
any  case,  in  the  annotated  bibliography,  guidance  is  given 
to  where  this  information  may  be  found. 

The  reader  will  notice  that  no  references  are  given  in 
the  text,  which  is  an  unusual  procedure  for  a report.  In 
adopting  this  procedure,  we  do  not  wish  to  deny  credit  to 
anyone,  rather  we  have  followed  regular  textbook  practice  of 
providing  a bibliography,  but  not  breaking  the  presentation 
with  references.  Partly  this  procedure  is  a matter  of  self- 
discipline;  sometimes  it  is  very  tempting  for  an  author  to 
try  to  avoid  a difficult  matter  by  slipping  in  a reference. 

Chapter  2 deals  with  some  fundamental  ideas  in  the 
energy  of  vibration  of  single  and  multi  degree  of  freedom 
systems.  The  storage  of  kinetic  and  potential  energy  by 
modes  in  free  and  forced  vibration,  and  the  decay  or  rate 
of  energy  removal  due  to  damping  are  discussed,  partly  as  a 
review  of  this  important  subject  and  partly  to  establish 
some  of  the  basic  vocabulary  of  the  work  to  follow.  Also, 
since  both  modal  and  wave  descriptions  of  vibration  are 
employed,  energy  storage,  dissipation,  and  flow  in  a wave 
description  are  reviewed  also. 

Chapter  3 is  concerned  with  developing  the  theory  of 
average  power  flow  (average  in  both  a temporal  and  ensemble 
sense)  between  single  and  multi  degree  of  freedom  systems. 
First,  power  flow  between  two  simple  resonators  is  calculated. 
The  important  ideas  of  a blocked  system  is  introduced  and 
power  flow  is  calculated  in  terms  of  both  blocked  and  coupled 
system  energies.  The  idea  of  averaged  modal  interaction  as 
a white  noise  source  is  also  introduced. 

The  latter  sections  of  Chapter  3 are  the  central  to  our 
discussion.  The  very  important  ideas  of  blocked  systems, 
ensemble  averages  of  system  interactions,  and  the  definition 
and  use  of  the  coupling  loss  factor  are  introduced  here.  The 
use  of  reciprocity  for  development  of  certain  useful 
relations  for  evaluation  of  the  coupling  loss  factor  is  also 
discussed.  The  chapter  ends  with  some  elementary  applications 
of  the  SEA  relations. 
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Chapter  4 completes  this  part  of  the  report  and  is 
concerned  with  the  problems  of  estimating  response,  based  on 
the  average  energy  distribution  calculations.  The  first  part 
of  this  problem  is  estimating  response  variables  of  more 
particular  interest  than  energy  — displacement,  stress,  etc. 
The  second  is  the  development  of  estimation  intervals  and 
their  associated  confidence  coefficients,  particularly  when 
statistical  analysis  of  variance  indicates  that  the 
standard  deviation  is  an  appreciable  fraction  of  the  mean. 
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CHAPTER  2 ENERGY  DESCRIPTION  OF  VIBRATING  SYSTEMS 


2.0  Introduction 


This  report  is  concerned  with  the  exchange  of  energy  in 
coupled  vibrating  systems.  This  chapter  introduces  the  use  of 
energy  variables  in  describing  the  vibration  of  systems. 

Common  ways  of  analyzing  such  systems  employ  modal  oscillator 
or  wave  descriptions  of  the  motion.  A major  goal  of  this 
chapter  is  to  show  how  energy  analysis  applies  to  both  of 
these  models  and  how  certain  relations  between  them  are  re- 
vealed by  the  use  of  energy  variables. 

The  chapter  considers  free  vibration  and  both  sinusoidal 
and  random  excitation  of  vibrating  systems.  This  is  a very 
large  topic,  but  we  are  only  inter' .sted  in  certain  aspects  of 
it.  We  specialize  our  interest  to  linear  systems  and  to  some 
interesting  measures  of  response  that  are  particularly  re- 
lated to  vibrational  energy.  These  measures  are  mean  square 
[ response,  spatial  and  temporal  coherence  and  admittance 

functions. 

The  chapter  begins  with  a discussion  of  the  energetics  of 
modal  resonators  in  free  vibration.  The  cases  of  sinusoidal 
and  random  excitation  are  studied  next.  The  role  of  damping 
in  causing  decay  of  free  vibration  and  bandwidth  in  forced 
vibration  is  of  particular  interest.  Systems  having  many 
modes  of  vibration  are  studied  next,  and  the  notions  of  modal 
density  and  average  admittance  are  introduced.  At  this  point, 
the  very  important  and  pervasive  idea  of  statistical  modeling 
of  real  systems  is  introduced  for  the  first  time. 

Vibrating  systems  modeled  as  collections  of  free  and 
forced  waves  are  introduced  next.  Certain  concepts  unique  to 
such  a description  like  energy  velocity  and  wave  impedance 
are  introduced.  Descriptors  that  are  common  to  wave  and  modal 
descriptions  such  as  mass  and  damping  are  also  discussed.  In 
paragraph  (2.4)  of  this  chapter,  we  point  out  some  relation- 
ships between  modal  quantities,  such  as  modal  r'ensity  and 
their  wave  descriptor  counterparts,  such  as  energy  velocity. 

The  goal  in  this  mode-wave  interplay  is  to  develop  a way 
of  thinking  about  vibrating  systems  that  allows  one  shift 
back  and  forth  between  these  two  viewpoints,  exploiting  the 
one  that  is  best  suited  to  the  problem  at  hand. 
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2.1  Hodal  Resonators 


To  begin  our  discussion,  we  examine  the  energetics  of 
the  simple  linear  resonator  shown  in  Fig.  2.1.  We  shall 
show  later  in  this  chapter  that  this  system  is  a useful  model 
of  the  dynamics  of  the  modal  amplitudes  of  multi-degree-of- 
freedom  (dof)  systems. 

The  dashpot  or  mechanical  resistance  R in  Fig.  2.1 
produces  a force,  — Ry,  opposite  in  direction  to  the 
velocity  y of  the  mass  M.  The  spring  or  stiffness  element 
K produces  a force  — Ky  opposite  in  direction  to  displace- 
ment y from  the  equilibrium  position  of  the  mass.  These 
forces,  in  combination  with  the  force  Mt)  applied  by  an 
external  agent,  results  in  an  acceleration  of  the  mass. 


•• 

Mt)  - Ry  - Ky  « My  , 


or,  more  conventionally, 


y + u)0ny  + y * Mt)/M  (2.1.1) 


where  u)Q  = /k/M  (natural  radian  frequency)  and  n = R/u)QM 
(loss  factor). 

Free  vibration  - No  Damping.  To  study  the  case  of  free 
vibration,  we  set  Mt)  * 6 in  Eq.  (2.1.1).  For  the  moment, 
we  also  neglect  damping  by  setting  n ■ 0.  We  then  have 


y + “Jy  = o#  (2.1.2) 


which  has  the  two  solutions  cos  wQt  and  sin  u)Qt.  The  general 
solution,  therefore,  is  u 


y + A cos  wnt  + B sin  ooftt  = C sin  (wnt+4>) 


(2.1.3) 


are  determined  if  y and  y are  known  at  any  time. 
The  kinetic  energy  of  the  mass  at  any  time  is 


KE  * | My2  = 1 MC2wJ  cos2  (wQt  + (j>)  (2.1.4) 


and  the  potential  energy  in  the  stiffness  element  is 


| PE  » i K y2  = i KC2  sin2(u)Qt  + <J>)  . (2.1.5) 

I 

[ The  sum  of  these  is 


E « KE  + PE  - i KC2 


which  is  time  independent  and  dependent  only  on  the  peak 
amplitude  of  vibration.  Since  we  have  an  isolated  system 
vibrating  without  damping,  it  is  evident  that  its  vibrational 
energy  should  not  vary  in  time. 

The  displacement  and  velocity  repeat  themselves  in  a 
period  l/fQ.  If  we  average  the  kinetic  energy  and  the 
potential  over  such  a period,  we  get 


<KE>^  = <PE>  * | KCJ  « | E (2.1.6) 


Thus,  the  time  average  kinetic  and  potential  energies  are  both 
equal  to  1/2  of  the  energy  of  vibration.  We  shall  make  con- 
siderable use  of  this  relation  in  the  work  that  follows. 
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y(t)  * C e”iwOnt  sin(u)dt  + *).  (2,1.9) 


In  this  case,  the  oscillation  occurs  at  radian  frequency  uj 
and  the  amplitude  of  the  oscillations  decreases  exponentially 
in  time  due  to  the  extraction  of  energy  by  damping.  If  the 
loss  factor  n is  0.5  or  smaller,  then  Wjj  is  very  nearly  equal 
to  uq  and  the  period  of  damped  oscillations  is  essentially  the 
same  as  that  for  undamped  oscillations. 

The  potential  energy  of  vibration  in  this  case  is 


PE  = | K yJ  » j KC2  e a)®ritsin2  (w^t  + 4>)  (2.1.10) 


The  kinetic  energy  is  a little  more  complicated.  It  is 


t. 


KE 


I Ki'.-j  e'“0nt 


“d 

- — cos  (»dt  + «) 

o 


+ i n Sin  (u»dt  + <j») 


(2.1.11) 


These  expressions  simplify  considerably  if  we  average  over  a 
cycle  of  oscillation,  neglecting  the  slight  change  in 
amplitude  in  this  period  due  to  the  exponential  multiplier. 
The  result  is 


<PE>^  * | KC2e“WoT)t  « <KE>^  - | <E>^  (2.1.12) 


In  free  damped  vibration,  for  which  n < 0.5,  we  have  the  same 
relations  between  kinetic,  potential,  and  total  energy  that  we 
obtain  for  undamped  vibration.  Note  that  for  both  damped  and 
undamoed  vibration 


<y2> 

* % 


<y a\  /<*> 


2 

0* 


(2.1.13) 


The  loss  factor  n is  simply  related  to  other  standard 
measures  of  damping.  For  example,  by  its  definition,  n*8 1/Q, 
where  Q is  the  resonator  quality  factor,  much  used  in 
electrical  engineering.  From  Eq.  (2.1.9),  the  amplitude  decays 
as  e-7rri  'T  and,  therefore,7rn  is  the  logarithmic  decrement 
(nepers  per  cycle).  Finally,  the  oscillations  cease  when 
n 4 2 so  that  n is  twice  the  critical  damping  ratio. 

From  Eq.  (2.1.12),  we  have 


<E> 

% 


E0e-wont 


(2.1.14) 


23 


A measure  of  damping  widely  employed  in  acoustics  is  the 
reverberation  time  TR,  which  is  the  time  required  for  the 
vibrational  energy  to  decrease  by  a factor  of  10" *.  Thus, 


e‘uon  tR  * 10“  • 


which  results  in 


(2.1.15) 


Sinusoidal  Forced  Vibration.  If  the  applied  force  i (t) 
in  Fig.  2.1  is  sinusoidal  at  radian  frequency  w,  then  it  is 
convenient  to  use  the  exponential  form 


A(t)  * |L|cos  (wt  +\|>)  * Re{  |L|exp[-iu)t-ii|>]  } (2.1.16) 


where  Re  (...)  means  "real  part  of".  The  dynamical  equations 
are  linear  and  one  can,  therefore,  consider  response  to  the 
actual  excitation  as  a combination  of  (complex)  response  to  the 
complex  excitation  and  its  complex  conjugate.  This  means  that 
we  may  express  the  response  variables  in  terms  of  complex 
exponentials  also. 

Accordingly,  without  loss  of  generality,  we  describe  the 
complex  force  and  velocity  as 


Mt)  + L e 


iwt 


y(t)  = V e 


iojt 


(2.1.17) 
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where  L and  V are  complex  numbers  of  the  form  |L  je”^ 

and  Jv|e“lCl,  respectively.  Differentiation  and  integration 
of  these  functions  is  particularly  simple 


(t)  s iu  L e = iui 


/■ 


i (t)dt  « (-iw)"1!,  e“iwt 


Substituting  for  l and  y in  Eq. 


* A/(-iw)  (2.1.18) 

(2.1.1),  we  obtain 


L * V(-iu)QM)  { (w/Wq-Wq/oj)  + in)  ? VZ  (2.1.19) 


where  Z is  mechanical  i**  .ance  of  the  resonator.  This  can 
also  be  expressed  by  an  admittance  Y * 1/Z, 


V/L  = Y = (w0n  M -i  (wM  - K/u)  r1 


(2.1.20) 


Since  the  magnitudes  of  L and  V are  the  peak  force  and 
”«locity,  we  can  graph  Eq.  (2.1.20)  as  shown  in  Fig.  2.2. 

The, average  power  supplied  to  the  resonator  by  the  source 
is  n *<£y>£.  When  variables  are  described  as  complex  variables, 
the  time  average  of  their  prc  \ct  is  easily  expressed  by  their 
complex  amplitudes,  as  follows. : 


n = <Jiy>t  = Re  (LV*)  = | | L j 2 Re(Y)* 


- | |v j 2 Re  (Z)  . 


(2.1.21) 
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where  ( )*  denotes  the  complex  conjugate.  Since 


Re  (Y)  * REU/Z)  = RE  (Z)/|  Z | 2 = wQn  M |Y|2,  (2.1.22) 


the  input  power  has  the  frequency  dependence  shown  in  Fig.  2.2. 
The  maximum  value  is  11*1/ 2 |l| 2/uipn  M * < £2>t/R,  the  dissipation 
of  a system  with  a resistance  only. 

At  frequencies  w^/k/m  , the  power  diminishes,  reaching 
1/2  its  value  when 


* V0'nM> uV/2,1-231 


which  occurs  when  w * wQ  ± u>Qn/2  (assuming  n < 0.3),  as 
shown  in  Fig.  2.2.  These  frequencies  are  both  the  limits  for 
half-power  and  the  boundaries  for  simple  forms  of  dynamical 
behavior.  When  u < wo(l-n/2),  the  admittance  is  adequately 
represented  by  neglecting  R and  M and  keeping  the  stiff- 
ness term  only 


Y * - i | [«  < u)0(l-n/2)J 


(2.1.24) 


we  denote  this  region  as  "stiffness  controlled". 

For  frequencies  w < cjgd+n/2),  the  admittance  is  usefully 
approximated  by  the  mass  term  only 


Y * i/wM  [w  > U)0(l+n/2)J  (2.1.25) 


The  region  is  called  "mass  controlled".  The  intermediate 
region  is  called  "damping  controlled".  These  regions  are 
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noted  in  Fig.  2.2  also.  These  simplified  limiting  forms  of 
behavior  of  the  resonator  are  of  great  importance  in  our 
consideration  of  systems  with  many  degrees  of  freedom. 

It  is  possible  to  discuss  sinusoidal  excitation  very 
extensively,  and  this  is  done  in  many  textbooks.  For  our 
purposes  here,  however,  we  will  only  note  some  simple 
relations.  First  note  that  the  relative  phase  of  the 
velocity  y with  respect  to  the  force  f is  the  phase  of  the 
admittance  function  Y,  which  is  the  same  as  the  phase  of 
Z*, 


Arg{Z*> 


tan  {u>  - 


w J/w}/w, 


= tan  1 1 2 (a)  - 


V/<V 


(2.1.26) 


It  is  clear  from  this  that  the  phase  is  changing  rapidly  as  the 
system  goes  through  resonance-the  more  quickly  the  smaller  the 
damping.  This  is  in  contrast  to  the  behavior  of  the  amplitude 
which  has  a horizontal  slope  at  w=Wq.  For  this  reason, 
resonance  can  frequently  be  more  accurately  determined  from 
examination  of  the  phase  rather  than  the  amplitude  of  response. 

The  second  point  is  concerned  with  mean  square  response 
at  resonance: 


<y2>  = | |V ! 2 = i |L|  2/u>Jn2M2  (o»=O)0)  (2.1.27) 


Since  y = (V/-iw)e  and  y =-iw  V e”2^,  then  when 


<y2>  = ^ |v|2/wj  = <y2>/wj 

<y2>  * j Wq |V | 2 = wj<y2>,  (2.1.28) 
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the  same  relations  as  found  for  free  vibrations.  These 
relations  do  not  hold  outside  the  damping  controlled  region, 
however,  but  they  do  hold  to  an  acceptable  degree  of  accuracy 

for  Ju)  - w0|  < 1/2  cu^n- 

Random  Excitation.  Most  of  the  applications  of  SEA  that 
concern  us  are  situations  in  which  the  excitation  by  the  load- 
ing environment  is  random.  We  should  emphasize  here,  however, 
that  the  critical  feature  of  SEA  is  that  we  assume  a 
statistical  model  for  the  system  being  excited , not  necessarily 
for  the  excitation.  Thus,  it  is  perfectly  proper  to  apply  SEA 
to  systems  excited  by  pure  tones  if  a statistical  model  of  the 
system  and  the  description  of  its  response  using  energy 
variables  are  appropriate.  Nonetheless,  the  existence  of 
random  excitation  generally  means  that  much  less  averaging 
of  system  parameters  is  necessary  and  consequently  less 
variability  of  response  from  the  calculated  mean  in  any 
particular  situation  may  be  expected. 

There  is  no  single  satisfactory  definition  of  a stationary 
random  signal  from  all  viewpoints,  but  one  that  has  particular 
appeal  from  an  experimental  viewpoint  can  be  readily  developed. 
Imagine  a filter  having  the  frequency  response  shown  in 
Fig.  2.3  and  that  the  load  function  l(t)  is  applied  to  the 
input  of  this  filter.  Now  let  the  bandwidth  Af  of  the  filter 
become  very  small.  If  the  mean  square  output  of  the  filter 
becomes  proportional  to  Af,  then  the  force  &(t)  is  random. 

Note  that  a pure  tone  does  not  satisfy  this  requirement  since 
the  mean  square  output  would  be  independent  of  bandwidth  as 
long  as  the  frequency  of  the  tone  were  in  the  pass  band  of 
the  filter.  A similar  statement  can  be  made  for  any 
deterministic,  periodic  signal. 

This  statement  not  only  supplies  a definition  of  a 
stationary  random  signal  that  fits  out  intuition  and  is 
mathematically  respectable,  it  also  provides  a direct 
indication  of  how  the  frequency  decomposition  of  a random 
signal  is  effected.  The  mean  square  (m.s.)  force  corresponding 
to  the  band  of  frequencies  is,  therefore: 


■ sl  « 


(2.1.29) 


Sj,  being  the  factor  of  proportionality.  It  is  usually  the 
case  that  this  proportionality  factor  depends  on  the  center 
frequency  f,  so  it  is  written  S^(f),  and  it  is  called  the 
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power  spectral  density  (psd)  of  the  random  variable  £(t) 


% 


I 

i. 

i 


Let  us  suppose  that  the  psd  of  £(t)  has  been  determined 
for  all  values  of  f and  the  result  has  been  plotted  as  in 
Fig,  2.4.  We  also  suppose  that  l (t)  is  applied  to  two  filters 
like  the  one  in  Fig.  2.3,  except  they  are  centered  at  fre- 
quencies f,  and  f2.  We  now  combine  the  output  of  the  two 
filters  to  yield 


<*2>  = (f Af  + (f2)  Af  , (2.1.30) 


since  the  mean  squares  of  two  time  functions  containing 
different  frequency  components  are  simply  added  together  to 
form  the  mean  square  of  the  sum.  Proceeding  in  this  way, 
the  mean  square  output  of  a filter  that  has  unity  gain  from 
frequency  f^  to  f2  is  simply 

r2 

<l*>  - / Sp (f)  df  , (2.1.31) 


or,  if  the  filter  has  a gain  G(f)  instead  of  unity,  the 
dependence  is 


<1 2 


S£(f)  G(f)  df. 


(2.1.32) 


The  total  unfiltered  m.s.  value  of  the  loading  function  is 
found  by  setting  f^  * 0 and  f2  = 00  in  Eq.  (2.1.31). 

Let  us  suppose  that  the  loading  force  in  Eq.  (2.1.1)  is 
a noise  excitation  having  a psd  S$,  that  is  constant.  Such  a 
noise  signal  is  called  "white  noise".  The  m.s.  force  pro- 
duced within  a very  narrow  frequency  band  df  is,  therefore: 


<l2> 


Af 


Spdf 


(2.1.33) 


The  m.s.  velocity  response  of  the  resonator  to  this  force  at 
frequency  f is  given  by 


<y2>df  = s£df  | y | 2 


(2.1.34) 


According  to  Eq.  (2.1.32),  therefore,  the  psd  of  y is  given 
(aside  from  a constant)  by  Fig.  2.2,  and  the  total  m.s. 
velocity  is  found  from 


(2.1.35) 


for  white  noise. 

It  is  clear  from  the  form  of  Fig.  2.2  that  most  of  the 
contribution  to  the  integral  in  Eq.  (2.1.35)  will  come  from 
the  region  denoted  "damping  controlled".  This  observation 
leads  to  a useful  concept  of  "equivalent  bandwidth",  A . 
This  is  the  bandwidth  of  a system  with  a rectangular  pals 
band  (as  in  Fig.  2.3)  that  has  a constant  admittance  deter- 
mined by  the  damping  along,  Y@q  = (wpnM)”1,  that  has  the 
same  response  to  the  white  noise  excitation  that  the  actual 
system  does.  Thus,  by  this  definition, 


<y2>  = S^Ae(w0nM)”2 


(2.1.36) 


s„ (w.nM) 


7n  *1^0 


-r 


d w 


l+(w2-w2)  2/n2u >2uJ 


This  integrand  is  simplified  by  changing  variables  to 
C=2  (u>— o)q ) /n<^Q  and  noting  that  the  greatest  contribution  to 
the  integrand  is  at  or  K-0,  Thus, 
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0 —00 


and , therefore , 


Le  * f n f0  ' (2.1.37) 


We  note  that  this  bandwidth  is  greater  than  the  "half  power" 
bandwidth  in  Fig.  2.2  by  the  factor  of  tt/2.  The  "replace- 
ment" of  the  resonator  by  this  equivalent  filter  becomes  a 
very  useful  approximation  in  many  situations. 

There  is  an  alternative  time  domain  definition  of  white 
noise  that  is  useful.  In  this  case,  we  consider  the  force  to 
be  made  up  of  a series  of  impulses  of  strength  ±a,  occurring 
randomly  in  time  as  shown  in  Fig.  2.5.  When  this  excitation 
is  applied  to  the  system  in  Fig.  2.1,  the  mass  is  given  a 
sequence  of  changes  in  velocity.  Let  each  change  in  velocity 
be  Av  « a/M.  These  occur  randomly  in  time  with  random  sign. 
The  corresponding  change  in  energy  of  vibration  is 


AE  * ^M{ (y+Av) 2-y2}*  M y Av  + ^M(Av)2 


(2.1.38) 


If  we  average  this  over  a sequence  of  impulses,  the  term 
<yAv>  will  vanish  since  a positive  Av  is  just  as  likely  as  a 
negative  one.  The  average  energy  increment  is,  therefore, 

1/2  M <Av> 2 = a2/2M. 

If  the  average  rate  of  impulses  per  second  is  v,  then  the 
power  fed  into  the  resonator  is  va2/2M.  But  this  must  also 
equal  the  dissipated  power  <y2>  M.  We  have,  therefore, 

an  estimate  for  the  m.s.  velocity. 


2 

<y2>  = v aV2u>0n  M2  « Sl  Ae  (w0nM)‘  , (2.1.39) 


31 


which  leads  to  SA(f)=2va2,  which  is  a constant.  More  im- 
portantly , however , this  relation  gives  us  an  additional 
physical  interpretation  of  the  spectral  density  function  for 
white  noise. 

If  we  return  to  Eq.  (2.1.36),  the  expression  for  m.s. 
displacement  is  simply 


<y 2 > 


s 

2tt  Sl 


(oi0nM) 


7 


dw 


u)2  (i+(w2-uj)  2/n2w2wJ] 


(2.1.40) 


and  if  the  same  assumptions  are  made  regarding  the  level  of 
damping  are  made  as  previously,  we  get 


<y2>  = <y2>  / co2 


(2.1.41) 


If  we  try  to  calculate  the  m.s.  acceleration  in  the  same  way, 
we  run  into  difficulty,  since  the  integral 


m 2 

<y  * > 


= i-  S 
2 r b 


(wQnM) 


/ 


max 


CO 


dco 


[ 1+  (w2-cj2)  2/n2w2w; 


(2.1.42) 


will  not  converge  as  f -►  «.  For  larne  co,  the  integrand  Is 
simply  wlr}2  • We  can  solve  the  problem  by  "subtracting  out" 
this  part,  and  obtain 


M j 

<y‘>  = 


co; 


<y  > + S, 


WM‘ 


(2.1.43) 
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where  the  first  term  represents  the  damping  controlled  accel- 
eration and  second  part  is  the  mass  controlled  acceleration. 

In  most  instances  of  random  excitation,  the  resonant  part  will 
dominate  and  we  can  simply  use 

<y2>  * o)q2  <y2>  (2.1.44) 

We  have  seen  that  free  vibration  and  the  resonant  control- 
led response  to  sinusoidal  and  random  excitation  all  produce  the 
same  relations  between  m.s.  displacement,  velocity  and  accel- 
eration. This  is  very  useful  to  us  since  it  allows  us  to  trans- 
form from  one  response  variable  to  smother  without  concern  for 
the  exact  nature  of  the  excitation  as  long  as  the  general  require- 
ments of  resonant  dominated  response  for  the  validity  of  these 
relations  are  met. 

2.2  Modal  Analysis  of  Distributed  Systems 

The  systems  of  interest  to  the  mechanical  designer  are 
much  more  complicated  than  a linear  resonator.  In  real 
systems,  the  stiffness,  inertia  and  dissipation  are  all  dis- 
tributed over  the  space  occupied  by  the  structure.  A displace- 
ment of  the  system  that  increases  the  potential  energy  is 
resisted  by  the  elastic  restoring  forces.  A rate  of  change 
in  the  displacement  is  resisted  by  the  damping  forces,  and 
these  forces,  along  with  the  loading  excitation  cause  acceler- 
ation of  the  mass  elements. 

If  we  represent  the  generalized  displacement  of  the 
system  oy  v,  then  an  equivalent  statement  of  the  above  is 

py  + ry  + Ay  = p,  (2.2.1) 

where  p is  the  mass  density,  r is  a viscous  resistance 
coefficient  and  A is  a linear  operator  consisting  of 
differentiations  with  respect  to  space.  In  the  case  of  a 
flat  plate,  for  example,  o is  the  mass  per  unit  area  of  the 
plate  and  A = B 74 , where  B is  the  bending  rigidity.  The 
use  of  simple  viscous  damping  is  a valuable  simplification  and 
does  not  affect  the  utility  of  our  results  as  long  as  the 
damping  is  fairly  small. 

When  this  system  is  bounded,  with  well  defined  boundary 
conditions,  the  solution  is  frequently  sought  by  expansion  in 
eigenfunctions  , which  are  solutions  to  the  equation 

“n  *n  (2-2'2) 
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where  the  functions  tj*n  satisfy  the  same  boundary  conditions 
that  y does,  and  the  quantities  are  determined  by  the 
boundary  conditions.  The  response  and  the  excitation  are 
then  expanded  in  these  functions 


y=  J Yn(t)  *n(x) 


p/p  * 1 Lit)  ip  (x)  (2.2.3) 

n n n 


If  we  multiply  Eq.  (2.2.2)  by  ip  and  integrate  wer  the 
region  of  the  structure  and  subtract™  from  this  the  same 
equation  with  the  indices  reversed,  we  get 


{V*n 


- W 


dx  * 


(wj-wj) 
n m 


/ 


i>  p (x)  ill 
vnr  ' ' yn 


dx. 


(2.2.4) 


When  n * m,  this  is  satisfied  in  a trivial  fashion.  When 
n^n,  there  are  certain  specific  (but  nevertheless  very  use- 
ful) conditions  under  which  the  differential  operator  has 
the  property  that 


/ <V'*n  - W ** 
will  vanish.  This  means  that 

/V  *n  dx 
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must  vanish,  a condition  that  is  referred  to  as  orthogonality 
of  the  eigenfunctions  with  a weighting  function  p. 

A convenient  normalization  of  the  amplitude  of  the  eigen- 
functions is  simply 


\l>  p ill 

nr  rn 


dx 


= <ib  > 

Mtrn 


*m,n 


(2.2.5) 


which  we  may  think  of  as  a mass  density  weighted  average  of 
the  product  . For  systems  in  which  the  mass  density 

is  uniform,  this  becomes  a simple  average  over  the  spatial 
coordinates. 

If  we  now  place  the  expansions  of  Eq.  (2.2.3)  into 
Eq.  i 2. 2.1),  we  obtain 


l 

n 


a 


n 


r • 

+ £ y 

P n 


n n 


u>fY)  * 


n 


p I L \lf 

M „ m Mn 
m 


(2.2.6) 


We  can  simplify  this  immediately  if  we  also  assume  that  r(x) 
is  proportional  to  p(x);  r = Ap.  We  can  do  this  on  two 
counts.  To  take  a somewhat  cynical  view,  since  we  have 
introduced  the  damping  in  a rather  ad  hoc  fashion,  we  can 
feel  free  to  configure  it  any  way  we  like.  The  other  basis 
for  the  assumption  is  that  research  studies  of  this  problem 
have  shown  that  the  consequence  of  this  assumption  is  to 
ignore  a degree  of  inter-modal  coupling  that  is  less 
significant  than  other  forms  of  coupling  that  we  will  be 
concerned  about. 


We  now  multiply  Eq.  (2.2.6)  by  (x)  and  integrate  over 
the  system  domain.  Using  Eq.  (2.2.5),  we  obtain 


M{Y  +AY  +w  2Y}  = L (t) 
m m m m m 


(2.2.7) 
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Thus,  each  modal  response  amplitude  obeys  the  equation  of  a 
linear  resonator  of  the  sort  discussed  in  paragraph  2.1.  This 
result,  in  conjunction  with  the  spatial  orthogonality  of  the 
mode  shapes  according  to  Eq.  (2.2.5),  leads  us  to  the  concept 
of  a complex  dynamical  system  as  a group  of  independent 
resonators  of  mass  M,  stiffness  w2M  and  mechanical  resistance 
MA.  m 


Since  there  is  a one-to-one  correspondence  between  the 
modes  and  the  wn's,  we  can  use  the  latter  to  keep  track  of 
the  modes.  An  example  should  be  helpful  here.  For  the  two- 
dimensional,  simply  supported,  isotropic  and  homogeneous 
rectangular  flat  plate  of  dimensions  A,  by  A,,  drawn  in 
Fig.  2.6,  the  mode  shapes  are  1 1 


2 sin 


sin 


n2irx2 


(2.2.8) 


and  w2  is 
n 


k2c2 


2 

= (k£+k2)K2c£  = kV2c2  (2.2.9) 


where  k is  the  radius  of  gyration  of  the  plate  cross- 
section,  C£  = /Yp/pm  is  the  longitudinal  wavespeed  in  the 
plate  material,  where  Y is  the  plate  Young's  modulus  and 
Pm  is  the  material  density  and  n^  and  n2  are  integers. 

Eq.  (2.2.9)  suggests  a very  convenient  ordering  of  the 
modes,  which  is  shown  in  Fig.  2.7.  By  inspection,  we  can 
see  that  each  point  in  this  "wave-number  lattice"  corresponds 
to  a mode.  Further,  the  distance  from  the  origin  to  that 
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point  will  determine  the  value  the  resonance  frequency  of 
the  mode  wn.  The  main  value  of  this  ordering  is  that  it 
enables  us  for  example  to  count  the  modes  that  will  resonate 
in  some  frequency  interval.  It  will  not  always  be  as  con- 
venient as  it  is  in  this  case,  but  the  ordering  by  some 
parameter  is  a necessary  condition  to  allow  us  to  do  the 
counting  at  all. 

When  the  ordering  indices  form  a lattice  as  shown  in 
Fig.  2.7,  then  each  lattice  point  corresponds  to  an  area 
AAk^irVAp,  where  A Xi“!Llli2  is  the  area  of  the  Plate*  As  we 
increase  the  wavenumber  from  k to  k + Ak,  we  include  a new 
area  1/2  trkAk.  On  the  average,  this  will  include  1/2  irkAk/AAfc 
new  modes.  Thus,  the  average  number  of  modes  per  unit 
increment  of  wavenumber  is 


n(k) 


TrkAk 

2AAkAk 


(2.2.10) 


which  we  may  call  the  modal  density  in  wavenumber. 

To  find  the  modal  density  in  frequency  n(w),  which  tells 
us  how  many  modes  on  the  average  are  encountered  when  we 
increase  frequency  by  1 unit,  we  use  the  relation 
n(w)Aw  * n(k)  A k,  or 


n (w) 


Ak  _ irk 
Aw  = 2cgAAk 


(2.2.11) 


where  we  have  used  the  result  from  elementary  physics  that 
Aw/Ak  is  the  group  velocity  for  waves  in  a system  that  has  a 
phase  velocity  w/k. 

For  a flat  plate,  the  group  velocity  c is  twice  the 
phase  velocity  c^  * w/k  * /wkc^  , so  that  the  modal  density 
in  cycles  per  second  (hertz)  is 


n(f)*n(w) 


dw 

d I 


2TTawA 

teg? 


(2  2.12) 
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where  we  have  used  k = h/2  /3 , the  radius  of  gyration  for  a 
homogeneous  plate  of  thickness  h.  Note  that  the  modal 
density  in  this  case  is  independent  of  frequency.  As  an 
example,  consider  a plate  that  has  an  area  of  10  ft.2  and  a 
thickness  of  1/8  in.  (approximately  10-2  ft.).  Then,  since 
c^  = 17,000  ft. /sec.  (for  steel  or  aluminum),  we  have 


n(f)  * (,oI)  U7,000)  ~ 0,1  mode/Hz  (2.2.13) 


For  this  plate,  one  mode  is  encountered  on  average  whenever 
the  excitation  frequency  is  increased  by  10  Hz.  A plate  with 
a larger  area  or  smaller  thickness  will  have  a modal  density 
that  is  greater  than  this. 

The  kinetic  energy  of  vibration  is 


5 / dx  p(lt)  ’ l/v  Sni.ltlin,tlM»W*,W 


- i M I Y2 (t)  , (2.2.14) 

i _ m 
m 


where  we  have  used  the  orthogonality  relation.  Thus,  the 
kinetic  energies  of  the  modes  add  separately  to  give  the  KE 
of  the  system.  Since  the  kinetic  and  potential  energies  of 
each  resonator  are  equal  at  resonance 


(<^n>  * “n  ' 


the  total  energies  of  the  modes  simply  add  to  form  the  total 
system  energy. 
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Response  of  System  to  Point  Force  Excitation.  We  now 
imagine  that  a point  force  of  amplitude  LQ  is  applied  at  a 
location  xg  on  the  structure.  If  we  solve  for  the  modal 
amplitudes  from  Eq.  (2.2.3),  we  get 


= 5 /p* m dlt  * L0  *m(xs,/M  ' 


(2.2.15) 


where  M is  the  system  mass.  If  the  excitation  and  response 
are  proportional  to  e 1U)t,  Eq.  (2.2.7)  becomes 


Mtw£~iw  V "w2)  Wn(xs)  ' 


(2.2.16) 


so  that  the  formal  expression  for  the  response  is 


y(x,t) 


L e'1"* 
o 

M 


WV10 

wl-w2-i  u wn 
n n 


(2.2.17) 


We  shall  look  at  some  ways  of  simplifying  this  complicated  re- 
sult. 

The  velocity  at  the  excitation  point  xs  is  -ituy(xs).  The 
ratio  of  this  velocity  to  the  applied  force  is  the  input  con- 
ductance of  the  system. 


-iwY  (x_ ,w) 
s 


*n(xs> 
n n 


5 G - iB  . 


(2.2.18) 
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g 


where  G,  the  real  part  of  the  sum  is  the  conductance  and 
B,  the  imaginary  part,  is  the  susceptance.  By  rationalizing 
the  complex  in  Eq.  (2.2.18),  we  get 


G 


l an(xs)gn(w) 


(2.2.19) 


where  an  * i(i£(xs)A>Mn  and  gn*(C2+l)”1,  where  5*(u>n-a))/wn.  The 
susceptance  is 


B * l a (x  )b  (w) 
n n s n 


(2.2.20) 


where  bn*£ (C2+l)-1 • In  deriving  these  relations,  we  have 
assumed  that  the  damping  is  small  enough  so  that  the  individual 
modal  admittance  an(9n~ibn)  is  quite  sharp  in  frequency. 

The  total  admittance  Y is  a rapidly  fluctuating  function 
of  frequency.  We  can  simplify  this  result  however,  if  we  con- 
sider averages  of  Y with  respect  to  the  variable  £.  Such 
an  average  is  appropriate  if  we  are  exciting  the  system  with 
a band  of  noise,  since  the  response  to  noise  having  a uniform 
spectral  density  from  to  u>2  is  the  same  as  the  average  of 
the  m.s.  response  to  a pure  tone  as  may  be  seen  by  referring 
to  Eqs.  (2.1.32)  and  (2.1.35). 

Another  possibility  is  to  assume  that  the  system  itself 
is  "random".  That  is,  that  the  exact  mode  shapes  and  resonance 
frequencies  are  not  known,  either  because  of  random  irregularities 
in  their  construction,  or  because  the  detailed  calculation 
procedures  are  not  accurate  enough  to  calculate  them.  In  this 
case,  we  assume  that  the  resonance  frequencies  u>n  are  uni- 
formly distributed  over  some  frequency  interval.  Such  an 
approach  is  one  ' ample  of  statistical  modeling,  a central 
theme  in  Static -real  Energy  Analysis.  In  the  work  that 
follows  we  take  this  latter  approach,  assuming  that  w is 
fixed  and  that  is  the  random  variable. 
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Aw, 


If  the  interval  of  resonance  frequency  uncertainty  is 
then  if  Aw  » j wn#  the  average  over  wn  will  give 


<g  > 


w 


n 


wn 

2aw 


-/’?%■  SB 


(2.2.21) 


which,  interestingly  enough,  is  simply  the  ratio  of  the  modal 
bandwidth  to  the  averaging  bandwidth.  In  the  interval  Aw, 
the  number  of  modes  that  can  contribute  to  the  average  is 
n (w) Aw.  Thus,  from  Eq.  (2.2.19),  the  average  conductance  is 


<G> 


Vys 


" un  .I”M 

2Aw  r'\  2 M 


(2.2.22) 


where  the  average  on  ip 2 is  over  the  mass  distribution  of  the 
system,  which  is  of  course  the  same  as  a spatial  average  for 
a uniform  mass  density. 

Eq.  (2.2.22)  is  a useful  and  general  result  for  multi- 
modal systems.  In  the  case  of  a flat  plate,  n(w)=A  /4tt  kc^, 

and  M=p  A , where  p_  is  the  surface  density  of  the  plate.  Then, 
SO  s 


<G>  = (8p  Kcj"1  (2.2.23) 

S A/ 


which  is  also  the  admittance  of  an  infinite  plate.  It  very 
often  happens  that  average  impedance  functions  of  finite 
systems  are  the  same  as  those  for  the  same  system  infinitely 
extended. 

The  average  susceptance  is  of  less  interest,  but  we  note 
that  when  the  modal  density  is  constant,  then  the  average  sus 
ceptance  will  vanish  because  the  integral  of  b„(£)  vanishes. 
When  the  modal  density  is  not  constant,  the  calculation  is 
more  complicated,  and  has  been  dealt  with  in  the  references. 
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Before  leaving  the  discussion  of  impedance,  let  us 
consider  the  same  problem  from  the  point  of  view  of  noise 
excitation.  If  the  force  in  Eq.  (2.2.15)  has  a flat  spectrum 
S&  over  the  band  Aw,  then  the  power  fed .into  any  one  mode 
can  be  found  from  the  dissipation  w0nM<y2>is,  according  to 
Eq.  (2.1.39),  just  S^n  (x_)/4M.  If  the  bandwidth  of  the  noise 
is  Aw  (radians/sec)  or  Aw/2ir  Hz  then  the  number  of  modes 
randomly  excited  is  nAw  and  the  power  input  to  the  system, 
averaged  over  source  location  is 


<n> 


ir 

1 


SjAw 

2tT 


* <f2>  <G> 


(2.2.24) 


where  again  <G>  » im/2M.  Thus,  we  can  view  the  conductance  as 
a measure  of  the  number  of  modes  that  are  available  to  absorb 
energy  from  a noise  source.  On  this  basis,  it  is  apparent 
that  there  should  be  a close  tie  between  the  average  con- 
ductance and  modal  density.  This  relation  can  in  fact  be 
exploited  as  a way  of  measuring  modal  density  when  the  modes 
are  so  closely  packed  or  the  damping  is  so  large  that  counting 
resonance  peaks  is  not  feasible. 


2.3  Dynamics  of  Infinite  Systems 


When  the  system  is  infinitely  extended,  then  an 
alternative  formulation  is  needed.  The  differential  equation 
governing  the  motion  is  still  given  by  Eq.  (2.2.1).  We  now 
assume,  however,  that  the  mass  and  damping  distributions  are 
uniform  (p  and  r * constant)  and  that  the  linear  differential 
operator  A is  a simple  polynomial  in  the  spatial  derivatives, 
A(3/dX£)  with  constant  coefficients.  With  these  assumptions, 
we  may  assume  a "wave"  solution  to  the  equation  for  unforced 
motion  in  the  form 


i()c.x  -wt) 

y - e 


(2.3.1) 
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With  this  substitution , the  equation  of  motion  becomes  the 
"dispersion  relation"  between  frequency  and  wave  number: 


-pw2  - iwr  + A(ik^)  = 0 


(2.3.2) 


Let  us  consider  some  simple  examples  of  Eq.  (2.3.2).  For 
the  undamped  string,  p = lineal  density,  r * 0,  and  A=-T(3/3x)2. 
Thus,  Eq.  (2.3.2)  becomes 


Tk2  = pw2,  or  k = ±w/c  , 


(2.3.3) 


where  c=/T/p  is  the  speed  of  free  waves  on  the  string.  In 
the  case  of  undamped  bending  motions  of  a thin  beam,  we  would 
have 


A (3/3x)  = B (3/3x)4 


(2.3.4) 


where  B = p c2<2  is  the  bending  rigidity  of  the  beam.  In  this 
case  the  dispersion  relation  is 


= ^2/<2c2  - (lo/C^)4 


(2.3.5) 


the  same  as  Eq.  (2.2.9).  The  parameter  Cw  is  the  phase 
velocity  for  bending  waves  on  the  beam.^  In  the  case  of  a 
two  dimensional  plate,  the  wave  vector  k in  Eq.  (2.3.1)  has 
the  two  components  (k^,k2).  For  either  the  beam  or  the 
plate,  the  phase  velocity  cb  is  a function  of  frequency 


cb  = /uj  * 


(2.3.6) 
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and  the  system  is  said  to  be  dispersive.  When  damping  is 
included,  the  propagation  constant  k is  complex  and  an 
attenuated  wave  results. 

Energy  variables  are  of  great  importance  in  infinite, 
free  wave  systems  just  as  they  are  for  finite,  modal 
descriptions.  In  the  infinite  system,  we  are  interested  in 
energy  density,  the  energy  of  vibration  per  unit  length, 
area,  or  volume  depending  on  the  dimensionality  of  the 
system.  It  is  shown  in  advanced  textbooks  that  for  the 
system  that  we  are  considering,  the  kinetic  and  potential 
energy  densities  of  free  waves  are  equal,  so  that  the  total 
energy  density  is  just  twice  the  kinetic  energy  density, 
which  is  p(9y/3t)2. 

The  intensity  I of  a free  wave  is  equal  to  the  power 
flowing  through  a unit  width  (or  area)  of  the  system  due  to 
that  wave  as  it  propagates.  If  this  power  flows  for  1 second, 
then  the  amount  of  structure  that  has  filled  with  energy  is 
numerically  equal  to  c„  = dw/dk,  the  energy  velocity.  If  C 
is  the  energy  density,  then  the  energy  that  passed  the  reference 
location  is  ffc=I.  From  the  dispersion  relation,  Eq.  (2.3.2) 
with  r=0,  the  energy  velocity  is 


c - Jr  - 777-  A*  (ik)  . (2.3.7) 

g dk  2poj 

Since  E * P (3v/3t) 2 , 

I - - y'  <y2>t  A'  (ik)  . (2.3.8) 


Let  us  consider  some  examples  of  the  use  of  Eq.  (2.3.8) 
for  some  familiar  systems.  In  the  case  of  the  string, 

A(ik)  = - T(ik)2  and  A' (ik)  * - 2T(ik).  Thus, 


String"  ’ F^'t  (’pc2)  2ik  = pc  <(9y/3t)2>t 

(2.3.9) 
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In  the  case  of  the  beam,  A *p  <2c|(ik)1'  and  A'  = 4p  x2c|(ik)J. 
Thus 


beam  2u 


(nV> 

\3ty  t 


4p  <2c|(ik) 


2pcb  < 


(l?)2>t 


(2.3.10) 


In  both  of  these  cases,  the  intensity  is  a mean  square  velocity 
of  motion  of  the  system  times  an  impedance  term  of  the  form 
pc,  where  p is  the  density  and  c is  a wavespeed. 

When  damping  is  included,  its  most  important  effect  is 
to  cause  the  propagation  constant  to  become  complex.  The 
new  dispersion  relation  is  found  by  the  substitution. 


k (cu) 


(2.3.11) 


Thus,  for  the  string, 


ik  x i-[l+(iV2Mx  i~  - *3-  x (2.3.12a) 

e s * e c *e  e (string) 


and  for  flexural  motion  of  the  beam  (k^  = oj/c^) 


ik,  x 
o 


4 


uj  (1+  (in/2) ) 
<cn 


(beam) 


(2.3.12b) 
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From  this,  it  is  clear  that  the  form  of  the  dispersion  relation 
will  affect  the  rate  at  which  the  wave  decays  in  space  as 
damping  is  added.  For  a non-dispersive  system,  the  attenuation 
is  2 nr)  nepers  or  27.3  n dB  per  wavelength.  It  is  clear  that 
if  we  used  the  energy  velocity,  the  attenuation  factors  in 
Eqs . (2. 3.12a)  and  (2.3.12b)  would  be  formally  the  same  since 
Cg  = 2 C5  for  the  plate. 

Most  of  the  structures  of  interest  to  structural  designers 
consist  of  segments  of  beams  and  plates,  so  that  we  shall  place 
most  of  our  emphasis  on  such  systems.  In  paragraph  2.2,  we 
studied  the  impedance  looking  into  a finite  plate,  and  found  that 
the  average  impedance  was  real  and  related  to  the  modal  density 
and  the  mass  in  a particularly  simple  fashion.  We  now  want  to 
carry  out  this  same  calculation  for  the  infinite  plate. 

The  two  dimensional  flat  plate  has  the  equation  of  motion 
Eq.  (2.2.1)  with 


The  point  force  of  strength  L«e~  w is  assumed  to  act  at  x * 0. 
The  problem  is  to  calculate  tne  motion  y at  x ■ 0 and  form  the 
ratio  of  [dy/dt)  to  the  force  to  find  the  admittance  of  the 
plate.  We  solve  the  problem  by  using  two-dimensional  Fourier 
transforms 


(2.3.13) 


with  similar  relations  between  y(x)and  its  transform  Y (k) . 

Since  p(x)  acts  only  at  x = 0 with  strength  L«,  the  second 
integral  above  is  simple  P(k)  = Lq.  Placing  the  transform  into 
the  equation  of  motion 
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(2.3.14) 


[-03 2 (l+in)P  + A(ik1#ik2)J  Y()c)  = LQ 


which  gives 


Y (0) 


± (2  1 K\ 

A (ik)  - 03 2 (l+in)P 


Since  there  is  no  azimuthal  dependence  in  the  (k.,k2)  integration, 
we  replace  the  area  element  by  dk^  dk2  -*>2irk  dk  and  have  a simple 
integration  over  the  magnitude  k. 


y (0) 


k dk  * 

ptc2cj[k4-k£(l+in)] 


--iS_  f"  dc i 

4ttpk2c2  J 0 £2-££(l+in) 


where  £ * k2  and  £b  = k£.  Thus, 


y (0) 


8np*.  2c2 


X 


d£ 

{ C-eb  ( 1+i n75 ) ] f 5+5b  fl+Tn/27 ] 


(2.3.16) 
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y (0) 


or. 


L0/(-iu» 


8pxc£ 


allowing  n to  vanish/  and  the  path  of  integration  is  taken 
as  in  Fig.  2.8.  Thus,  the  ratio  of  velocity  -iw  y(0)  to  the 
force  Fq  is  the  infinite  plate  admittance 


* (8  pkCjT1  (2.3.17) 


which  is  the  impedance  that  we  found  for  the  finite  plate  when 
an  average  was  taken  over  modal  response  and  source  location. 

A mean  square  force  <12>  applied  at  a point  on  an 
infinite  plate  will,  therefore,  inject  an  amount  of  power 


11  in  * <*’>  (2.3.18) 


into  the  plate.  This  power  will  propagate  at  the  energy 
velocity  outward  from  the  point  of  excitation  as  a circular 
wave.  When  a boundary  is  encountered,  a reflection  occurs  and 
the  energy  propagates  unimpeded  until  another  boundary  is 
encountered.  The  average  distance  traversed  between 
reflections  is  called  the  "mean  free  path",  d,  and  is  given 
by  d = TrAp/P,  where  P is  the  perimeter  of  the  plate  and  Ap 
is  the  area  of  the  plate.  The  attenuation  rate  (in  space) , 
according  to  Eq.  (2.3.12b)  is  wn/2Cg  nepers  per  unit  length, 
or  4.3  4o)ndB  per  second  (c_  is  the  energy  velocity  defined  by 
Eq.  (2.3.7)).  g 

For  systems  in  which  the  energy  is  contained  in  the 
propagation  of  free  waves,  the  damping  is  frequently  expressed 
in  terms  of  the  time  (referred  to  as  the  reverberation  time  Tr) 
required  for  the  energy  to  decay  by  60  dB.  Thus,  4.34  air;  TR  = 60 
or, 
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(2.3.19) 


T 


R 


2.2 

Fn' 


a result  that  we  found  earlier  for  single  dof  and  modal 
resonators.  Measurement  of  decay  rate  or  reverberation  time 
is  a commonly  useo  procedure  for  determining  loss  factor. 

The  decay  rate  procedure  applies  equally  well  whether  our 
model  of  the  system  is  one  of  a collection  of  modes  or  a group 
of  waves  rebounding  within  the  system  boundaries. 

In  Chapter  3#  we  will  be  concerned  with  the  dynamics  of 
interacting  systems.  As  a preview,  1st  us  consider  here  the 
interaction  between  the  single  dof  system  and  a finite  plate 
shown  in  Fig.  2.9.  We  wish  to  calculate  the  mean  square 
velocity  of  the  resonator  as  a result  of  its  attachment  to 
a plate  which  is  vibrating  randomly. 


The  resonator  consists  of  a mass  M,  stiffness  K,  and 
dashpot  resistance  R,  configured  as  shown.  A diffuse 
reverberant  vibrational  field  (equal  wave  intensity  in  all 
directions)  is  assumed  to  exist  on  the  plate,  resulting  in  a 
transverse  velocity  v.  At  the  point  of  resonator  attach- 
ment xg,  the  plate  velocity  is  v_.  The  velocity  of  the 
resonator  mass  is  vM.  The  reaction  force  l resulting  from 
compression  of  the  spring  K is  due  to  a difference  in 
velocities  ve  and  vu: 

5 M 


*/ 


(VS-VM> 


dt 


M 


dVM 

~3t 


+ R v 


M 


(2.3.20) 


The  velocity  vg  is  equal  to  the  velocity  that  would  exist 
if  there  were  no  resonator  v,  less  the  "reaction"  or  "induced" 
velocity  vr  that  is  induced  by  the  force;  vr  = 1<G> , where 
<G>  is  the  plate  admittance  (8p~  KCf,)-1.  In  using  <G> , we 
assume  that  the  modal  density  of  the  plate  np  is  sufficiently 
high  so  that  several  plate  modes  resonate  within  the  combined 
equivalent  bandwidth  of  plate  and  resonator  modes: 


5 u^p  + nQ)  np  » 1. 

We  now  differentiate  Eq.  (2.3.20)  with  respect  to  time  to  obtain 


49 


d2v  dvM 

M + R + K vM=  K v - K Jt  <G>  (2.3.21) 

at2  dt  M 


Substituting  for 


l from  Eq.  (2.3.20),  we  get 


d2vM  dvM 

^7"  + J*0(n0+ncoup)  “dt  + w0 {1+n0ncoup)vM*w0V  (2,3,22) 


where  uQ  = /k7m,  rig  = R/wQM  and  ncoup  is  the  combination  of 
of  parameters  wflM  <G>.  Eq.  (2.3.22;  is  the  equation  for  the 

response  of  a resonator  to  a random  base  excitation  velocity 
v.  The  effects  of  the  plate  are  expressed  by  modified 
damping  and  resonance  frequency  of  the  resonator.  If  v is 
assumed  to  have  a flat  spectrum  <v2>/Aw  over  a bandwidth  Aw, 
then  from  the  discussion  following  Eq.  (2.1.35),  the  response 
of  the  resonator  is 


"V  ■ 
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<V2> 

Aw 


(2.3.23)  j 


or, 


M <VM>= 


n 

coup 

n0  + ^coup 


M < v 2 > 

_E 

n Aw 
P 


I 

(2.3.24)  ! 


The  term  Mp  <v*>/npAw  is  simply  the  vibrational  energy  of 

the  plate  divided  by  the  number  of  modes  that  resonate  in  the 
band  Aw.  The  ratio  represents,  therefore,  the  average  energy 


per  mode  of  vibration  within  that  band.  Since  the  ratio 
Hcoup/fncoup  + no)  is  always  less  than  or  equal  to  one, 

Eq.  12.3.24)  says  that  the  average  energy  of  vibration  of  the 
resonator  cannot  exceed  the  average  modal  (resonator)  energy 
of  the  plate.  These  resonator  energies  will  in  fact  be  equal 
when  the  loss  factor  due  to  the  coupling  between  the  resonator 
and  the  plate  Hcoup  (also  called  the  coupling  loss  factor)  is 
large  compared  to  the  internal  damping  of  the  resonator.  This 
tendency  toward  energy  equality  is  an  example  of  "equipartition 
of  energy",  a principle  that  is  well  known  to  students  of 
statistical  mechanics. 


2.4  Modal-Wave  Duality 


In  paragraphs  2.2  and  2.3  of  this  chapter,  we  have  dis- 
cussed modal  and  wave  descriptions  of  finite  systems.  These 
two  ways  of  describing  the  motion  are  in  most  cases  fully 
equivalent  to  each  other,  but  that  does  not  mean  that  the 
choice  of  description  is  arbitrary.  Certain  aspects  of 
structural  vibration  are  much  more  readily  interpreted  in  one 
description  than  in  the  other.  The  effect  of  damping  at  the 
boundaries  of  a plate  is  better  described  by  wave  reflection 
processes.  The  spatial  variations  in  vibration  amplitude  are 
better  described  by  the  modal  analysis.  But  we  must  emphasize, 
that  it  is  always  possible,  at  least  in  principle,  to  arrive 
at  the  same  conclusions  by  either  approach. 

We  have  already  demonstrated  an  important  example  of 
modal-wave  duality,  the  average  point  conductance  of  finite 
plates.  We  have  found  that  the  conductance  computed  by 
averaging  over  modes  of  a finite  plate  is  the  same  as  that 
found  by  considering  only  waves  radiated  outward  by  driving 
an  infinite  plate  at  one  point.  The  fact  that  these  waves 
reflect  from  the  boundaries  and  may  return  to  the  drive  point 
is  presumed  to  be  unimportant  in  affecting  the  drive  point 
impedance  because  they  will  have  random  phases,  or  be  in- 
coherent with  the  excitation,  particularly  if  the  drive  point 
is  assumed  to  be  randomly  located  over  the  surface  of  the 
structure. 

Another  aspect  of  this  duality  is  in  the  description  of 
damping.  We  found  in  discussing  resonators  in  paragraph  (2.1) 
that  the  damping  was  closely  associated  with  resonator  (or 
modal)  bandwidth  and  decay  of  vibrations  in  time.  In  dis- 
cussing waves,  in  paragraph  (2.2)  we  found  that  the  damping 
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was  related  to  the  spatial  decay  of  energy.  However,  if  this 
decay  in  space  were  related  to  the  tine  required  for  the 
energy  to  propagate  from  one  point  to  the  other,  the  tine  decays 
in  both  descriptions  were  related  to  the  loss  factor  in  the 
same  way.  Thus,  decay  rate  of  energy  forms  a link  between  the 
two  descriptions  even  though  modal  bandwidth  is  a very  dif- 
ficult concept  to  explain  by  a wave  analysis  and  spatial  decay 
of  vibration  is  an  equally  difficult  concept  to  explain  using 
a modal  description. 

In  the  chapters  that  follow,  particularly  in  examples 
dealing  with  applications  of  SEA,  we  find  that  certain 
quantities  that  enter  the  modal  description  for  example,  are 
equivalent  to  other  quantities  in  the  wave  description.  For 
example,  modal  energy  is  usually  equivalent  to  a spatial 
energy  density.  For  example,  in  a plate,  the  energy  is  ps<v2>, 
whereas  the  average  modal  energy  for  the  same  plate  is 
(4tt  kc£)  pg  <v2>.  In  a sound  field,  the  energy  density  is 
<P2>/pc2,  and  the  average  modal  energy  for  the  same  sound 
field  is  (2ir2c3/w2) <p2>  p c2,  where  c is  the  speed  of  sound 
and  p is  the  pressure  fluctuation. 

Other  equivalences  exist  between  coupling  loss  factors, 
appropriate  for  modal  systems,  and  junction  impedances, 
appropriate  for  wave  descriptions.  An  example  of  this 
equivalence  was  used  in  the  example  of  the  combined  plate- 
resonator  systems  of  Fig.  2.9,  where  we  found  r)c  woM  <G>* 

We  shall  encounter  more  examples  of  this  in  the  chapters  to 
come. 

In  the  remainder  of  this  section,  we  shall  discuss  a very 
important  aspect  of  wave-mode  duality.  This  is  the  problem  of 
the  "coherent"  and  "incoherent"  wave  fields  in  a plate.  In 
paragraph  2.3  we  began  this  discussion,  but  we  wish  to  delve 
further  into  the  matter  at  this  point.  Let  us  begin  by  con- 
sidering the  wave  description  first. 

Let  us  suppose  that  a point  force  located  at  position  xg 
on  a very  large  plate  excites  the  plate  in  a bandwidth  Ato. 

If  the  rms  force  is  L,  then  the  power  injected  into  the  plate 
at  xg  is  simply 


nin  * l2y»  " L‘<G>  (2.4.1) 

where  <G>  is  given  by  Eq.  (2.3.17).  This  energy  flows 
radially  outward  from  this  point,  so  that  at  a distance  r 


f 
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from  x , the  mean  square  velocity  is  found  from  the  intensity 
relation 


P cr 


<v*> 


2irr  * n 


in* 


(2.4.2) 


This  is  the  velocity  that  would  exist  on  an  infinite  plate 
without  damping.  If  damping  is  present#  the  energy  is  re** 
duced  by  the  factor  exp  I-nwr/c  1 # so  that 


<vf>« 


Jln  i -*">r/cg 
2ir  p c„  r 9 


(2.4.3) 


This  part  of  the  plate  motion  is  called  the  "direct  field" 
of  the  source  and  has  a "geometric  attenuation"  of  3 dB  per 
double  distance  and  an  attenuation  due  to  damping  that  in- 
creases linearly  with  r.  The  direct  field  is  the  only  con- 
tribution to  plate  vibration  if  the  plate  is  infinite  in 
extent  or  if  n is  so  large  that  the  vibration  has  nearly 
ceased  by  the  time  the  direct  wave  reaches  the  boundary  of 
the  plate. 

When  the  energy  in  the  direct  field  reflects  from  the 
boundary#  then  if  the  boundary  is  not  perfectly  regular#  it 
is  frequently  assumed  that  coherence  with  the  direct  field 
is  lost#  particularly  as  the  number  of  reflections  mounts. 

The  motion  vR  corresponding  to  this  reflected  energy  is 
called  the  "reverberant"  field  and  is  determined  by  the  power 
dissipated  by  it: 


nR  * nin  exp  I-wnd/2CgJ  * M <vR>  wn  (2.4.4) 


where  d is  the  mean  free  path.  This  leads  to  an  expression 
for  the  reverberant  velocity 


<v|>  * n.  exp  l~wnd/2c  ]/u>r)M. 
it  in  ^ 


(2.4.5) 
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Since  the  direct  and  reverberant  velocity  components  are 
assumed  to  be  incoherent , the  total  mean  square  velocity  of 
vibration  is  obtained  by  simple  addition, 


<v2>  « <vJ?>  + <v2>. 

D K 


(2.4.6) 


Near  the  point  of  excitation,  the  direct  field  will  dominate, 
but  at  large  distances,  the  reverberant  field  will  dominate. 
The  "boundary"  between  the  direct  and  reverberant  fields  is 
defined  as  the  distance  rD  for  which  these  fields  are  equal 


s j£n  M 
rD  2ir  p cg 


(2.4.7) 


The  total  mean  square  velocity  pattern  is  shown  in  Fig.  2.10. 


We  can  also  analyze  this  situation  by  a modal  description. 
From  paragraph  2.2,  the  response  of  a plate  to  a point  force  is 


V 


i u L, 


M 


Z 

m 


(x)  tii  (x  ) 
Mtr  ' ymx  s' 

w2  (1+in )-<*>* 
m 


(2.4.8) 


Let  us  assume  that  the  plate  is  simply  supported  so  that  the 
mode  shapes  are 


, 2 ik.x.  -ik.x^  ik.x?  -ik^x? 
^m(x)  W 4 ? (e  ”e  > <e  “e  ) 


(2.4.9) 

There  are  16  terms  in  this  product,  each  representing  a plane 
wave  having  one  of  the  force  wave  vectors  shown  in  Fig.  2.11. 
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The  phase  factors  are  the  "dot  products"  of  these  wave 
vectors  with  the  four  position  vectors  shown  in  Fig.  2.12. 
These  position  vectors  are  the  difference  vectors  between 
the  point  source  and  the  observation  point  and  its  three 
images  in  the  coordinate  axes. 

If  the  drive  point  is  excited  at  frequency  u>,  then  be- 
cause of  the  modal  response  bandwidth,  we  will  assume  that 
all  modes  within  a band  Aw  = irun/2  will  be  excited,  and  that 
modes  outside  this  bandwidth  are  not  excited.  As  we  sum  over 
the  circle  of  excited  modes  in  k-space  (the  quarter  circle  of 
Fig.  2.7  becomes  a full  circle  because  of  the  addition  of  the 
images  of  k in  Fig.  2.11),  some  of  the  terms  will  fluctuate 
wildly  in  phase  as  we  sum  over  the  indices  m^,  while  others 
will  combine  because  the  phase  variation  between  them  is 
small.  The  smallest  ghase  variation  will  occur  for  the 
smallest  vector,  x - xs,  as  seen  in  Fig.  2.12.  We  shall 
assume  that  only  these  terms  contribute  to  the  "coherent" 
part  of  the  summation  in  Eq.  (2.4.8). 

We,  therefore,  replace  Eq.  (2.4.9)  by  its  coherent  part 


Vx)  K(xs}  “ \ ei^  (2.4.10) 


where  R = x - x , and  £ now  varies  over  the  complete  circle  in 
Fig.  2.11.  If  the  phase  variation  in  going  from  one  *attice 
point  (allowable  wave)  to  another  when  summing  over  m is 
less  than  ir/2,  then  the  summation  in  Eq.  (2.4.8)  may  be  re- 
placed by  an  integration  over  angle  in  k-space. 

iu)L>Q  exp  wm  * R cosQ) 

v(x)  * ' “W  £ (l+in/iTJ- Tui”+ii) Ci'+in/J) T <2-4-u> 


where  9 is  the  angle  between  it  and  R. 

Let  us  note  that  M = p A„  = ptt2/AAjc.  We  can  now  write 
AAj<  = k Ak  A6,  and  the  summation  in  Eq.  (2.4.11)  becomes 


55 


v(x) 


i«L 

iaL0  exp  (-~  R cos 6 ) k Ak  A6 

^ i 11+ W5>  > !«„■*■«  (1+in/i) } 


i“l0  ( k dk 

7^  J 0 <<V~“  (Uin/2)  )T^-wTTITIri/;n 

r2n 

X / exp (i—  R cos  6)  de.  (2.4.12) 

J 0 


The  second  integral  is  simply  JQ  (k  R) , while  the  first 
integral  is  approximately 


v(x) 


i Lq  j?  k dk  JQ(k  R) 
WpJ^  wm-w(l+  in/2) 


Vf> 


(2.4.13) 


The  condition  that  the  phase  between  terms  vary  less 
than  tt/2  is  most  severly  tested  when  t and  $ have  the 
same  direction  (cos0»l).  The  number  or  resonantly  excited 
modes  is  tt/2  conns  . The  average  distance  between  these 
inodes  in  k-space,  therefore,  is 


2irk/  ~ o)nng  « 4/ngncb, 

and  our  phase  requirement  is 

4R/nsncb  < tt/2;  R < ' (2.4.14) 
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which  is  consistent  with,  but  more  stringent  than. 

Eg.  (2.4.7). 

The  Bessel  function  J.fwR/c.)  in  Eg.  (2.4.13)  has  the 
asymptotic  form 


/UR  tt  » 

cos  (r—  - j) 

cb  4 


so  that  4 he  mean  sguare  velocity  in  the  coherent  waves  varies 
as  1/R,  wiich  is  the  same  as  our  result  for  the  "direct  field" 
in  the  wave  model  of  the  plate  response. 

The  incoherent  part  of  the  velocity  field  is  found  from 
Eg.  (2.4.8)  by  adding  the  mean  sguare  values  of  modal  response 
incoherently.  The  mean  sguare  value  of  each  term  is 


1 . ZO  _1_  Tr  _J3 

wJn2  2 unM  * 1 M 


- 7 i*n  G/wn 


which  is  the  same  as  the  result  in  Eq.  (2.4.5)  if  we  assume 
that  very  little  dissipation  of  the  direct  field  occurs  before 
the  first  reflection  from  a boundary. 

We  have  shown  here  that  the  direct  and  reverberant  fields 
that  arise  naturally  in  the  wave  analysis  have  their  direct 
counterparts  in  the  "coherent"  and  "incoherent"  components  of 
the  summations  of  the  modal  description.  Such  correspondences 
are  very  useful  in  that  they  allow  us  to  interpret  some  of  the 
phenomena  described  by  a wave  picture  in  terms  of  their  effects 
on  a model  analysis.  For  example,  the  transmission  of 
vibrational  energy  through  a connection  between  two  plates  will 
impose  a degree  of  coherence  between  the  modes  of  the  receiving 


plate.  But  if  the  bulk  of  the  vibrational  energy  may  be 
judged  to  be  in  the  reverberant  field,  we  may  treat  the  modes 
of  the  "source"  and  "receiving"  structures  to  be  incoherent 
with  each  other. 
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FIG  2.3 


FREQUENCY  RESPONSE  OF  IDEAL  RECTANGULAR  FILTER 


on* 


COORDINATES  OF  RECTANGULAR  PUTE 


FIG  2.T 

WAVE  NUMBER  LATTICE  FOR  RECTANGULAR  SUPPORTED  PLATE 


FIG  2.8  \ 

CONTOUR  FOR  EVALUATION  OF  INTEGRAL  IN  EQUATION  1.3.16  j 
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FIXED  surface 


INTERACTION  OF  SINGLE  RESONATOR  A NO  FINITE  PLATE 


? 

FIG  Z,  10 

DIRECT  AND  REVERBERANT  FIELDS  ON  ELASTIC  PLATF 


FIG  2.11 

CONSTRUCTION  OF  4 WAVE  VECTORS 
CORRESPONDING  TO  A SINGLE  MODE 


FIG  2.12 

CONFIGURATION  OF  SOURCE  POINT  AND  THE 
OBSERVATION  POINT  AND  ITS  IMAGES 


CHAPTER  3.  ENERGY  SHARING  BY  COUPLED  SYSTEMS 


3.0  Introduction 


In  SEA,  we  are  principally  concerned  with  systems  that 
may  be  subdivided  into  subsystems  that  are  directly  excited, 
and  other  subsys terns  that  are  indirectly  excited.  For 
example,  a vehicle  excited  by  a turbulent  boundary  layer  may 
be  thought  of  as  a single  system  excited  by  random  noise, 
but  we  are  more  likely  inclined  to  treat  it  as  a collection 
of  subsystems.  One  of  these  will  be  the  exterior  shell  of 
the  vehicle,  which  is  directly  in  contact  with  the  turbulence 
A second  system  will  include  interior  structures,  such  as 
bulkheads,  shelves,  etc.  which  are  excited  through  their  con- 
nections to  the  exterior  shell.  In  this  case,  such  a break- 
down is  immediate  because  of  the  nature  of  these  structural 
elements  and  their  configuration. 

When  the  subsystems  have  been  defined,  we  find  that  the 
effect  of  their  connections  is  to  provide  coupling  between 
them.  This  coupling  is  the  mechanism  by  which  vibrational 
energy  is  transferred  from  directly  excited  to  indirectly 
excited  systems.  If  a modal  description  has  been  applied 
to  the  subsystems,  then  intermodal  coupling  forces  result 
from  the  connection.  If  a wave  description  has  been  used, 
the  coupling  is  more  naturally  represented  by  impedances  at 
the  junction  between  the  structure.  In  certain  cases,  one 
finds  it  convenient  to  use  a mixed  description,  one  of  the 
subsystems  being  treated  as  a wave  system  and  the  other  as  a 
modal  or  resonator  system.  It  is,  therefore,  quite  important 
that  we  be  able  to  shift  our  viewpoint  between  these 
descriptions  as  the  need  arises. 

An  important  key  to  use  of  SEA  is  the  definition  of  sub- 
systems. This  is  a part  of  the  modeling  process  that  we 
shall  deal  with  later  on.  Even  for  "obvious"  choices  of 
subsystems,  the  breakdown  has  an  arbitrary  element  to  it. 

In  many  cases,  structural  elements  are  the  subsystems,  but  in 
others,  modal  (or  wave)  types  may  be  the  key.  The  shell  of 
a vehicle  has  modes  that  represent  mainly  transverse  or 
flexural  vibration  and  modes  of  longitudinal  or  torsional  in- 
plane motion.  It  may  be  necessary  to  treat  transverse  and 
in-plane  modes  as  two  separate  subsystems,  particularly  if 
their  excitation  by  the  environment  or  their  interaction  with 
the  rest  of  the  structure  is  markedly  different. 
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We  begin  the  chapter  by  discussing  the  interaction  of 
(sub)  systems  described  in  modal  terms.  Single  resonator 
systems  are  dealt  with  first,  and  some  important  basic 
theories  of  SEA  are  derived.  Then,  multi-modal  systems  are 
studied.  Various  expressions  for  the  energy  flow  are 
presented  in  this  section. 

In  paragraph  3.1  of  this  chapter,  we  study  interactions 
of  systems  using  a wave  description.  Impedances  and  energy 
densities  become  the  important  interaction  and  dynamical 
quantities  here.  The  definition  of  transmissibility  co- 
efficients is  also  made,  a parameter  that  is  widely  used  in 
acoustical  analysis  of  interacting  systems.  Equivalences 
between  the  parameters  of  modal  and  wave  descriptions  are 
then  drawn.  We  also  show  how  reciprocity  arguments  may  be 
used  to  develop  response  estimates  which  result  in  the  same 
results  as  previously  found.  Reciprocity  methods  are  very 
powerful  in  obtaining  equivalences  between  response  para- 
meters in  mechanical  systems  and  one  of  the  goals  of  this 
chapter  is  to  show  how  they  fit  into  the  general  SEA  frame- 
work. 

We  complete  the  chapter  with  the  analysis  of  some  fairly 
simple  systems  that  have  practical  importance.  This  serves 
to  illustrate  the  methods  developed  in  the  chapter,  and  also 
to  show  how  we  begin  to  develop  a catalog  of  SEA  parameters 
as  a result  of  the  applications  that  have  been  made  of  SEA. 


Energy  Sharing  Among  Resonators 


We  began  our  discussion  of  vibration  in  Chapter  2 with 
a single  dof  system.  Although  quite  artificial,  we  were 
able  to  develop  basic  ideas  from  the  1 dof  system  that  were 
applicable  to  much  more  complex  systems.  This  approach  will  be 
followed  in  the  present  chapter  also.  We  shall  study  the 
energy  interactions  of  the  system  shown  in  Fig.  3.1,  con- 
sisting of  two  simple  linear  resonators  coupled  by  con- 
servative elements.  We  shall  find  that  much  more  complex 
systems  can  be  represented  by  this  system  in  the  following 
sections. 


The  equation  of  motion  of  mass  My  is  given  by  calculating 
all  forces  on  it,  including  those  arising  from  the  motion  of 
mass  M2. 
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In  this  instance,  this  is  best  done  by  using  the 
Lagrangian.  The  kinetic  energy  of  the  system  is: 


KE 


- I M1  % + J m2  * f "c  '*1  + 


The  potential  energy  of  the  system  is: 

PE  * | Kiy|  + | K2y|  + \ Kc(y1-y2)! 

The  equation  of  motion  in  the  absence  of  velocity 
dependent  forces  are  given  by 


d 3KE  9 PE 


Li 


(3. 


Using  Eq.  (3.1.1)  in  Eq.  (3.1.2)  results  in  velocity 
independent  equations  of  motion.  In  addition,  we  include 
the  damping  force  Ry  and  the  gyroscopic  coupling  force  G£, 
to  obtain 


<»i  + ? v yi  + W'W  yi 

■ • 7 Mcy2  (3. 

tM2  * 7 V y2  + V2+,K2+Kc>  ^2 

= L +K  y1-GJ1  - i M yx 


.1) 


.2) 


.3a) 


(3.1.3b) 


i 


These  equations  clearly  display  the  effect  of  motion  in 
one  system  causing  forces  (and  power  input  to  or  absorption 
from)  the  other. 

In  application  to  real  systems,  the  question  invariably 
arises t "Wnat  are  the  uncoupled  structures?"  From  a 
mathematical  viewpoint,  we  may  obtain  uncoupled  equations  hy 
"clamping"  system  (2)  [setting  y2  * Y2  “ Y2  * ® in  E<3*  (3.1.3a)) 
and  by  clamping  system  (1)  in  Eq.  (3.1.3b).  An  alternative 
is  to  reduce  the  coupling  parameter  to  zero,  i.e.,  setting 
Kc  * Mp  * G * 0 in  both  equations.  Either  is  possible,  but 
note  that  the  effective  mass  and  stiffness  of  each  resonator 
will  change  in  the  process  of  decoupling  if  we  use  the  latter 
approach.  Since  we  do  not  want  to  change  the  resonator  we 
deal  with  in  the  process  of  removing  the  coupling,  we  prefer 
the  "clamping"  definition  of  the  uncoupled  system. 

The  clamping  or  blocking  concept  is  also  appealing  from 
an  experimental  viewpoint.  A welded  junction  between  a pipe 
and  a plenum  chamber  for  example  is  not  readily  identified  in 
terms  of  the  coupling  parameters  of  Eq.  (3.1.3).  Neverthe- 
less, we  can  imagine  reducing  the  displacements  of  one  element 
to  zero  by  suitable  clamping  arrangement,  even  though  it  might 
be  a very  difficult  thing  to  do  from  a practical  viewpoint. 

If  stress  variables  are  used  then  the  systems  are  decoupled 
by  being  -,at  apart.  Such  a thought  experiment  is  of  particular 
relevance:  in  the  event  that  one  has  a finite  element  computer 
model  of  each  subsystem  since  the  decoupling  can  be  done  by 
some  simple  instructions.  In  this  monograph  we  will  con- 
sistently define  the  decoupled  forms  of  Eq.  (3.1.3)  or  its 
counterparts  to  be  found  by  setting  the  response  variables  of 
the  other  systems  to  be  zero. 

Eq.  (3.1.3)  displays  more  symmetry  when  rewritten  as 


(3.1.4a) 


y2  + A2y2+W2y2+  x^yi+yyi'Kyil  ■ l2  (3.1.4b) 
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Ai«Ri/wi (Hi+Mc/4) (Ki+Kc) / (M^l^/4) 

y* (Mc/4 ) (Mx+Mc/4 ) “ * (H2+Mc/4 ) “ * 

y-G/  (M;l+Mc/4  ) * , (M2+Mc/4  ) “ i , 
k-Kc/  (M1+Mc/4)  * (M2+Mc/4)i 

and  A - (M1+M<3/4 ) ^ (M2+Mc/4 ) ^ . 


Note  that  the  equations  could  be  made  completely  symmetric 
(that  is,  A would  vanish)  by  further  defining  variables 

yi“yi(Mi+Mc/4)*‘  Also  note  Ai"Li/tt!i+Mc/4)* 

In  Chapter  2,  we  were  interested  in  the  energy  of 
vibration  of  a single  resonator  excited  by  "white"  or  broad 
band  noise.  In  this  Chapter  we  are  concerned  with  the 
response  of  two  resonators  when  both  fci  and  J l?  are  independent 
white  noise  sources.  Since  Eqs.  (3.1.4)  are  linear,  such  an 
analysis  is  fairly  straightforward,  but  it  is  instructive  to  go 
through  it  in  some  detail.  Since  we  want  to  apply  these 
results  to  energy  flow  problems,  we  shall  be  particularly 
interested  in  calculating  the  power  flew  between  the 
resonators . 

Before  we  treat  Eqs.  (3.1.4)  in  all  of  their  detail, 
however,  let  us  consider  a particularly  simple  cases  that 
for  * A,  wi«=a)2*!wnf  A»l,  and  y=*y«0.  Thus,  we  have  two 

identical  resonators  with  stiffness  coupling  only.  If  we 
add  and  subtract  Eqs.  (3.1.4),  we  obtain, 


z±  + Az,  + n*  z ^ = g^ 


(3.1.5a) 


Z _ + Az«  + Q j;  z_  as  g _ 

2 2 ‘2  2 y 2 


(3.1.5b) 
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when 


zryl+y2'  z2=yl"y2'  VV‘2'  g2“W 


, fl|=w£+IC  . 


From  elementary  dynamics,  we  know  that  if  we  start  this 
system  vibrating  without  damping  or  other  excitation  present, 
that  the  general  free  vibration  response  of  either  mass  is  of 
the  form 


y^ttJsA  sin  (fi^t+ij^)  + B sin  (ftjt+ijij)  • (3.1.6) 

If  A and  B are  nearly  equal  in  magnitude  (this  will 
occur  if  one  mass  is  held  fixed  while  the  other  is  displaced 
and  then  released)  then  the  motion  of  either  mass  is  of  the 
form 


yx(t) 


* C sin  (-£  t+ifi _)  sin 
w0  3 


(o»0t+^b)  , 


(3.1.7) 


which  represents  an  oscillation  that  is  modulated  by  an 
envelope  that  varies  at  radian  frequency  k/w0.  Indeed,  this 
beating  phenomenon  is  one  of  the  best  known  and  most  remark- 
able features  of  coupled  oscillators. 

If  the  resonators  have  damping  A,  however,  then  they 
will  have  a decaying  motion  represented  by  an  additional 
"modulation"  of  Eq.  (3.1.7)  by  the  term  e~At/2  [see 
Eq.  (2.1.9)].  If  the  damping  is  very  light,  so  that  k/wq>>A, 
then  the  beating  phenomenon  will  still  be  apparent  but  ir  the 
damping  is  very  strong  compared  to  the  coupling,  so  that 
A>x/(jjq , then  the  energy  of  vibration  will  be  dissipated  be- 
fore the  beating  oscillations  can  take  place.  We  shall  find 
that  quite  parallel  conclusions  may  be  drawn  when  the 
resonators  are  excited  by  a white  noise  source. 

Let  us  now  return  to  Eqs.  (3.1.5)  and  set  &2=0.  Then 
by  examining  cyi^  we  can  determine  the  power  that  is  flow- 
ing into  the  indirectly  excited  resonator.  When  I2-O , then 
91  = 92  an&  t*le  responses  y^,  y2  will  be  identical,  except 
that  their  response  curves  will  be  shifted  as  shown  in 
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Fig.  3.2.  When  there  is  no  coupling,  then  z^  = z£ 
and  #2  * 1/2  (ii  - ^2)  * 0.  When  the  coupling  is  such  *hat 
</o>g<<A , then  tne  spectrum  of  Is  ve*y  nearly  the  dif- 
ference between  the  spectra  S£.  and  as  shown  in 

Fig.  3.2.  The  response  y2  is  weak  in  this  case,  very  broad 
band  but  with  a slight  double  peak  in  its  spectrum  that 
would  show  a slight  beating  effect.  The  response 
y2  * 1/2  (zi+Z2)*zi  has  the  spectrum  *ii”*'*r 

When  the  coupling  is  increased,  the  two  spectra  shift 
farther  apart,  and  the  difference  spectrum  is  increased,  as 
shown.  This  means  that  <$2>  is  increasing  with  the  increase 
in  coupling.  The  spectrum  of  £2  is  even  more  strongly 
double  peaked,  and  fluctuations  at  the  difference  frequency 
</u>q  (bfats)  will  become  more  pronounced.  Also,  since 
z^  and  Z2  are  no  longer  identical,  <y^>  will  be  diminished. 

Finally,  when  the  coupling  is  strong  enough  so  that  the 
resonances  have  shifted  beyond  the  damping  bandwidth 
(k/wq>a) , the  response  #1  and  yo  are  statistically  independent 
since  they  do  not  share  common  frequency  components. 
Accordingly,  ^yf>  « <$-£>  = <if>/ 2.  Thus,  each  resonator  dis- 
sipates 1/2  of  the  power  supplied  by  the  white  noise  source 
&1.  This  power  depends  only  on  the  mass  of  resonator  "1" 
and  not  on  the  coupling  [see  discussion  in  connection  with 
Eq.  (2.1.39)].  Now  the  "beating"  phenomenon  will  be  quite 
strong  implying  an  oscillation  of  energy  in  the  system. 

This  flow  is  in  addition  to  the  steady  power  flow  from 
resonator  "1"  to  resonator  "2"  required  to  supply  the  power 
dissipated  by  resonator  £<yj>>. 

Let  us  now  return  to  the  more  complex  system  represented 
by  Eq.  (3.1.3)  and  (3,1.4).  We  assume  that  both  Lj  and  L2 
are  independent  white  noise  sources.  The  power  supplied  by 
che  source  is  <L^  y^>,  or 


;L1y1>  * (m1+^mc)  <y1y1>  + R1<yi>  + (W^i^i 


- Kc  G (3.1.8) 


Since  y£>  = y|>  = 0 for  stationary  processes,  the 
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first  and  third  terms  on  the  right  hand  side  of  Eq.  (3.1.8) 
vanish.  The  term  R<f£>  represents  dissipation  by  the  damper 
Ri  and  the  other  terms  represent  power  flow  into  the  coupling 
elements. 

In  a similar  fashion,  the  power  supplied  by  source  2 is 


<L2y2>  * R <y|>  + Kc  <yxy2>  ♦ G <y1y2>  + jMc  <yxy2> 


(3.1.9) 


Using  the  fact  that 


! 


A ** 

'Jt  *lV  ■ + * 0 


and 


« ••  •• 

'at  ^1^'  ■ + '^V  * 


we  can  add  Eqs.  (3.1.8)  and  (3.1.9)  to  obtain 


<L1y1>  + <L2y2>  « ^ <y[>  + R2  <y*> 


(3.1.10) 


which  states  that  the  total  power  supplied  is  dissipated  in 
the  two  damping  elements;  the  coupling  elements  in  the  form 
presented  are  non-dissipative. 

We  also  note  from  Eq.  (3.1.8)  or  (3.1.9)  that  the  power 
flow  from  resonator  1 to  resonator  2 is  simply 


ni2  " Kc 


<y2Yi>  - G <y2yi>  " 


Our  next  task  is  to  evaluate  these  averages  in  terms  of  the 
spectral  densities  of  the  sources  and  the  system  parameters. 
This  evaluation  may  be  done  in  the  time  or  frequency  domain. 
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In  the  frequency  domain  calculation,  we  introduce  the  com- 
plex system  frequency  response  N)  defined  as  the  response 

yp  to  a force  Lq  » e“i«t.  Thus,  all  time  derivatives  in 
Eq.  (3.1.3)  are  replaced  by  ♦ -iu>,  and  we  have 


Hu(w)  « (H1+Jmc)“1  (-u2-iwA2+aj|)/D 


(3.1.12) 


H2^(w)  * A (H^+^#SC)"1  (w2u+iwy  + <)/D 


(3.1.12b) 


H12(u)  « (M2+Jmc)“1A“1  (ai2y-iu)Y  + <)/D  (3.1.12c) 


H22(w)  - (M2+^Mc)’1 (-w2-iwA1  + w2)/D  (3.1.12d) 


where  D is  the  system  determinant  given  by 


w'*(l-y2)  + iw3  (A^+A2)-w2  (u^+a)2+A^A2+Y2+2u<) 


ioj  (A^w2+A2wp  + (u}2(jj|-ic2 ) . 


(3.1.13) 


Using  these  transfer  functions,  the  various  second 
moments  of  the  dynamical  variables  can  be  generated  by 
integration  over  frequency  from  -®  to  + «■: 


<y[>  = SLj  f |HUI»  da  + f |h12| 


dw  (3.1.14a) 


<^2>  = S 


L2  f I H22  ^ 2 dU)  + SLX  f lH21^  d“  (3*1-14b) 
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<y*>  ■ sLl_/  ““IHjjI2  dm  + SL^  /L1 Jh12 I 2 du  (3.1.14c) 
m % J oi2|H22j2  du  + SL  y*a)2|H21|2  du  (3.1.14d) 

= SLl  f H21  dw  + SL2  y “iw  H12  H22dw 


(3.1.14e) 


<Vi>  - \ J “ 


H11  H21  dw  + SL2  7 t0*  H12  H22  dw 


(3. 1.14f ) 


H11  K21  a“  + SL 2/iu3  Hh  H22  d" 


(3.1.14g) 


The  final  three  relations  in  Eqs.  (3.1.14  e,  f,  g)  are  the 
terms  we  need  to  determine  power  flow,  the  first  four, 

Eqs.  (3.1.14a,  b,  c,  d)  are  needed  in  evaluating  the  energy  of 
vibration.  In  these  equations  we  have  used  the  "two-sided" 
spectrum  which  has  a range  from  <d=-°°  u=+<»  rather  than 

the  more  "physical"  one-sided  spectral  density  in  cyclic 
frequency  f wnich  has  a range  from  f = 0 to  f = «.  Since  the 
integral  over  the  range  of  each  must  produce  the  same  mean 
square  value,  i.e., 


<y  2> 


S (w)  d to  = / S (f)  df, 
y Jo  y 


(3.1.15) 
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the  relation  between  them  must  be 


\ 

I S(f)  « 2 S (w) 

i I J 


dw 

ar 


4ir  Sy  (w)  . 


(3. 


I 


k 


We  will  not  bother  to  evaluate  the  integral  of 
Eqs.  (3.1.14)  here,  since  they  are  of  standard  form.  Even 
so,  the  form  of  the  relations  is  quite  complicated.  They 
are  found  to  be: 


ttS, 


<y2yl>  = 


w 


ttS. 


W |MC) 


(3.: 


A1A2  [u  (AjWj+Ajwj;)  + y(u)*-w£)  + tc  (A j+Aj)  ] 
X T 

? <M1+  K)  (MC+  iV  d 

»; 


<y2yi> 


ttS, 


W?V 


ttS. 

VV?V 


A^  [Y(A1tf*+A2wp  - k(<*>2-w£)] 

<V  Jmc)4  (M2H  Jmc)4  d 


r 
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k 


.16) 


.17a) 


. irST  irST 

- ) ^ + **2  \ 

1 vK  vk  ( 


M I (A^Wj+AjWp  (u)J+Wj+A1A2+y2+  2yic)-  (A^+Aj)  (w£w|-ic2)] 

(Mi+  k>*  'V  K>i(1-«,2)  d 


“2WSL, 


“2,SL. 


\i  (A1o)2+A2<*)2) 


Vr®  Wc  (Mi+rc>  (m^)  d 


(3.1.17b) 


A1A2 


VVlV  WfV  (M2W  d 


X j-~jt  (A^2  +A2wP  (Wj+Wj+A^  +y2  + 2y<) 


- (A^+Aj)  (a)2u)2-<2  )+<  (A1a)2+A2co2 ) 
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Wj+jM.)  ’ (Mj+Jm0)4  (1-yM  d 


“j”SL, 


Wc 


VK 


ylA^jj+AjW2] 


y— — i i T 

A**  \ 7 /ii  % \ I 


(Mi+rc)  (m2+imc)  d 


(3.1.17c) 


The  quantity  d that  appears  in  these  equations  is  given  by: 


d = Aj  A2  [ (w£-w2) 2 + (A^+A2)  (A^u2+A2wp  ] + y2  (A^io2+A2io£)  2 


+ (y2+2yK)  (A1  + A2)  (A1w2+A2wp  + k2  ( A^+A2 ) 


(3.1.18) 
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These  expressions  are  now  placed  into  Eq.  (3.1.11)  to 
evaluate  the  power  flow.  Incredibly,  they  simplify  enormously 
(although  they  are  still  relatively  complicated) . The  result 
of  the  substitution  is 


ttS, 


ni2  = A 


VMiW 


*SL 

4 

VV? 


(3.1.19) 


where 


A.A,  ( 

A a < y2  [A.wli+A-wII+A.  A,  (A.w^+A-w?] 

(1-y 2 )d  ( 1221121221 


+ (y2+2y<)  +a2U)i^+ic2  ^Ai+A2^| 

(3.1.20) 


and  d is  given  by  Eq.  (3.1.18). 

Comparing  Eqs.  (3.1.4)  to  Eq.  (2.1.1),  we  see  that 
and  A2=w2n2  » where  the  n's  are  th§  loss  factors  of 
each  uncoupled  resonator.  Thus,  when  Xj55  *2  = 0,  the  term  is 


ttS 


Lx(w) 


*Lx(f) 


A1  <Ml+fV 


4o,ln1»'l+J«c) 


- <?!>  <vK> 

(3.1.21) 


according  to  Eqs.  (2.1.39)  and  (2.1.40),  the  energy  of 
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resonator  "1". 


Similarly 


ttS 


L2(o3) 


a2  (m2+^c) 


’ <VK)  <;2> 


(3.1.22) 


t 

I 


is  the  energy  of  resonator  "2"  when  Thus,  Eq. 

(3.1l9)  states  that  the  power  flow  from  resorator  "1"  to 
resonator  "2"  is 


(a)  directly  proportional  to  the  difference  in  decoupled 
energy  of  the  resonators,  where  decoupling  is  de- 
fined by  constraining  the  resonator,  the  energy  of 
which  is  not  being  evaluated 


(b)  since  the  quantity  A is  positive  definite,  the 

average  power  flow  is  from  the  resonator  of  greater 
to  lesser  energy 


(c)  the  quantity  A is  also  symmetric  in  the  system 

parameters  so  that  an  equal  difference  of  resonator 
energies  in  either  direction  (l-*-2  or  2-*l)  will 
result  in  an  equal  power  flow,  we  may  say  that  power 
flow  is  reciprocal. 

The  concepts  of  power  flow  between  resonators  being 
proportional  to  the  difference  in  uncoupled 
energies  is  useful.  When  we  are  exciting  a 
structure  with  noise  so  that  the  input  power  is 
known,  then  since  Sp/M  is  proportional  to  input 
power,  this  interpretation  of  Eq.  (3.1.19)  is  the 
appropriate  one.  However,  we  often  do  not  know  r he 
input  power,  we  only  know  the  m.s.  response  values 
<y^>  and  <y 2 > by  measurements  taken  or.  the  system 
as  it  vibrates  in  its  coupled  state.  We  need  co 
relate  Eq.  (3.1.19)  to  these  valuesalso.  To  do 
this,  we  require  calculations  of  < Yl>r  <yy>  ' 

<y2>  and  <y‘.  from  Eqs.  (3.1.14).  The  ^results 

of  the  calculations  are  as  follows: 
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l<WJM*CftOT5S 


* ]A3wi[  (u>?-w2)2 

(Mj+Jjy2  (w2w2-tc2)d  I 2 2 12 


+ (Aj^+Aj)  (A1w|-»-A2a)|)  ] + u2  (A^j+AjUp  ] 


+ (y2+  2uk)  [w2  ( A1+A2 ) ] - K2  [ (A^+Aj)  (A2“2u2) 


+ (A^2  +A2a2)] 


{M1+K){M2+K)(uH“k2)  d 


X |u2[uToj2^  (A^w2  + A2ap]  + (y2+2vik)  [u2u>2  ( A^+Aj)  ] 


+ k2  (A^A2(A^+A2)  + (A^a^+A2w2) ]|  (3.1.23a) 


(M1+^C)Z  {1"y2)d 


A2  [ (w2-w2) 2+(A1+A2)  (A1a2+A2a2)) 


*•  ia2(a2(A^+A2)  + (A2-2w2)  (A^u)2+a2u)£)] 
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+ (y2+  2ytc)  + icMAj+Aj] 


T r )y2  [A. w^+A^oj? 

(Ml+|Mc)(M2+|Mc)(l-p2)d  ) 1221 


+ A^A2  (A1oj2+A2u^)  ] + (y 2 + 2y<)  [A1u2+A2o)2+k2  (A1+A2)  ] 


(3.1.23b) 

The  expressions  fox  <y2>  and  <y2>  are  found  by  revising  sub- 
scripts in  these  two  equations. 

, We  note  from  Eqs,  (3.1.23)  that  in  coupled  vibration 
<y2>  / w£<y2>  so  that  a detailed  equality  of  potential  and 
kinetic  energy  has  been  lost.  Nevertheless,  by  direct 
substitution. 


(Ml+iWc)  <^1>  ’ (M2+*Mc)  <^2>  " 'VKc)<yi>"(K2+Kc) 


"y 2 > = "V  [Ki)i+(vi2)(iiu2+Vi)] 


A2(M2+iMc> 
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(3.1.24) 


Thus,  we  arrive  at  the  remarkable  conclusion  that  the  power 
flow  is  proportional  to  the  difference  in  the  actual  kinetic 
energy,  or  potential  energy,  or  total  energy.  This'  is  a 
very  useful  result  because  it  allows  us  to  calculate  power 
flow  from  input  power  to  the  system  or  from  measured  vibration. 
We  are  assured  that  we  will  get  the  correct  answer  whether 
we  base  our  calculations  on  the  uncoupled  systems  or  on  the 
actual  vibrational  energy  of  the  coupled  systems. 

If  the  average  energy  of  resonator  "1"  is  and  the 
average  energy  of  resonator  "2"  is  Eo t then  the  power  flow 
between  them  is,  according  to  Eqs.  (3.1.19)  and  (3.1.24) 


"i2  - B(ErE2> 


(3.1.25) 


where 


B= 


y 2 [A1(i)2+A2wi+AiA2  (Aj_w2+A2a,P  ) + (y2+2]jk)  (A1(j2+A2w2)+k2  (A-^+A^ 


(1-y2)  [ (<d2-<d2)  2+(A1+A2)  (A1u2+A2o)2)  ] 


(3.1.26) 


Let  us  suppose  that  resonator  "2"  has  no  direct  excitation,  so 
that  =0.  Then,  the  power  dissipated  by  the  resonator  must 

equal  t£e  power  transferred  from  ”1"  to  "2"?  n2  diSs~n12*,'A2E2 
= B(E1-E2)  or, 


B 

V5 


(3.1.27) 
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This  relation  shows  that  the  largest  value  E2 
can  have  is  which  occurs  when  the  coupling 
(determined  by  B)  is  strong  compared  to  the 
damping 

When  the  two  resonators  are  identical  and  have 
stiffness  coupling  only,  then  from  Eg.  (3.1.26) 


2 Ak  ^ tc  ^ 

B = — = — — (3.1.28) 

4AZ0)q  2Au)q 

and, 


E-  k2/2A2uj2 

E1  1+k2/2A2'.)2 


which  shows  that  e2'*'ei  as  wqA/ic-^0  ar.d  E2+0  as  k-»-0, 
which  is  consistent  with  our  earlier  discussion  of 
this  system. 

We  can  now  make  some  additional  points  to  those  already 
made  in  the  discussion  following  Eq.  (3.1.22).  They  are: 

(d)  the  power  flow  is  also  proportional  to  the  actual 
vibrational  energies  of  the  system,  the  constant 
of  proportionality  being  B,  defined  by  Eq.  (3.1.26) 

(e)  the  parameter  B is  positive  definite  and  symmetric 
in  system  parameters;  the  system  is  reciprocal  and 
power  flows  from  the  more  energetic  resonator  to 
the  less  energetic  one 

(f)  if  only  one  resonator  is  directly  excited,  the 
greatest  possible  value  of  energy  for  the  second 
resonator  is  that  of  the  first  resonator. 
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These  six  points  form  the  basis  for  our  studies  of  energy 
flow  in  mechanf^al  and  acoustical  systems,  since  the  cal- 
culation of  power  flow  between  resonators  according  to  the 
scheumes  indicated  in  Eqs>  (3.1.19)  and  (3.1.25)  will  be 
fundamental  to  our  study  of  more  complex  systems.  Vie  will 
have  to  make  some  additional  assumptions  as  we  proceed,  but 
we  could  not  begin  without  these  (conceptually)  simple 
relationships. 

Before  we  leave  this  subject,  let  us  develop  an 
additional  result  of  modal  interaction  that  will  be  useful 
to  us  in  dealing  with  multi-dof  systems.  We  derive  the 
average  power  into  resonator  2 when  the  only  active  source  is 
Li  and  wt  is  allowed  to  vary  randomly  over  an  interval  Aw 
that  includes  ui2.  To  do  this  calculation,  we  see  from  Eq. 
(3.1.19)  that  for  any  particular  value  of  w^,  the  power  flow 
is  11x2  = A irS^  /A i (M^+1/4MC)  and  that  an  average  over 

requires  that  we  evaluate  <A>W  . As  long  as  the  interval 
Aw  is  not  too  large  compared  io  the  values  of  ui^  and  u>2, 
we  can  replace  both  by  w,  the  center  frequency  of  the 
averaging  band,  everywhere  except  in  the  (wf-<D2)2  term  in  d, 

since  this  is  the  term  that  is  sensitive  to  the  relative 
values  of  and  w2  . We  then  form  the  average 


<i>  = _L_  < i _L_  L.  I .J_ 

° U).  A,A_  (w2-u)|)  2+ui2A2  A,A0  Aoi  2 0)2T 

1 12  12  01,12 


(3.1.29) 


where  we  have  used  the  result  leading  to  Eq.  (2.1.37),  and 
have  defined 


A= (A1+A2) [1+ (y 2w2+y2+2uk+k2/w2 ) /A^ A2 ] 1 . 


(3.1.30) 
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We  see  tha*  A is  very  nearly  the  total  damping  of  the  two 
resonators  when  the  coupling  is  very  weak  but  is  more  com- 
plicated for  stronger  coupling.  The  average  power  flow  then, 
is 


and 


with 


<11, = 
12  w 


L2 


(3.1.31) 


= S 


hi 


M2+  K 

v ¥c 


2Aw 


(1  + X) 


£ 


(3.1.32) 


X = [p2w2+(yz+2MK)+icz/ti)z]/A^A2 


(3.1.33) 


and  we  have  assumed  p2<<l  and  A/w<<1. 

The  average  interaction,  therefore,  is  equivalent  to 
a white  noise  source  acting  on  resonator  2 of  spectral 
density  S£  . This  effective  spectral  density  is  pro- 

portional to  St. , to  the  ratio  of  the  effective  bandwidth 
KA2/2  of  resonator  2,  to  the  averaging  bandwidth  Aw,  to  the 
mass  ratio  of  the  systems,  and  to  the  strength  of  the  system 
coupling  as  measured  by  X(l+X)-*.  The  replacement  of  modal 
interactions  averaged  over  system  parameters  with  white  noise 
sources  is  an  important  device  in  the  development  of  SEA 
methods. 
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3.2  Energy  Exchange  in  Mul ti -Degree-of-Freedom  Systems 


We  now  consider  two  subsystems  that  are  connected 
together*  We  combine  the  analysis  of  distributed  systems 
in  paragraph  2.2  with  that  for  two  single  resonators  in 
paragraph  3.1  to  develop  a theory  of  multi-modal  interactions. 
The  actual  system  that  we  are  interested  in  is  the  situation 
in  Fig.  3.3a,  in  which  the  two  subsystems  are  responding  to 
their  own  excitation  and  the  interaction  "forces".  Each  of 
the  sub-systems  vibrates  "on  its  own"  when  the  other  system 
is  blocked  or  clamped  as  shown  in  Figs.  3.3b,c,  and  can  be 
analyzed  according  to  the  methods  of  paragraph  2.2.  Thus, 
the  equations  for  the  blocked  subsystems  are: 


(b)  + (r1/p1)y<b,+ 


Alyi 


(b) 


= pi/pi  i=l,2 


(3.2.1) 


As  in  Paragraph  2.2,  the  operators  have  eigenfunctions  and 
eigenvalues 


A • ip  • 
l yia 


= „2 


CO  7 V . 

ia  ia 


with,  of  course 


<{p . ip.  „>  = 6 0 

ia  iB  afi 

Pj 


(3.2.2) 


(3.2.3) 


The  boundary  conditions  satisfied  for  \J>.  include  the  clamped 
condition  on  subsystem  1. 

We  now  suppose  that  the  spectral  densities  of  the  forces 
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tia<t)  ’ ia 


is  flat  (white)  over  a finite  range  of  frequency  Aw,  and 
that  within  this  band  there  are  Nj.  « Aw  nodes  of  each 
subsystem,  where  n^(w)  is  the  modal  density  of  the 

ifc^  subsystem.  The  modes  for  these  subsystems  may  be 
illustrated  as  shown  in  Fig.  3.4a.  Each  mode  group 
represents  a model  of  the  subsystem.  In  the  applications 
of  SEA,  this  model  has  some  very  particular  properties,  which 
we  now  list  and  discuss: 

1*  Each  mode  is  assumed  to  have  a natural 

frequency  w^a  that  is  uniformly  probable  over 
the  frequency  interval  Aw.  This  means  that 
each  subsystem  is  a member  of  a population  of 
systems  that  are  generally  physically  similar, 
but  differ  enough  to  have  randomly  distributed 
parameters.  The  assumption  is  based  on  the  fact 
that  nominally  identical  structures  or  acoustical 
spaces  will  have  uncertainties  in  modal  para- 
meters, particularly  at  higher  frequencies. 

2.  We  assume  that  ev?ry  mode  in  a subsystem  is 
equally  energetic,  and  that  its  amplitudes 

Yia  (t)  ~f  piYi*i«dxi/M 


are  incoherent,  that  is, 


<YiaYi8\  " 


This  assumption  requires  that  we  select  mode 
groups  for  which  this  should  be  approximately 
correct,  at  least,  and  is  an  important  guide 
to  proper  SEA  modeling.  It  also  implies 
that  the  excitation  functions  L,  are  drawn 
from  random  populations  of  functions  that 
have  certain  similarities  (such  as  equal  fre- 
quency and  wavenumber  spectra)  but  are 
individually  incoherent. 

3.  As  a matter  of  convenience „ we  may  also  assume 
that  the  damping  of  each  mode  in  a subsystem  is 
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the  same.  This  is  not  essential , but  it 
greatly  simplifies  the  formalism  and  tends  to  be 
nearly  true  for  reasonably  complex  subsystems. 

The  conditions  just  described  are  the  basis  for  the 
word  "statistical”  in  SEA.  The  concept  of  systems  derived 
by  selection  of  individual  modes  and  modal  excitations  from 
random  populations  is  of  greater  importance  than  the  use  of 
random  (or  noise)  excitation.  Indeed,  we  can  use  SEA  to 
good  accuracy  in  situations  in  which  the  excitation  is  a 
pure  tone  if  there  are  a sufficient  number  of  modes  in  inter- 
action to  provide  "good  statistics".  We  shall  see  an  example 
of  this  in  Chapter  4. 

We  now  "unblock"  the  system  and  consider  the  new 
equations  of  motion,  which  are 


[i*j;  i» j = 1,2]  (3.2.4) 


One  now  expands  these  2 equations  in  the  eigenfunctions 

iK  (x, ) to  obtain 
lot  i 


£«w«aVWa*c.oV  (3-2-5a) 


a a a 


W,  Ko<,V  <3-2-5b> 

a 
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where  we  now  reserve  the  indices  a for  subsystem  1 
and  o,t,...  for  subsystem  2,  and  we  have  also  set  A]*rj/p£. 
The  mass  of  subsystem  i is  in  these  equations.  The 
coupling  parameters  are 


ao  i 


y12*xl'x2*xl^atxl^o  tx2*xi*l  dx1#etc. 


boundary 


(3.2.6a) 


a a 


boundary 


y21fxl(x2*  #x2^0tfxi'x2)^o(x2)dx2'  etC* 

(3.2.6b) 


where  the  integrations  are  taken  along  the  boundary  between 

the  subsystems  and, there fore,  over  the  same  range  of 

xi , x2  in  both  integrals.  Conservative  coupling  requirements 

are  met  by  u12”1‘n'  Y12“Y21'  and 

k12  =k21  ’ 

The  coupled  systems  may  now  be  represented  as  shown  in 
Fig.  3.4b,  with  the  interaction  lines  showing  the  energy 
flows  that  result  from  the  coupling.  Suppose  we  are 
interested  in  the  energy  flow  between  resonator  (mode)  a of 
subsystem  1 and  mode  a of  subsystem.  2.  We  showed  at  the  end 
of  paragraph  3.1  that  energy  modal  pair  interaction,  when 
averaged  over  the  ensembie  of  systems  will  act  like  a white 
noise  generator.  Thus,  modes  a and  a have  energies  Ea  and  Ea 

(yet  to  be  determined)as  a result  of  these  noise  generators. 
Further,  the  modal  energies  of  subsystem  1 modes  are  all 
equal,  so  that  Ea  * - const,  and  Ea  * E2  * constant.  (Note 

that  E,  and  E,  are  modal  energies).  Under  these  circumstances 
we  have  for  tne  inter-modal  power  flow 


<B  > (E.-E-) 

a0  wa’wo  1 2 
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f 


(3.2.7) 


according  to  Eg.  (3.1.25).  The  average  value  of  is  found 

in  a manner  completely  analogous  to  the  result  in  Eg.  (3.1.29) 
to  be 


TT  12 

<Baa>  * 2 Aw  <^>ol,o  (3.2.8) 

V“o 


where  \ is  defined  by  Eg.  (3.1.33)  and  the  average  is  taken 
with  respect  to  freguencies  wa,w^. 

The  total  power  flow  from  all  modes  of  subsystem  1 
to  mode  o of  subsystem  2 is,  therefore. 


"l.a  \<Baa>  WE2> 


Finally,  the  total  power  from  subsystem  1 to  subsystem  2 is 
found  by  summing  over  the  N2  modes  of  subsystem  2; 


JI 


12 


nin2  (ErV 


7T 


Aw 


y2w2 


+ (v2+2ytc)  +k2/w2 


(ErE2> 


(3.2.9) 


where  6 w^  = Aw/N,  is  the  average  freguency  separation 
between  modes. 

Several  interesting  features  may  be  seen  in  this 
result.  The  first  is  that  the  power  flow  is  proportional 
to  the  bandwidth  Aw,  a result  that  agrees  with  our 
intuition.  Secondly,  the  power  flow  is  proportional  to  the 
difference  in  average  actual  modal  energies  of  the  two 
subsystems,  which  by  now  we  also  expect.  We  recall  from 
Fig.  3.2  and  Eg.  (3.1.7)  that  jc/w  is  the  freguency  shift 
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produced  by  stiffness  coupling  alone.  We  may  suspect  (and 
more  detailed  analysis  confirms)  that  the  numerator  of  the 
fraction  in  Eg.  (3.2.9)  is  the  square  of  the  frequency 
deviation  produced  by  the  coupling.  The  ratio  of  this  to 
the  square  of  the  geometric  mean  of  the  average  modal 
spacings  of  the  two  systems  is  a measure  of  the  strength  of 
the  power  flow. 


Let  us  define  the  total  energy  of  the  subsystems  by 


El,tot  and 

VE2,tot/N2' 


E-  . . . Then,  clearly  since  E,=E,  . ,/N,  and 

4|tOt  1 1 ftOt  1 


> N.N, 
12  a a 12 


1,  tot 
N, 


2,  tot 
N- 


= WT) 


12 


El,tot  " N2  E2,tot 


(3.2.10) 


where  r|i2S<Baa>N2//°J  * V7e  define  r,2l”Nlr|12/,,N2 ' then 


n 


12 


w (n12El,tot"n21E2,tot) 


(3.2.11) 


The  quantities  m2  ahd  noi  ar®  called  the  coupling  loss 
factors  for  the  systems  1,  2.  By  introducing  them,  we  have 
lost  the  symmetry  of  the  energy  flow  coefficients  of  the 
preceding  equations,  but  the  coupling  loss  factor  has  strong 
physical  appeal.  Notice  that  wr}]_2  Ei  tot  represents  the 

power  lost  by  subsystem  1 due  to  coupling,  just  as  the 
quantity  E^  tQt  represents  the  power  lost  to  subsystems 
1 by  damping,  as  measured  by  the  damping  loss  factor  n =A/u>. 
Of  course,  if  system  1 is  connected  to  another  system  with 


91 


energy  E,  . , it  will  receive  an  amount  of  power  E,  • 

The  bisiHgfationship  2'tot  21 


nln12  * n2n21 


(3.2.12) 


is  very  useful  in  practical  situations  of  power  flow  calcul- 
ations . 


We  may  now  use  these  results  to  calculate  the  still  un- 
known system  energies.  Consider  the  system  of  Fig.  (3.4b)  to 
be  more  simply  represented  as  shown  in  Fig.  3.5.  Then  the 
power  flow  equations  for  systems  1 and  2 respectively  are 


nl  ,in“ni  ,diss+ni2*w  [r)lEl , tot+n12El , totfn21E2 , tot1 


(3.2.13a) 


n2,ins*II2,diss+II2l“altri2E2,tot+n21E2/tot"ri12El,tot] 


(3.2.13b) 


First,  consider  the  case  where  only  one  of  the  systems 
is  directly  excited  by  an  external  source.  Set  IT-  . =0,  and 
accordingly,  from  Eq.  (3.2.13b)  we  get  ' 


E2,tot  = n12  = n21 

El,tot  n2+n21  ni  n2+n21 


(3.2.14) 


Which  is  the  multi-modal  equivalent  of  Eq.  (3.1.27).  A 
solution  of  Eq.  (3..'»13)  in  terms  of  both  power  sources  is 


E1  * ("i,in  'Vn* 


n2,inn21)/  “D 


(3.2.15a) 
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(3.2.15b) 


®2  = {I12,in(nl+n12)  + ni,inn12)/wD 


D m (n^-r  n12  ^n2+r,21^  “n12n21  * 


(3.2.15c) 


The  coupling  loss  factor  or  its  various  equivalent 
expressions  is  a measure  of  inter-modal  forces  at  the  system 
junction,  averaged  over  frequency  and  over  the  modes  of  the 
interacting  systems.  Sometimes  this  calculation  may  be 
carried  out  directly,  but  more  often  it  is  found  useful  to 
adopt  another  view  of  the  interaction  - that  developed  by 
the  impingement  of  waves  on  the  boundary.  Consequently,  we 
now  proceed  to  a discussion  of  the  coupling  of  wave  bearing 
systems . 


3.3  Reciprocity  and  Energy, Exchange  in  Wave  Bearing  Systems. 


In  Paragraph  3.2,  we  described  the  vibrations  of  the 
systems  that  we  are  examining  in  terms  of  modes.  An 
alternative  formulation  of  the  problem  in  terms  of  waves 
exists  that  sometimes  offers  both  conceptual  and  practical 
computational  advantages  over  a modal  analysis.  This 
approach  was  touched  on  in  paragraph  2.3,  wherein  we  found 
the  dispersion  relation  between  wavenumber  and  frequency 
in  the  form  of  Eq.  (2.3.2).  We  examined  the  admittance 
function  for  a flat  plate,  both  finite  and  infinite,  to 
a point  force  transverse  to  the  plane  of  the  plate. 

In  paragraph  2.3,  we  also  discussed  the  vibration  of 
a simple  resonator,  shown  in  Fig.  2.9  as  a result  of  its 
attachment  to  a plate  undergoing  random  vibration.  We 
found  that  the  energy  of  the  resonator  was  given  by 
Eq.  (2.3.24),  which  is  identical  in  form  to  Eq.  (3.2.14), 
where  (the  resonator  is  subsystem  2 and  the  plate  is 
subsystem  1)  , one  has  n =1  (one  resonator  mode)  , and 

El,tot=Mp<vp>  and  ncoup  - r>21‘  Thus,  the  quantity  ncoup 
that  we  calculated  on  the  basis  of  an  interaction  of  the 
resonator  with  many  modes  of  the  plate  is  just  the  coupling 
loss  factor  of  Paragraph  3.2.  But,  it  is  also  the  same  as 
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would  be  calculated  on  the  basis  of  the  resonator  inter- 
acting with  an  infinite  plate  [See  Eg.  (2.2.23)].  Consequently, 
it  is  often  advantageous  to  develop  the  equations  for  coupl- 
ing loss  factors  on  the  basis  of  interaction  with  an  infinite 
system,  since  such  calculations  may  be  much  simpler  than  they 
would  be  on  a modal  basis. 

He  should  also  make  another  point  about  our  use  of  the 
term  "blocked"  as  it  applies  to  both  interacting  modal  and 
wave  systems.  He  define  the  blocked  system  to  be  that 
system  that  results  when  the  other  system  has  a vanishing 
response  for  all  time.  Thus,  for  example,  when  a sound 
field  described  by  its  pressure  response  interacts  with  a 
structure  described  by  its  displacement,  the  "blocked" 
sound  field  is  that  field  that  occurs  when  the  structure  is 
absolutely  rigid.  The  "blocked"  structural  vibration  is 
that  vibration  which  occurs  when  the  pressure  vanishes,  i.e., 
the  in  vacuo  condition  of  the  structure.  We  would  ordinarily 
think  of  such  a vibration  as  "free",  but  we  must  distinguish 
between  the  "blocked"  or  "uncoupled"  interaction,  and  our 
ordinary  ideas  of  free  or  constrained  boundary  conditions. 

The  interpretation  depends  on  the  variables  that  we  use  to 
define  "response"  in  the  adjoining  system. 

He  saw  in  Eq.  (2.4.9)  and  in  Fig.  2.11  that  we  could 
form  four  simple  waves  out  of  an  eigenfunction  of  the 
rectangular  plate.  Although  such  a replacement  is  not 
universal,  it  is  useful  to  think  of  it  as  a typical  example 
of  wave-mode  duality.  The  wave  number  lattice  for  normal 
inodes  such  as  shown  in  Fig.  2.7  may  be  re-interpreted  as 
a distribution  of  waves  by  reflecting  the  lattice  points 
in  the  various  coordinate  axes  (or  planes)  of  the  wave 
number  space.  Incoherence  of  modes  means  incoherence  of 
waves  in  different  directions.  Equality  of  modal  energy 
is  translated  into  equality  of  intensity  for  wave  groups 
contained  by  equal  angles  of  spheres  or  circles  in  k-space 
as  shown  in  Fig.  3.6.  A cylinder  has  a more  complicated 
locus  for  bending  modes  than  that  shown  in  Fig.  3.6  but  in  this 
case  also,  the  equivalence  between  waves  and  modes  is 
essentially  as  shown  here. 

The  arguments  for  coherence  and  power  are  based  upon 
the  decomposition  of  modes  into  waves  as  developed  in 
Paragraph  2.4.  The  relation  between  modal  energy  and 
- power  is  obtained  by  simply  noting  that  equal  modal  response 
amplitudes  will  result  in  equal  amplitudes  of  each  wave  in 
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the  decomposition  of  Eq.  (2.4.9).  Secondly,  if  the  modal 
responses  are  incoherent,  the  wave  motions  for  adjacent 
points  on  the  wave  number  circle  in  Fig.  3.6  will  be  in- 
coherent. There  is  still  the  possibility  that  there  will 
be  appreciable  coherence  between  waves  in  different 
directions  at  any  point,  since  there  is  functional  relation- 
ship between  the  phase  of  any  one  wave  component  of  a mode 
and  the  3 (for  a two-dimensional  mode)  other  components. 
However,  as  noted  in  the  discussion  preceding  Eq.  (2.4.10), 
tills  phase  coherence  is  quite  small  unless  the  source  and 
observation  points  are  close  together.  If  the  source  may 
be  assumed  to  be  randomly  located  over  the  surface  (or 
volume)  of  the  system,  even  this  coherence  will  vanish. 

We  imagine,  therefore,  a diffuse  wave  field  (in  the 
SEA  sense)  possessing  energy  distributed  uniformly  over 
the  frequency  interval  Au,  The  energy  density  (Energy  per 
unit  "area")  of  this  field  is  a£  . The  various  parcels 
of  energy  are  carried  by  the  waves  with  an  energy  velocity 
Cg  as  defined  by  Eq.  (2.3.7),  and  the  power  is  distributed 
among  the  various  directions  according  to  an  intensity 


d(Ai)  - a£  Ci  d (si)  dft/n. 

9 ^ 


(3.3.1) 


where  d(Ai)  is  the  intensity  of  wave  energy  in  the  interval 
of  directions  dft,  D(ft)  is  a weighting  function  that 
essentially  measures  the  distribution  of  area  in  the  in- 
terval between  k,  k + Ak  as  a function  of  ft  in  Fig.  3.6(b), 
and  ftt  is  the  total  range  of  the  angle  ft.  One  must  have 
<D(ft)>o  s li  and  it  often  occurs  that  D (ft) = 1 (this  is 
called  "isotropy"). 

Reciprocity.  There  is  a general  principle  that  applies 
to  systems  composed  of  linear,  passive  and  bilateral  elements, 
called  reciprocity,  that  is  quite  useful  in  our  discussions 
of  wave  interactions.  If  we  imagine  for  the  moment  that  the 
system  is  cut  up  into  a very  large  number  of  tiny  masses, 
springs  and  dashpots  (dampers)  these  descriptors  mean 
the  following: 

- (a)  linear:  the  mechanical  response  of  each  element 

is  directly  proportional  to  the  force  (including 
its  sense)  causing  the  response 
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(b)  passive:  the  only  sources  operative  are  those 

explicitly  set  aside  as  sources  in  the 
equations  of  motion;  no  element  of  the  system 
can  generate  energy. 

(c)  bilateral:  in  transferring  forces  from  one 

neighbor  to  the  next,  a reversal  of  roles 
between  the  neighbors  as  to  force  interactions 
will  result  in  an  exact  reversal  of  the 
relative  motions.  Notice  that  gyrators  are 
excluded  from  the  system  at  this  level,  even 
though  the  interaction  of  system  modes  may  end 
up  as  a gyroscopic  interaction. 

The  statement  of  reciprocity  may  be  made  with 
references  to  Fig.  3.7.  If  at  any  two  terminal  pairs  of 
a reciprocal  system  we  generate  a drop  i at  terminal  pair 
1 and  measure  a flow  U through  a "wire"  connecting  the 
terminal  pair  2,  then  if  a drop  p*  is  applied  to  terminal 
pair  2,  a short  circuit  flow  v'  will  be  developed  at 
terminal  pair  1.  The  statement  of  reciprocity  is 

v'  U 

pr  l (3.3.2a) 


The  only  restriction  is  that  p',  U and  A,v'  be  con- 
jugate variables,  that  is,  their  time  average  product 
equals  the  power  flow  appropriate  to  that  set  of  terminals. 
Clearly 


(3.3.2b) 


and  if  p’  and  1 are  noise  signals  with  identical  spectral 
shapes,  then 


i 


<y*  2>  = <ua> 
<P,2>  <lZ> 


(3.3.2c) 
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Similar  statements  are  easily  derived  for  the  systems 
show.*  in  Figs.  3.8  and  3.9.  They  are 


<p2>  x <V  2> 
<v2>  <u,2> 

for  system  B,  and 

<J2>  B <V  2> 
<v2>  <p,2> 


(3.3.3) 


(3.3.4) 


for  system  C.  These  various  statements  have  their  applications 
depending  cn  whether  it  is  more  natural  for  us  to  apply 
forces  or  velocities  at  one  location,  and  measure  forces  or 
velocities  at  the  second  location. 


Systems  Connected  at  a Point.  As  a simple  example  of 
reciprocity  consider  the  force  exerted  on  a point  clamp  at 
the  edge  of  a plate,  as  shown  in  Fig*  3.10,  due  to  vibration 
of  the  plate.  Let  £(t)  be  a band  limited  noise  load  applied 
to  an  arbitrary  point  xa  on  the  surface  of  the  plate,  and  let 
be  the  load  that  results  on  the  clamp.  The  load  Jt  will 
inject  an  amount  of  power  <£2>Gp  into  the  plate  (Gp  is  the 
mechanical  conductance  for  the  plate)  and  this  will  result 
in  a m.s.  plate  velocity  <v^>  given  by 


<v2>  * <£2>G7<i)n  M 
P P P P 


(3.3.5) 


where  np  is  the  loss  factor  of  the  plate  modes  and  Mp  is  the 
plate  mass.  Now  we  say  that 


<i2>  = r <v2> 


(3.3.6) 
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where  r is  the  parameter  that  we  seek. 

In  the  reciprocal  situation , we  drive  the  clamp  with 
a prescribed  noise  velocity  v'ij(t)  that  has  the  same  band 
limited  spectrum  that  £(t)  has.  This  velocity  source 
generates  an  amount  of  power  <v'|>  R^,  (R.  is  the  input 
resistance  at  the  edge  of  the  plate)  which  results  in  a plate 
m.s.  velocity  <v'p>  given  by 


<v,2>  = <v,2>  R^wripM,  (3.3.7) 


If  we  now  apply  the  reciprocity  statement  of  Eq.  (3.3.4), 
or  Pig.  3.9,  we  obtain  <i2>/< A* > * <v^2>/<v^2>,  or 


r 


VG 


1,< 


(3.3.8) 


Let  us  now  extend  this  example  by  connecting  a second 
plate  to  the  first  at  the  edge  point  under  discussion, 
requiring  that  they  move  together  in  a transverse  way  but 
(to  keep  the  matter  simple)  no  moment  is  transmitted.  Again, 
we  let  £i(t)  be  the  band  limited  noise,  that  generates  a 
m.s.  velocity  on  plate  1 given  by  <v£>.  We  are  interested 
in  finding  the  velocity  <v2>*r'  <v?>  that  results  from  this 
motion.  Clearly  1 


<v|>  ■ T'  <v£>  » r'  <fc,2>  (3.3.9) 


If  we  now  apply  the  load  ££(t)  (same  spectrum  as  A)  at  the 
position  where  V2  was  measured,  we  will  inject  an  amount  of 
power  <Jt2,2>  G2  « Into  plate  1 through  the  junction  (call 
it  lU^)  and  a certain  amount  will  be  dissipated.  From  our 
earlier  definitions  of  loss  factors  and  coupling  loss  factors, 
this  fraction  is 


<*J2>  G2 


,»  n21/(n21 
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(3.3.10) 
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This  power,  transmitted  into  plate  1 will  produce  a m.s. 
velocity 


<vj2>  = n^/^n^. 


(3.3.11) 


Applying  the  reciprocity  condition. 


r' 


°2.»  n21 
°1,-  Vn21 


> 


<v2> 


(3.3.12) 


If  the  conditions  are  such  that  65  « * (ir/2(n2/M2)  and 
Gi,®  =»  (tt/2)  (n1/M1) . Then  Eq.  (5l3.ll)  can  Be  rewritten 
as* 


m2<v|>  m1<v|>  n21 

n2  ~ Vn21 


(3.3.13) 


which,  of  course,  is  the  same  result  as  Eq.  (3.2.14).  In  a 
sense,  this  offers  an  alternative  derivation  of  the  basic 
power  flow  relations  of  SEA  and  it  avoids  discussion  of 
either  modes  or  waves.  Of  course,  to  evaluate  * and 
G2  a.  , one  must  have  either  a modal  or  wave  model.  Also, 
to  evaluate  H2i/  one  must  either  do  a modal  analysis, 
following  the  chain  from  Eqs.  (3.2.4,  5,  6),  through 
Eq.  (3.2.9)  and  (3.2.10),  or  proceed  by  wave  analysis. 

We  can  relate  the  power  flow  from  plate  1 to  plate  2 
to  the  blocked  force  given  by  Eqs.  (3.3.6)  and  (3.3.8), 
by  subdividing  the  situation  shown  in  Fig.  3.12a  with  the 
two  problems  shown  in  Fig.  3.12b.  The  actual  force  1 
that  the  function  applies  to  plate  1 is  equal  to  the  blocked 
force  £3,  less  the  force  induced  by  the  velocity,  Jti*vZ^. 

Since  our  convention  is  that  upward  forces  and  velocities 
are  positive,  the  upward  velocity  v imposed  at  the  boundary 
will  result  in  upward  force  vZ,  on  the  edge  of  plate  1.  The 
motion  v will  also  produce  an  upward  force  VZ2  on  plate  2 
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and  a conisquent  downward  reaction  force  on  the  edge  of 
plate  1.  Consequently^ 


l ■ - vZj  * 4g  + vZ^ 


or 


<i2>  « <v2>  (Zj+Zjl2  (3.3.14) 


We  now  relate  <i£>  to  <v|>  by  Eq.  (3.3.6).  We  can 
relate  <v2>  to  <v|>  by  a relation  similar  to  that  of 
Eq.  (3.3.7); 


<v|>  * <v2>  Rj/  o)M2n2  • (3.3.15) 

We  may  then  combine  Eqs.  (3.3.6,  8,  14,  15)  to  obtain 


2 R1R2  \ 1 

" |Z1+Z2|Mwn2  (3-3'16) 

In  many  ways,  Eq.  (3.3.16)  is  the  multi-dof  equivalent 
of  Eq.  0.1.19)  since  M<v|>  is  determined  in  Eq.  (3.3.16)  by 
the  power  injected  into  system  1,  and  the  power  flow  co- 
efficient in  the  relation 


M2<v2>  m1<v1>  ( 
n2  " nl  ( 


a12 


_2 

irwn^ 


R1R2 


'Zl+Z2 


(3.3.17) 
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is  the  "coupling  loss  factor"  for  the  uncoupled  system 
energies.  Tt  satisfies  the  same  reciprocity  relation  that 
doe 8 - namely 


nla12  * n2a21 


(3.3.18) 


The  quantity  aij?  will  take  on  various  forms  depending 
on  the  nature  of  the  interacting  systems,  and  its  relation 
to  will  depend  on  the  system.  We  shall  come  back  to 

this  in  paragraph  3.4. 

The  result  of  Eq.  (3.3.17)  for  the  system  in 
Eq.  3.11  is  typical.  Sinca  n^.2  must  ultimately  be  deter- 
mined by  aio,  the  coupling  loss  factor  is  usually 
representable  in  terms  of  certain  junction  impedances  or 
averages  over  such  impedances.  If  we  had  allowed  moment 
constraints  between  the  systems,  then  our  impedances  would 
have  become  matrices  involving  both  force  and  moment  terms, 
but  the  structure  of  the  result  would  remain  the  same.  The 
impedances  that  enter  Eq  . (3.3.17),  Z^R^-iXi  and 
Z2=R2-iX2  may  be  evaluated  on  either  a modal  or  wave  basis. 

As  we  saw  in  Chapter  2,  there  are  circumstances  in  which  a 
junction  impedance  can  be  quite  simply  calculated  on  a wave 
basis  if  resonance  frequencies  are  averaged  over  a band  of 
frequencies  and  if  the  drive  point  is  also  randomly  located. 
Thus,  one  important  use  of  wave  analysis  in  SEA  is  in  the 
calculation  of  junction  impedances. 

Systems  connected  alon^  a Line.  The  coupling  loss 
factor  Is  fairly  asy  to  derive  when  the  systems  are  con- 
nected at  one  or  more  points.  Another  relatively  tractable 
situation  is  when  the  connection  is  along  a line  (for  two 
dimensional  systems)  or  along  a surface  (for  three  dimen- 
sional systems),  if  the  dimensions  of  the  line  (or  surface) 
are  large  compare^  to  the  length  of  a free  wave  in  the 
system.  In  this  case,  a set  of  waves  in  the  interval  d ft, 
shown  in  Fig.  3.13  will  be  partially  transmitted  and  partially 
reflected  at  the  junction.  In  acoustics,  the  ratio  of  the 
transmitted  to  incident  power  is  called  the  transmission 
coefficient  f(ft)  and  will  in  general  depend  on  the  direction 
ft  that  the  incident  waves  are  travelling  with  respect  to 
the  bounding  line  or  surface. 
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By  our  hypothesis,  therefore,  the  total  transmitted 
power  will  be 


(3.3.19) 


where  d(AI)  is  given  by  Eq.  (3.3.1),  ftinc  is  the  range  of 
angles  of  propagation  that  impinge  upon  the  boundary,  and 
Lp(ft)  is  the  projection  of  the  boundary  length  (or  area) 
presented  to  the  waves  travelling  in  direction  ft.  Since 
n^2  is  the  transmitted  power  due  to  energy  > where 

A1  is  the  area  (or  volume)  of  system  1,  we  have 


n12 


12 

to  A^  A^ ^ 


L (ft)D(ft)dft. 
p (3.3.20) 


For  example,  in  a 2-dimensional  isotropic  system,  let 
the  transmitting  boundary  be  a straight  line  of  length  L 
and  ft  be  the  angle  of  wave  incidence,  as  shown  in  Fig.  3.14. 
Then  D(ft)  * 1,  ft.  * ir,  ft^  = 2t 


n12 


c L 

a-: 

2tiu>A1 


tt/2 

T (ft) 

ti/2 


cosft  dft  * 


c L 


<T (ft)C0Sft>fi 


(2-dim' 11  (3.3.21) 


Of  course,  the  transmissibility  must  be  known,  but 
this  information  is  usually  available  from  existing  sources, 
or  is  calculable.  In  the  3-dimensional  acoustic  case,c*Cg* 
speed  of  sound?  A^V^,  the  room  volume  and  L^AW,  the  area 
of  the  adjoining  partition.  Also  ftt=4ir  (steradians)  and 
ft^=2n  steradians.  Then, 


- ^L 
^12  " 2wV, 


<t  (ft) cos^>  ^ [3-dim'l] 


(3.3.22) 
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where  <j>  is  the  angle  between  the  wave  vector  and  the 
normal  to  the  panel,  and  -10  log  <t (fljcosij)^  is  the 
ordinary  transmission  loss  (TL)  of  building  acoustics. 

In  this  situation,  therefore,  we  would  be  advised  to  evaluate 
n12  from  existing  data  sources  on  wall  transmission  loss. 

3.4  Some  Sample  Applications  of  SEA. 


Here,  we  present  some  applications  of  SEA  that  demon- 
strate its  usefulness  in  relatively  simple  situations  that 
are  of  practical  interest.  The  emphasis  here  is  in  the 
use  of  relations  developed  in  paragraph  3.3,  rather  than  in 
development  of  the  theory.  Of  course,  every  new  application 
will  generate  formulas  or  data  for  the  parameters  that  we 
are  interested  in,  such  as  modal  densities,  coupling  loss 
factor,  and  damping  loss  factor. 

Resonator  Excited  by  a Reverberant  Sound  Field.  As 
the  first  example  we  consider  a resonator  formed  by  the 
small  piston  of  mass  Mq  in  the  wall  of  a room  as  shown  in 
Fig.  3.15.  The  piston  is  supported  by  a spring  of  stiffness 
Mwq  and  a dashpot  having  a mechanical  resistance  UqHqMq. 

We  shall  call  the  room  system  "R"  and  the  resonator  system 
"0".  Then,  according  to  Eg.  (3.2.14) 


fo  c _1 n0R 

ER  " VnOR 


(3.4.1) 


where  <0R  is  the  coupling  loss  factor  from  the  resonator  to 
the  room.  We  evaluate  n0R  by  using  IJq.  (3.3.17).  The 

impedance  ZR  is  the  radiation  impedance  of  the  piston 
"looking  into"  the  room 


R , - io  M , 
rad  rad 


(3.4.2) 


These  quantities  are  well  known  for  the  piston  radiator.  The 
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impedance  "looking  into"  the  resonator  was  given  in  Eq. 
(2.1.19), 


zo  * VoMo  ■ iw  Moa-wo/a,2) 


(3.4.3) 


thus,  averaging  a (defined  by  Eq.  (3.3.17))  over 
the  band  Aw  gives  UK 


wft<a„>  A 
0 OR  Aw 


2Aw  ( ; 

11  \ 


R jW-nAM- 
rad  000 


0 ' " 0‘"0+Rrad‘i“  (M0+Mrad)  +i“o V“ 1 * 


w n p 

0 0Krad 

Vo  M0+Rrad 


(3.4.4) 


and  substituting  wQ  <a>  into  Eq.  (3.3.16)  , 


M 


0 


<v0> 


2al 

pc: 


rad 

Vo‘VRrad 


(3.4.5) 


where  we  have  used  n2*nA»  Comparing  Eqs..  (3.4.5)  and  (3.4.1), 
we  see  that  J 


nOR  * Rrad/“M0  ' (3-4-6) 

The  coupling  loss  factor  for  the  piston  resonator,  therefore, 
is  determined  by  the  radiation  resistance  of  the  piston 
"looking  into"  the  room.  When  Rra(j>>wonoM0  * the  ener9y 
of  the  resonator  will  equal  the  average  umodel  energy  of 
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the  sound  level 


Two  Beams  in  Longitudinal  Vibration.  As  a second 
example,  consider  the  problem  o t two  long  beams  connected 
by  a fairly  weak  spring,  as  shown  in  Fig.  3.16.  If  beam  1 is 
directly  excited  to  energy  Mi  <v|>  and  beam  2 is  only 
excited  through  the  connection,  then  by  our  assumptions, 
the  energy  flow  will  be  fairly  weak  into  beam  2.  We 
may  use  the  result  expressed  by  Eq.  (3.3.16)  diroctly, 
where  we  define  system  1 to  the  beam  1 and  the  spring, 
and  system  2 to  be  beam  2 alone.  Then,  the  average 
mechanical  impedance  looking  into,  beam  2 is 


(P2c2A2) 


1 

7 


P2C2A2*Z2“R2 


(3.4.7) 


The  impedance  looking  into  system  1 is  complex 


Z 


1 


Rl^iVa))  iK  (K/u)) 2 

' » *1,00 


(3.4.8) 


Using  Ri~Re(zi) 


and  K/w  « Ri>00,  R2  a> 


we  get 


M2<v2> 

n„ 


Ml<Vl> 


z 4K2/m2 (3.4.9) 

7rwnn  P1C1A1»P2C2A2 


In  this  case,  ai2“T1i2  (n'odal  energy  of  system  2 <<  modal 
energy  of  system  1)  and  we  can  say 


8 K2Ao2 

'112  " if  P1C1A1.P2C2A2 


(3.4.10) 
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Note  the  similarity  between  this  result  and  that  of  Eq. 
(3.1.28). 


Two  rooms  coupled  by  a limp  curtain.  As  a final 
example  in  this  section/  we  consider  an  application  of 
Eq.  (3.3.22).  We  want  to  evaluate  ni2  for  the  system  shown 
in  Fig.  3.17.  Two  rooms#  of  volume  and  V2  are  separated 

by  a "wall"  that  has  a mass  density  m and  area  Aw.  The 
average  transmission  coefficient  for  this  situation  is  known 
to  be  approximately 


<T  (ft)  COS<J>>  ^ * 3 


(3.4.11) 


Thus 


n 


12 


p*ciAw 


(3.4.12a) 


The  energy  expression  Eq.  (3.2.14)  in  this  case  becomes  simply 


<P,>  V' 


<p{>  V 


n12 

W • 


(3.4.13) 


Where  n2  * 2.2/fT^,  T^  is  the  reverberation  time  of  room 
2#  and 

p2C3A 

n « 2 . (3.4.12b) 

21  6w3m2V2 

Equation  (3.4.13)  predicts  that  p2  * p|  when  losses  in  V2 

due  to  transmission  of  sound  through  the  wall  exceed  the 
losses  due  to  sound  absorption  in  room  2. 
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FIG  3.1 

TWO  LINEAR  RESONATORS  COUPLED  BY  SPRING,  MASS, 
AND  GYROSCOPIC  ELEMENTS 


FIG  3.2 

SPECTRAL  DENSITY  OF  INDIRECTLY  EXCITED  RESONATOR 
VELOCITY  AS  A FUNCTION  OF  THE  DEGREE  OF  COUPLING 
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(0)  (b)  (c) 


FI6  3.3 

( o ) SYSTEM  OF  INTEREST jTWO  SUBSYSTEMS  INDEPENDENTLY  EXCITED 
AND  RESPONDING  BOTH  TO  EXCITATION  AND  COUPLING  FORCES 
(b)  SYSTEM  1 VIBRATING  IN  RESPONSE  TO  ITS  OWN  EXCITATION, 
SYSTEM  2 IS  BLOCKED,  (c)  SYSTEM  2 IS  VIBRATING  IN  RESPONSE 
TO  ITS  EXCITATION,  SYSTEM  1 IS  BLOCKED. 


SYSTEM  1 SYSTEM  2 


(o)  BLOCKED 


FIG  3,4 

(a)  SHOWS  MODAL  SETS  WHEN  THERE  IS  NO  INTERACTION,  AS 
ACHIEVED  BY  CONDITIONS  SHOWN  IN  FIG  3.3  (b,c). 
lb)  SHOWS  MODE  PAIR  INTERACTIONS  THAT  OCCUR  WHEN  ACTUAL 
CONDITIONS  ARE  OBTAINED  AT  THE  JUNCTION 
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^12 

^2, tot 

« 1 I 

• 

N 

2 

1 1,  «iM 


’ 2.4i«« 


FIS  3.5 

AN  ENERGY  TRANSFER  /NO  STORAGE  MODEL  THAT  CAN  REPRESENT 
THE  SYSTEM  SHOWN  IN  FIG  3.4  (b)  WHEN  THE  SEA  MODEL  IS  USED 
THE  DEFINING  PARAMETERS  ARE  THE  MASSES  MitM2  ; 

THE  LOSS  FACTORS  if, , i*2  ; THE  MODAL  DENSITIES  iyN,/A« 
n£  *N2/A«  AND  THE  COUPLING  LOSS  FACTOR  ij12 


MODAL  RESPONSE 
IS  INCOHERENT; 
EQUALLY  ENERGETIC 


WAVES  ARE 
INCOHERENT; 
POWER  - i(\ 


ALLOWED 

WAVES 


(q)  MODAL  LATTICE 


(b)  WAVE  DISTRIBUTION 


FIG  3.6 


EQUIVALENCE  OF  MODAL  AND  WAVE  COHERENCE  ASSUMPTIONS  IN  SEA 
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IF  THE  SYSTEM  A IS  RECIPROCAL, THEN  IF  t AND  p'  ARE  'PRESCRIBED 
DROP’  SOURCES  AT  FREQUENCY  w,THEN  ONE  HAS  U/A  > v'/p' 


SYSTEM  B IS  THE  SAME  AS  SYSTEM  A WITH  THE  TWO  TERMINAL 
PAIRS  “OPEN"  THE  RECIPROCITY  STATEMENT  IN  THIS  CASE  IS  p/v  « Jt'/U' 


SYSTEM  C IS  THE  SAME  AS  SYSTEM  A WITH  TERMINAL  PAIR  I 
OPEN  - CIRCUITED  AND  TERMINAL  FAIR  2 SHORT-CIRCUITED.  THE 
RECIPROCITY  STATEMENT  IS  NOW  U/v«X>p'. 
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PLATE  DRIVEN  BY  POINT  LOAD  NOISE  SOURCE  WITH 
SECOND  POINT  ON  PLATE  EDGE  FIXED  WITHOUT  MOMENT  RESTRAINT 


THE  PLATE  IN  FIG.  3.10  IS  NOW  JOINED  TO 
ANOTHER  PLATE  AT  ONE  POINT  ALONG  THE  EDGE 
WITH  NO  MOMENT  COUPLING  OR  RESTRAINT 


ill 


FIG  3.12 

SKETCH  OF  2-PLATE  SYSTEM  OF  FIG  3.  II  SHOWING  HOW  THE 
ACTUAL  INTERACTION  IS  DEVELOPED  AS  A COMBINATION 
OF  2 IDEALIZED  PROBLEMS 


FIG  3.13 

ENERGY  FLOW  TOWARD  SYSTEM  JUNCTION  IN  THE  ANGULAR  INTERVAL  dfl 
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/ 


*«•**??*- 


\ 


FIG  3. 15 

A PISTON  RESONATOR  IN  THE  WALL  OF  A REVERBERATION 
CHAMBER  INTERACTING  WITH  THE  CONTAINED  SOUND  FIELD 
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FIG  3.16 


TWO  BEAMS  IN  LONGITUDINAL  VIBRATION 
COUPLED  BY  A SPRING  OF  STIFFNESS  K 


TWO  ROOMS  COUPLED  BY  LIMP  WALL  OF  SURFACE  DENSITY  m 


CHAPTER  4.  THE  ESTIMATION  OF  RESPONSE  IN  STATISTICAL  ENER6Y 
ANALYSIS 


4.0  Introduction 


In  the  preceding  chapters,  we  have  shown  how  the  energy 
of  vibration  may  be  estimated  in  an  average  sensed  based  on 
system  parameters  that  are  fundamentally  descriptive  of  energy 
storage  and  transfer  capabilities  of  the  systems.  The  use  of 
energy  quantities  has  the  great  advantage  that  sound,  vibration 
and  other  resonant  systems  all  are  described  by  the  same 
variables.  Consequently,  our  viewpoint  lias  great  generality 
which  is  useful  when  we  must  deal  with  systems  of  such  com** 
plexity. 

This  generality,  however,  has  its  shortcomings.  We  must 
remember  that:  our  motivation  for  these  studies  is  usually  not 
directed  to  the  energy  of  oscillation,  as  such.  Systems  do 
not  make  too  much  noise,  or  fatigue,  or  have  component  mal- 
function because  they  contain  too  much  energy,  but  because 
they  move,  are  strained  or  are  stressed  too  much  or  too  often. 
Consequently,  we  must  interpret  the  energy  of  oscillation  in 
terms  of  dynamical  variables  of  engineering  interest  such  as 
velocity,  strain,  pressure,  etc.  We  have  already  done  this  to 
a degree  in  earlier  chapters,  particularly  in  paragraphs  3.4. 

In  paragraph  4.1  we  shall  expand  on  this  important  aspect  of 
response  estimation. 

The  calculations  and  formulas  of  energy  prediction  in 
Chapter  3 are  directed  toward  predicting  average  energies, 
the  average  taken  over  the  set  of  systems  forming  the  ensemble 
of  "similar"  systems  constructed  according  to  the  SEA  model 
described  in  paragraph  3.2.  Since  we  have  a population  of 
systems,  the  energy  actually  realized  by  any  one  system  (the 
one  sitting  in  the  laboratory,  for  example)  will  not  be 
precisely  equal  to  the  average  energy  as  calculated  by  the  SEA 
formulas.  We  can  get  an  idea  of  how  much  our  system  may  deviate 
from  the  average  system  by  calculating  a standard  deviation 
(s.d.)  of  the  system  energy.  This  has,  in  fact,  been  done  for 
some  cases,  and  we  can  give  some  general  ideas  for  estimating 
the  variance  (square  of  the  s.d.)  for  practical  situations. 

The  calculation  of  response  variance  is  discussed  in 
paragraph  4.2. 
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The  calculation  of  standard  deviation  is  reassuring  to  one 
using  the  mean  energy  as  an  estimate  only  if  the  s.d.  is  a 
small  friction  of  the  mean.  Then  we  know  that  most  realizations 
of  the  system  (again  including  the  one  sitting  in  the 
laboratory)  will  have  a response  that  is  close  to  the  mean. 

But  what  if  the  standard  deviation  is  equal  to  the  mean,  or 
greater?  In  this  event,  the  probability  that  any  one  real- 
ization is  close  to  the  mean  will  be  small  and  the  mean  value 
loses  its  worth  as  an  estimate. 

In  such  a situation,  it  is  common  to  calculate  the 
probability  that  a realization  will  occur  within  some  stated 
interval  Of  response  amplitudes.  The  probability  associated 
with  this  "confidence  interval"  is  called  the  "confide  ice 
coefficient".  In  order  to  carry  out  these  calculations,  how- 
ever, we  need  a probability  distribution  for  the  energy.  This 
distribution  is  introduced  and  the  estimation  procedures  are 
derived  from  it  in  paragraph  4.3. 

The  distributions  used  for  generating  confidence  intervals 
tend  to  be  fairly  accurate  of  the  "middle  range"  of  response 
levels,  but  are  generally  fairly  inaccurate  in  predicting 
response  having  a very  small  probability.  If  we  have  a special 
interest  in  either  the  greatest  or  smallest  response  levels 
that  may  occur,  we  usually  must  take  a different  approach. 
Paragraph  4.4  describes  the  estimation  of  the  probability  of 
very  large  values  (having  small  probability)  in  a particular 
situation.  It  may  be  argued  that  the  problem  treated  in 
paragraph  4.4  is  not  properly  an  SEA  problem,  because  it 
deals  with  the  coherent  properties  of  modal  response.  Never- 
theless, the  entire  approach  of  describing  modal  responses  by 
their  m.s.  values  and  adopting  a statistical  view  of  the 
population  of  systems  being  considered  is  very  much  within 
the  SEA  framework. 


4. 1 Mean  Value  Estimates  of  Dynamical  Response 


We  can  make  useful  estimates  of  the  mean  square  amplitude 
of  response  based  on  the  energy  calculations  in  Chapter  3. 

We  first  consider  systems  like  those  discussed  in  paragraphs 
2.3  and  3.4,  in  which  only  one  mode  of  one  system  is  inter- 
acting with  a group  of  modes  of  another.  This  one  mode  may 
actually  be  a single  dof  system  like  that  shown  in  Fig.  2.1 
or  it  may  be  a mode  that  we  single  out  for  particular 
examination  within  a multi-modal  system.  We  may  do  this  by 
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restricting  the  excitation  bandwidth  so  that  this  is  the  only 
system  mode  that  is  resonant. 


Single  Mode  Response 


Suppose  that  the  multimodal  system  1 is  excited  by  noise 
sources  of  bandwidth  Aw.  System  1 has  total  energy  Ex  tot 
and  Nx  inodes  that  resonate  in  this  bandwidth.  System  z has 
only  one  mode  with  mode  shape  According  to  paragraph 

3.2/  the  energy  of  this  mode  will  be  given  by 


n21 

n2+n21 


(4.1.1) 


where  n21  may  be  determined  from  averages  over  interaction 
parameters  as  in  Eq.  (3.2.8)  or  from  boundary  impedances  as 
in  paragraph  3.3.  Let  us  assume  that  this  has  been  done,  and 
we  now  want  to  find  the  actual  response  at  location  X2. 

From  Eq.  (2.2.3)  we  have 


y2(x,t)  = Y2(t)  ^2^  ' 


(4.1.2) 


with  the  general  result  that 

E 2m  f dx  p<y2>t  - M2<Y^(t)>fc  = M2w|  <Y|(t)>t 

(4.1.3) 

where  Mp  is  the  mass  of  system  2 and  wp  is  its  resonance  fre- 
quency (which  might  just  as  well  be  taken  as  w,  since  u>2  is 
within  Aw  by  hypothesis).  Thus 
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<y|(x/t)>t  = —r l—  $2  f*> 


, El,tot  ^21 
N1  n2+n21 


(4.1.4) 


Thus,  our  statistical  model  of  the  system  does  not 
restrict  us  in  any  way  from  determining  the  spatial  dis- 
tribution of  response.  Of  course,  the  normalization  of  the 
<J>2  function  still  ensures  that 


<y|> 


P,t 


£ 2 / w Z^2' 


In  Fig.  4.1  we  show  the  form  of  the  distribution  of  y\  for 
a typical  mode  shape  of  a beam.  The  distribution  of  *<yjj>2 
which  would  be  proportional  to  the  surface  strain  in  the 
beam,  is  also  sketched. 


Multi-Modal  Response.  When  there  are  several  resonant  modes 
(N2T  of  (indirectly  excited)  system  2 in  the  band  Aw,  then 
Eq.  (4.1.1)  has  the  form 


E2, tot  _ _ El,tot  n21 

N2  " 2 " N1  V^l 


(4.1.5) 


and  the  m.s.  response  of  the  system  in  time  and  space  (or 
mass)  is 


<y«> 


2 P,t 


= E 


2 , tot 


/w*M. 


(4.1.6) 
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which  we  may  regard  as  our  estimate  of  the  system  response. 
Note,  however,  that  if  such  an  estimate  were  made  in  the 
situation  of  Fig.  4.1,  it  would  he  rather  poor,  since  the 
actual  value  of  <y2>t  only  equals  the  estimate  <y2>p  t at  a 
finite  number  of  points  and  oscillates  around  the  estimate 
without  converging  on  it. 

In  the  case  of  two  or  three  dimensional  systems  with  many 
modes  of  vibration,  the  function  tends  to  converge  on 

the  estimate  of  Eq.  (4.1.6)  in  many  important  instances.  To 
demonstrate  this,  we  consider  a rectangular  supported  plate 
which  has  modes 


ip  (x)  «*  2 sin^  sink2  x2 


(4.1.7) 


where 


k1  = dIi(/i1,k2ssa27r/i2jai  °2  9 


The  modes  that  resonate  within  the  interval  Ak  * Au)/c„  are 
shown  in  the  hatched  region  of  Fig.  4.2. 

As  a result  of  our  model  of  the  system  interaction, 
spelled  out  in  paragraph  3.2,  we  can  make  the  following 
observations  about  the  response  of  system  2. 


(1)  The  response  may  be  written  in  the  form 


y2(x,t)  * l Ya(t)  ifig  (x) 

o 


(4.1.8) 


(2)  The  modal  response  amplitudes  Yg  (t)  are  incoherent; 


<Ya(t)Yr(t)>-  <Y*(t)>t«0#T 


(4.1.9) 
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(3)  The  inodes  of  system  2 are  equally  energetic 


<Ya(t)>t  * <Y2 (t)  >t  * E2/bi2H2 


(4.1.10) 


Thus,  we  can  say  that 


<Y*  (x,t)>t»  Z <*a(t)yVr  (t)>t^CT  (x)  Vx) 
a t 


<y2(t)>t  1 * ' 

0 0 


(4.1.11) 


where  the  allowed  values  of  o are  determined  by  the  modes  of 
system  2 that  resonate  in  Aw. 

The  result  of  Eq.  (4.1.11)  is  quite  general.  For  the 
specific  two-dimensional  system  at  hand,  we  can  write, 
therefore , 


<y2>t/<y^>,_  p = 4N2*  Zsin2k1x1  sin2k2x2 


■ N2l  l (l-cos2k1x1) (l-cos2k2x2) 


(4.1.12) 


We  evaluate  the  sum  by  replacing  it  with  an  integral  over  <p 
in  the  hatched  region  of  Fig.  4.2.  We  can  see  that  if  dN2 
is  the  number  of  modes  in  the  rectangle  bounded  by  Ak  and  d$, 
then  dN-j/N2  ■ 2d<p/n,  since  the  modes  are  uniformly  distributed 
in  <p.  We  may  then  write 
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<y2>t/<y2>t,p "f f q^2  i-cos (2Joc1*in^  )cos  (2kx2cos$) 


+ cos  (2kx1sin^)cos  (2kx2cos^)l  d$ 


* 1-J0(2kx1)-J0(2kx2)  + J0(2kr) 

(4.1.13) 


where  r » yjx |+x|  . We  have  evaluated  the  integral  for  the 
region  near  x^,  X2  = 0.  Since  the  sum  in  Eg.  (4.1.12)  is 
unaffected  by  the  substitution  -x. , x2*ll'2m,x2,  the  tota* 

pattern  is  as  shown  in  Fig.  4.3. 

Thus,  we  can  estimate  a spatially  varying  temporal  mean 
square  response  of  a system  even  with  rather  extensive 
assumptions  regarding  the  system  interactions.  Patterns  of 
the  kind  shown  in  Fig.  4.3  have  been  developed  for  a variety 
of  systems,  and  are  available  from  the  literature.  They  are 
usually  of  interest  if  one  wishes  to  avoid  locating  an  item 
or  making  measurements  at  a point  where  the  response  is 
excessively  high  or  otherwise  unrepresentative. 

Wave  estimates.  The  use  of  wave  models  for  SEA  response 
estimates  extends  the  power  of  SEA  greatly.  We  can  develop 
response  patterns  like  the  one  shown  in  Fig.  4.3  equally  well 
using  a diffuse  wave  field  model  like  that  discussed  in 
paragraphs  2.4  and  3.3.  The  more  important  application  of 
wave  notions  employs  the  use  of  average  impedance  functions 
at  system  boundaries  for  the  evaluation  of  couplirig  loss 
factors  as  discussed  in  paragraph  3.3. 

As  an  example,  consider  the  system  shown  in  Fig.  4.4a, 
a beam  cantilevered  to  a flexible  plate.  Along  the  line  of 
connection  of  these  systems  there  will  be  very  little  dis- 
placement (we  shall  assume  there  is  none)  because  of  the  very 
high  in-plane  impedances  of  these  systems.  Nevertheless, 
power  will  be  transmitted  by  torques  (moments)  and  rotational 
motion  of  the  junction.  In  order  to  evaluate  the  coupling 
loss  factor,  we  must  find  the  power  transferred  between  the 
systems . 
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The  "decoupled"  systems  are  produced  by  eliminating  the 
rotation  of  the  contact  line  as  shown  in  Fig.  4.4b.  The  modes 
of  the  beam  for  this  condition  are  readily  found  - the  beam 
has  clamped  - free  boundary  conditions.  The  plate,  however, 
has  a very  complicated  boundary  condition  - the  boundary 
conditions  around  its  ouv  itr  edges  and  clamped  along  a finite 
line  in  the  interior!  It  would  require  a complicated  computer 
routine  to  find  its  modes. 

Nevertheless,  we  can  invoke  two  principles  from  our  work 
thus  far  that  allow  us  to  solve  this  problem.  The  first  is 
that  of  reciprocity  of  the  coupling  loss  factor 


Vpb 


Vbp 


(4.1.14) 


which  allows  us  to  use  n^p*  the  beam  to  plate  coupling  loss 
factor  (which  is  easier  to  calculate)  to  evaluate  npb*  Secondly, 
the  moment  impedance  "looking  into"  the  plate  may  b€  evaluated 
by  considering  the  plate  to  be  infinite.  This  has  the  effect 
of  eliminating  reflected  energy  in  the  plate  from  returning 
to  the  junction.  In  a manner  completely  analogous  to  the 
results  in  paragraphs  2.2  and  2.3  for  the  admittance  for  a 
transverse  force,  the  moment  impedance,  averaged  over  the 
location  of  the  junction  of  the  plate  and  over  resonance 
frequencies  of  the  interacting  systems  in  the  band  Au>,  is  given 
by  the  infinite  system  impedance,  which  is  known. 

In  paragraph  3.3,  we  obtained  a general  result  that,  in 
terms  of  actual  boundary  impedances,  the  power  flow  between 
systems  is  given  by 


^bp  Eb,tot 


(4.1.15) 


where 


u a.  , 
hp 


2 ¥p 

nnb  jzb+zpl * 


(4.1.16) 
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On  the  other  hand,  we  have  the  general  result  that  (see 
Eg.  3.2.10) 


Nb 

nbp  * wnbp^Eb,tot  " N~  Ep , tot^ 


(4.1.17) 


It  must  be  emphasized  that  Eqs.  (4.1.15)  and  (4.1.17)  are 
fully  equivalent  and  apply  to  all  conditions  of  system  size, 
damping,  etc.  In  general  and  aup  are  different  functions. 
However,  if  the  receiving  system  (the  plate  in  this  instance) 
is  made  very  large,  while  its  loss  factor  rip  is  held  fixed, 
the  term  Nb  Ep  tot/Np  will  vanish  in  Eq.  (4.1.17).  If  we  use 
the  infinite  plate  impedances  in  Eq.  (4.1.16), 


2 W“* 

&bp  * nwnb  |zb+2p<«.)  | = 


(4.1.18) 


The  relation  between  n^p  and  ab  is  then  readily  found. 
Since  the  input  power  to  the'system  pis  independent  of  the 
junction, 


ni,in  = *lhXot  A1  = El,totAl+wa12  El^tot 


(4.1.19) 


we  have 


“"bp  Eb,tot  ' “''bp 


1 - 

Ab 


?(b) 

Jb,  tot 


= u)a,  eP^]_  , 

bp  b, tot 


(4.1.20) 


and,  consequently, 


nbP  = abP  (1  • “abp/Ab>'‘ 


(4.1.21) 
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When  we  deal  with  a single  beam  mode,  then  Ab=  Rfc/Mb  and 
the  average  of  abp  over  the  band  is 


< 


abp>Aw 


w> 
<v  v 


vi. 

«b 


(4.1.22) 


Substituting  into  Eg.  (4.1.21)  then  gives 


wnbp  * Rp(0#)/Mb 


(4.1.23) 


On  the  other  hand,  when  the  modes  of  the  beam  are  dense  so 
that  n^A^  »1,  the  beam  "looks  infinite",  and  one  has 


M»bp  , 2 Rb‘-)RP<“>  1 << 

~o~  " |2b(«)+Zp(")  T Vb 


and  referring  to  Eq.  (4.1.21),  we  get 


(4.1.24) 


nbp  + abp  ” irn^ 


V^V00) 

|2b(«)+zp(fl<»)| 


(4.1.25) 


Where  we  have  used  Zp(»)*Rp(«)-i  Xp(«)  to  emphasize  that 

these  are  input  impedance  functions  for  the  infinitely 
extended  system.  Extending  the  receiving 
the  value  of  abD,  but  does  not  affect  the 
since  nbp  depends  on  a mode-mode  coupling 

<Ba0>  in  the  form  [see  Eq.  (3.2.10)] 


system  affects 
value  of  tv  , 
coefficient 
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(4.1.26) 


“”bp  ’ <Bbp>  "p  4“ 


As  the  plate  gets  larger,  the  coupling  <Bbp>  diminishes 
inversely  with  plate  area  since  the  quantities  *m2,y2»k2 
are  inversely  proportional  to  the  mass  of  the  receiving 
system  according  to  Eq.  (3.1.4)  (or  plate  area)  and  np  is 
proportional  to  plate  area.  Thus,  while  Eq.  (4.1.16)*applies 
to  any  situation,  Eq.  (4.1.21)  is  only  correct  when  infinite 
system  parameters  are  used  for  the  receiving  system. 

In  the  system  of  Fig.  4.4,  the  infinite  system  junction 
impedances  are 


Zb  ' pbCbSbKb°fb(1+i)  (4.1.27) 


and 

Z*J  - 16p  K2cVu)(l-ir)  (4.1.28) 

P s P P 


where  p^,  Cfc,^  are  the  material  density,  longitudinal  wave- 
speed,  and  radius  of  gyration  for  the  beam  with  the  same  para- 
meters for  the  plate  with  a "p"  subscript.  is  the  cross- 
sectional  area  of  the  beam,  c^  = ^w<bcb  is  flexural 
phase  speed  on  the  beam  and  ps  the  mass  per  unit  area  of  the 
plate.  The  quantity  r is  a moment  susceptance  parameter 
that  depends  on  the  shape  of  the  junction.  If  we  assume  that 

Zp  >>  Zt,  then  Rp|Zpj  2 +Gp  , the  plate  moment  conductance 
and  the  result  is  independent  of  r.  With  this  assumption 
and  assuming  the  beam  and  plate  are  cut  from  the  same 
material,  we  get 


v * v/u' 


(4.1.29) 
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a surprisingly  simple  result.  We  have  also  used  the  modal 
density  of  the  beam 


nb(w)  * i/2ircfb 


(4.1.30) 


The  coupling  loss  factor  from  plate  to  beam  is,  according  to 
Eq.  (3.2.12) 


npb  ■ Vbp/np  ' (4.1.31) 


where  * A /4tmc  c . 

P P P P 

Using  Eq.  (4.1.29),  we  obtain  an  estimate  for  the  average 
response  of  the  beam  for  bands  of  noise, 


<vb> 


<v2> 


Jsb-.. 

W 


(4.1.32) 


In  Fig.  4.5  we  show  a comparison  of  this  simple  theoretical 
result  for  the  two  limiting  experimental  conditions,  first 
for  and  then  for  nbp< <r>b  • Modal  energy  equipartition 

between  the  beam  and  plate  would  be  expected  to  obtain  in  the 
first  instance,  but  not  in  the  second. 

Strain  Response.  As  an  illustration  of  how  energy 
response  may  be  expressed  in  strain  or  stress  of  the  system, 
we  have  already  shown  that  the  pressure  within  a sound  field 
is  given  by 


J— ^ * <v2>  M (4.1.33) 

po°2 
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where  <v2>  is  the  m.s.  velocity  of  the  fluid  (air)  particles. 
Using  PqC2  * yPQf  this  expression  may  also  be  written  as 


<P.22_  * <(££)%  * — 
(YP0) 2 p0  c2 


(4.1.34) 


where  6p/po  is  the  volumetric  strain  (dilatation)  of  the 
fluid  since  yPQ  is  the  volumetric  stiffness  (bulk  modulus) 
of  an  ideal  gas.  The  equation  states  that  the  m.s.  strain  equals 
the  m.s.  mach  number  (ratio  of  particle  velocity  to  sound 
speed)  of  the  particles. 

Quite  a similar  result  also  obtains  for  plates.  In  a 
plate  of  thickness  h,  the  strain  distribution  across  the 
plate  thickness  is  given  by  (see  Fig.  4.6) 


e(z)  “IT  em  (4.1.35) 


where  is  the  strain  at  the  free  surface  of  the  plate, 

where  it  is  a maximum.  The  strain  energy  density  of  a layer 
of  plate  material  dz  thickness  is  1/2  dze2(z)  YQ,  where  Yfl 
is  Young's  modulus,  so  that  the  total  strain  energy  density 
is 


2k2 

h 


v*£> 


(4.1.36) 


If  this  is  equated  to  the  kinetic 
we  get 


energy  density 


h<v2> 


9 
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<e1 2>  “ — — *3  - (homogeneous)  (4.1.37a) 

*k  ci  °i 

where  we  have  used  <2=h2/12  for  a homogeneous  plate.  For  a 
sandwich  plate  with  all  its  stiffness  at  the  surface, 
x=h/2  and  one  has 


<ei>  * <v2>/c2  (sandwich).  (4.1.37b) 

m *. 

The  results  of  Eq.  (4.1.37)  are  entirely  analogous  to  that  of 
Eq.  (4.1.34)  and  allow  us  to  develop  mean  square  strain  estimates 
from  energy  (or  velocity)  estimates. 


4.2  Calculation  of  Variance  in  Temporal  Mean  Square  Response 


In  paragraphs  3.3,  3.4  and  4.1  we  have  shown  how  average 
energy  estimates  of  system  response  may  be  developed  and 
how  m.s.  response  estimates  are  derived  from  the  energy.  We 
have  seen,  however,  in  Figs.  4.1,  4.3  and  4.5  that  the  average 
energy  may  not  be  a good  estimate  of  m.s.  response  at  a 
particular  location  on  the  system  or  in  a frequency  band  in  which 
the  number  of  modes  is  too  low.  In  this  section,  we  shall 
develop  estimates  for  the  deviation  of  realized  response  from 
the  mean  in  terms  of  the  standard  deviation  (square  root  of 
variance)  from  the  mean  square  response.  Note  that  the 
variance  that  we  are  discussing  is  from  one  "similar"  system 
to  another,  or  from  one  location  to  another,  not  temporal 
variance.  If  we  are  dealing  with  noise  signals,  we  assume 
that  the  averaging  time  is  sufficiently  long  to  remove 
variations  due  to  this  cause. 

The  variance  in  response  between  the  mean  and  any  member 
of  the  population  is  produced  by  four  major  effects. 

(1)  The  modal  energies  of  the  directly  excited  system 

may  not  be  equal,  because  of  the  spatial 
dependence  of  the  actual  sources  of  excitation 
and  the  internal  coupling  of  the  blocked  system 
may  not  be  great  enough  to  ensure  equipartition. 
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(2)  The  actual  number  of  resonant  interacting  modes 
in  each  realization  of  the  coupled  system  will 
vary. 

(3)  The  strength  of  the  coupling  parameters  will 
fluctuate  from  one  realization  to  another  because 
of  slight  changes  in  modal  shape  at  the  junction 
of  the  systems. 

(4)  The  response  at  the  selected  observation  position 
will  vary  because  the  observation  position  is 

, randomly  located  and  because  the  mode  shapes  change 
from  one  realization  to  another. 

The  way  in  which  the  variance  of  re  ponse  depends  on 
these  four  factors  is  not  obvious,  and  indeed,  has  not  been 
worked  out  for  many  cases  of  interest.  There  is  a con- 
siderable  area  of  interesting  research  work  thit  needs  to 
be  done  in  analyzing  variance  of  interacting  systems.  We 
should  also  remark  here  that  we  are  required  to  use  modal 
analyses  to  calculate  variance.  The  wave-mode  duality  is 
useful  for  mean  value  estimates,  but  a wave  analysis  of 
variance  by  its  nature  disregards  spatial  coherence  effects 
that  are  essential  to  the  calculation  of  variance. 

Modal  Power  Flow  and  Response.  To  calculate  variance 
in  the  modal  energies,  we  return  to  the  system  of  Fig.  3.4. 
The  power  flow  from  all  of  the  resonators  of  system  1 to 
resonator  a is  (system  2 has  no  external  excitation) 


"lo  = ZA 

a ao  a 


(b) 


ZB 

a a a 


<vv 


(4.2.1) 


where  A~a  is  given  by  Eq.  (3.1.20)  and  Ba0  is  given  by 
Eq.  (3.1.26)  with  appropriate  subscripts  entered.  If  the 
damping  of  the  "o"  is  h0,  then  the  dissipation  of  energy  in 

this  mode  is  and  we  have 

o o 

Ea  SEib‘  'Ac</V  - Wa„>  <VSBac>*‘-  <4-2-2’ 
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If  ve  are  interested  in  the  velocity  response,  then  we  saw 
from  Eq.  (4.1.10)  that 


<va>t 


- <VM2,«'§(x2) 


(4.2.3) 


where  x2  is  the  observation  position.  Thus,  the  total  m.s. 
response  of  system  2 observed  at  position  x2  is  just 


<V2^X2,t)>t  “ M Za, 


<x2>  i E!a)Aaa/A0 


(4.2.3a) 


* k;  lfc  <x2>(«.Bao  ) /fa+fao)  (4-2-3b> 


In  writing  Eq.  (4.2.3)  we  have  now  relaxed  the  assumption 
of  Chapter  3 that  the  modes  of  system  1 are  equally  energetic, 
but  we  continue  to  assume  that  modal  responses  are  temporally 
incoherent. 

The  four  factors  listed  above  that  contribute  to  variance 
are  specifically  shown  in  Eq.  (4.2. 3) • The  variance  in  modal 
energies  (1)  is  in  E^b)  or  Ea.  The  coupling  variance  strength 
(3)  is  represented  by  the  quantities  Aaa  or  Baa,  and  the 
observation  variance  (4)  is  incorporated  in  \p^.  The  variation 
in  number  of  interacting  modes  (2)  is  represented  by  the 
summation  in  a.  Our  problem  is  to  find  the  variance  in  this 
sum  on  a knowledge  (or  reasonable  assumptions)  about  the 
statistics  of  the  terms  entering  the  equation. 

To  simplify  matters  a bit,  let  us  first  concentrate  on 
the  single  mode  response  <v*>t  so  that  only  a single  sum  over 

a is  involved.  The  appearance  of  B^  in  the  denominator  of 
Eq.  (4.2.3b)  is  awkward  from  an  analytical  point  of  view,  so 
that  we  would  prefer  to  use  Eq.  (4.2.3a).  However,  since  the 
coupling  factors  also  enter  the  denominator  of  Aa0  (tut  not  of 
Baa) , we  cannot  avoid  this  problem.  Another  reason  for  pre- 
ferring Eq.  (4.2.3a)  is  that  we  might  be  able  to  estimate  the 

statistics  of  E^  better  than  Ea,  since  the  former  does  not 
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contain  contamination  by  coupling  effects.  Nevertheless, 
we  are  essentially  forced  to  use  Eq.  (4.2.3b)  since  it  alone 
allows  us  to  use  an  important  result  in  the  statistics  of 
summed  random  processes  that  is  discussed  in  the  following 
paragraphs . 

There  are  two  situations  in  which  Eq.  (4.2.3b)  can 
simplify  sufficiently  to  allow  one  to  make  calculations. 

The  first  is  where  |B  >>A  . In  thi®  event,  we  have 


<vj>t  - <x2)/M2)  SB^SB  . (4.2.4) 

a a 


In  this  case,  the  energy  of  the  q mode  is  a weighted  average 
of  the  energy  of  the  modes  of  system  1 with  weighting  co- 
efficients baa  * Bda/5Baa  • Notice  that  this  condition,  which 
is  essentially  an  equipartition  condition,  does  not  require 
that  Ba0»Ao,  and  is  thus  a much  weaker  requirement  on  the 
coupling  than  when  one  has  two  single  dof  systems  interacting. 

The  second  situation  for  which  one  can  calculate  variance 
is  when  Aa>>5Bcta'  In  th-s  case 


<v0>t  - <*o/M2V  iVoo  «-2-5> 


We  note  that  Baa  is  given  by  Eq.  (3.1.26)  with  appropriate 
subscripts.  We  will  write  it  here  in  the  form 
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(4.2.6) 


where  Aaa  is  given  by  Eq.  (3.1.33)  and  £aa=2  (u)0-u>a)/ (Aa+Aa)  , 
in  a manner  analogous  to  the  development  in  paragraph  2.2. 

If  we  now  assume  A^A^®  const  (all  modes  of  system  1 have 
equal  damping),  then  according  to  Eq.  (2.2.21), 
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(4.2.7) 


<B  >w 
a a a 


AlAo 

<Xac>a,o  TEu 

1 a 


which  reproduces  Eq.  (3.2.8) 

In  the  present  instance,  we  are  interested  in  the 
variance  of  the  sum  in  Eq.  (4.2.5).  The  assumption  Aa«Ai  has 
the  effect  of  causing  each  of  the  functions  (^(J+I)“1  to 
have  the  same  shape  as  a function  of  frequency.  We  can  plot 
Eq.  (4.2.5)  then  as  shown  in  Fig.  4.7.  It  appears  as  a sum 
of  pulses,  of  strength 
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(4.2.8) 


If  we  assume  that  the  resonance  frequencies  are  uniformly 
probable  over  the  interval  Aw  and  are  independent  of  each 
other,  then  the  standard  deviation  of  Eq.  (4.2.5)  may  be 
shown  to  be 


a\=5nl(u)  (W 
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This  is  the  "important  result"  referred  to  above.  Consequently, 
using  Eq.  (4.2.7)  in  the  average  of  Eq.  (4.2.5),  we  get 
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showing  the  explicit  dependence  of  variance  on  uncertainty  in 
the  modal  energy  of  the  directly  excited  system  and  in  the 
coupling  parameters.  This  relation  also  shows  the  effect  of 
the  number  of  overlapping  modes  in  reducing  the  ratio  of 
variance  to  the  square  of  the  mean. 

Finally , if  we  refer  to  Eq.  (3.2.10)  we  see  that  the 
quantity  Ni<Bgto>ot  i wnoi*  The  standard  deviation  in  the  sum 
£Ba0  will  tell  us  something  about  the  uncertainty  in  the 
coupling  loss  factor  r?2i*  If  for  the  present  discussion  only, 
we  define  WTI21  * IBQ0  then  m^  * <Ti21>a  *s  the  coupling  loss 

factor  we  have  bein  dealing  with.  Its  variance  is  then  given 
by 
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(4.2.11) 


Returning  to  Eq.  (4.2.5),  if  the  observation  position 
X2  is  located  randomly  in  space,  then  we  have  additional 

factor  of  <\p*>  /<^2>2  in  Eq.  (4.2.10)  since  must  be 
regarded  as  a random  amplitude  factor.  If  we  assume  that 
system  1 is  excited  by  a randomly  located  point  force  at  x^, 
then  according  to  Eq.  (2.2.19),  the  energy  of  the  «th  mode 
will  vary  as  t|»2(x^).  If  we  assume  that  the  systems  are  joined 
at  a point  xj  “in  system  1 coordinates  and  in  system  2 
coordinates,  then  according  to  Eqs.  (3.2.6)  and  (3.2.9),  we 
may  expect  Xaa  to  vary  as  ^2  (xj)  ^J(xJ)  . Thus,  the 
greatest  variation  that  we  can  expect  in  the  response  of  a 
single  mode  of  system  2 is 
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(A1+A2) 


2 


(4.2.12) 


where  i/'i  and  4>2  are  typical  node  shapes  of  system  1 and 
system  2 respectively.  We  may  expect  that  in  any  particular 
system,  the  variance  will  be  less  than  that  implied  by 
Eg.  (4.2.12). 


Now  suppose  that  we  have  a group  of  modes  of  system  2 
represented  by  a modal  density  n2  (w) . We  have  a combination 
n2Au)  independent  response  functions,  each  of  which  has  a 
ratio  of  variance  to  square  of  mean  given  by  Eg.  (4.2.12). 
The  total  ratio  of  variance  to  square  of  mean  for  multi- 
modal response  of  system  2 is,  therefore, 
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(4.2.13) 


Notice  that  this  result  is  completely  symmetric  in  the  pro- 
perties of  system  1 and  2.  We  should  expect,  therefore,  that 
for  systems  excited  and  observed  at  point  locations,  the 
variance  in  response  will  be  independent  of  which  system  is 
used  as  a "source"  and  which  is  used  as  a "receiver". 


It  will  be  apparent  to  the  reader  that  a substantial 
number  of  assumptions  have  been  made  in  order  to  get  the 
result  expressed  by  Eq.  (4.2.13).  Most  of  these  assumptions, 
with  notable  exception  of  the  assumptions  of  uniform  damping, 
Aa=Aj,  and  Aa=A2,  would  have  the  effect  of  reducing  response 
variance.  We,  therefore,  may  reasonably  regard  Eq.  (4.2.13) 
as  a conservative  estimate  of  the  variance  in  the  work 
that  follows  in  that  it  estimates  more  variance  than  may,  in 
fact,  occur. 
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4.3  Calculation  of  Confidence  Coefficients 


We  see  from  Eq.  (4.2.13)  that  increasing  the  number  of 
modal  interactions  by  increasing  the  combined  modal  bandwidth 
ir/2  (Aj+A2)  , and  increasing  the  number  of  modes  observed  by 

increasing  the  excitation  bandwidth  Au  will  both  cause  a 
reduction  in  variance  of  response.  In  this  section,  we  show 
how  the  variance  estimates  may  be  used  to  predict  the 
fraction  of  measurements  that  will  fall  within  a certain 
interval  of  values,  generally  based  on  the  predicted  mean. 

The  Probability  Density.  If  the  total  probability  den- 
sity <f»  (0 ) for  the  observed  response  <v2>t=0  were  known , then 

the  probability  that  0 would  fall  within  the  "confidence 
interval"  01<0<02is  simply 


where  the  (fiducial)  probability  . CC  is  called  the  confidence 
coefficient.  Since  we  do  not  know  what  the  function  <f>  (0 ) is, 
this  may  seem  rather  useless.  We  do  know,  however,  that  0 
has  positive  values  only,  and  we  can  estimate  its  mean  mQ 
and  its  standard  deviation  a ^ from  paragraph  4.2. 

A probability  density  that  meets  our  requirements  and 
has  tabulated  integrals  is  the  "gamma  density" 


<j>(0)  = e^exp  <-0/A)/Awr  (y) 


(4.3.2) 


where  y =m2/c2  and  A=o2/me  and  T(y) 

If  we  change  variables  to  y = 0me/o2  , 
becomes 


is  the  gamma  function, 
then  Eq.  (4.3.1) 
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where  y(\i,B)  is  the  incomplete  gamma  function  defined  by 


Y(M/B) 
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(4.3.4) 


and  is  tabulated  in  standard  mathematical  handbooks. 

Confidence  Intervals.  Firsts  let  us  consider  a simple 
" exceedance ■ type  of  confidence  interval.  That  is,  search  for 
a value  of  6,,^  5 rm0  such  that  the  probability  (CC)  that  any 
realized  value  of  6 < 0max  is  known.  The  equation  to  be 
solved  then  is 


_i 

CC  - r (p)  Y (p,rp) 


(4.3.5) 


and  the  result  is  a line  of  constant  CC  as  a function  of  r 
and  p.  The  solution  to  Eq.  (4.3.5)  as  found  numerically  is 
graphed  in  Fig.  4.8. 

From  Fig.  4.8  we  see  that  if  the  ratio  of  the  standard 
deviation  to  the  mean  (as  computed  from  Eq.  (4.2.13)  for 
example)  is  equal  to  unity,  then  we  may  expect  a measured 
response,  selected  at  random,  to  have  a value  less  than  2.5dB 
more  than  the  mean  in  80%  of  the  cases.  The  measurement  should 
be  less  than  5 dB  more  than  the  mean  in  95%  of  the  cases  and 
less  than  7 dB  more  than  the  mean  in  99%  of  the  cases.  In 
many  situations,  this  is  precisely  the  way  that  we  would  want 
the  estimate  stated,  particularly  if  we  were  concerned  about 
the  response  exceeding  some  prescribed  value  that  might  lead 
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to  component  failure/  or  excessive  speech  interference,  or 
acoustical  detectability. 

Also  note,  that  since  the  calculation  of  a§/m|  is  to  be 
carried  out  in  frequency  bands  A u>,  or  for  individual  modes 
that  resonate  at  different  frequencies,  the  form  of  the  con- 
fidence interval  is  that  of  a spectrum.  One  would  usually 
plot  an  "average  estimate",  based  upon  the  methods  of 
section  4.1  as  in  Fig.  4.5,  then  plot  a second  spectrum 
determined  by  10  log  r for  each  frequency  band.  One  would 
then  look  upon  the  second  spectrum  as  an  "upper  bound"  for 
the  data,  with  a degree  of  confidence  determined  by  the 
selected  confidence  coefficient. 

As  a second  example  of  a confidence  interval,  consider 
one  that  brackets  the  calculated  mean.  Such  interval  might 
be  stated  as  "in  the  500  Hz  octave  band,  I estimate  that  in 
95%  of  the  cases,  the  measured  acceleration  level  is  within 
±5  dB  of  the  average  estimate  of  -10  dB  re  lg."  The  "cases" 
referred  to  here,  of  course,  are  reported  measurements  with 
the  loading  and  observation  points  varied  consistent  with 
the  calculations  of  mean  and  variance,  and  for  a population 
of  systems.  Obviously,  one  rarely  has  this  population  to 
deal  with;  there  is  only  one  system  in  the  laboratory,  or  at 
most  2 or  3.  The  system  that  one  has  might  be  very  good  or 
very  poor  from  a viewpoint  of  variance,  but  we  must  hope 
that  the  estimate  for  the  population  will  be  useful  in  deal- 
ing with  it.  This  situation  is  directly  analogous  to  gambl- 
ing. One  only  has  a particular  realization  of  all  possible 
hands  of  cards  before  him,  but  the  pattern  of  betting,  cal- 
culation of  odds  and  strategy  is  made  on  the  basis  of  all 
possible  hands  and  draws. 

The  evaluation  of  this  "bracketing"  interval  is  made  by 
setting  62*1^10  and  8i=me/r  in  Eq.  (4.3.3).  The  equation  to 
be  solved  in  this  case  is 


CC  * r (p)  (y  (y,M  r)-y(p,p/rU. 


The  solution  to  Eq.  (4.3.6)  is  shown  in  Fig.  4.9.  For 
example,  the  value  of  normalized  variance  that  would  produce 
the  hypothetical  estimation  interval  of  the  preceding 


paragraph  is  a2/ m2  = 0.2.  An  interval  of  this  type  is  very 
useful  in  determining  the  value  of  the  mean  alone  as  an 
estimate  of  response.  Since  one  can  usually  only  make  measure- 
ments of  structural  vibration  and  acoustical  noise  to  a 
repeatability  of  1 or  2 d8,  then  if  an  estimation  interval  of 
± 1 dB  has  a high  confidence  coefficient  (80%,  say)  then  the 
mean  value  is  a very  good  estimate.  We  see  from  Fig.  4.9 
that  this  requires  that  a2/m2<0.1. 

Example.  As  an  illustration  of  the  discussion  in 
paragraph  4 . 2 and  4.3,  we  consider  the  system  shown  in 
Fig.  4.4  excited  at  a point  and  observed  at  a point.  We 
derived  an  average  energy  estimated  for  this  system  in 
paragraph  4.1,  which  we  may  use  to  obtain  m0 . We  now  apply 
Eq.  (4.2.13)  to  obtain  an  estimate  of  variance  and  then  use 
Figs.  4.8  and  4.9  to  determine  confidence  intervals. 

First,  the  beam  (System  2)  is  a one-dimensional  system 
with  modes  of  the  form  sin  kx,  and  the  plate  (System  1)  is 
two  dimensional  with  modes  of  the  form  sin  k^xi  sin  k2X2.  If 
we  regard  the  locations  x,X]_,x2  as  uniformly  located  over 
the  structural  surfaces,  we  find  that 


<t|>*>  <(sinkx)*>  , 

« 5 — a — (one-dimensional) 

<(sinkx)2>2  z (4.3.7a) 

< (sin  kXj^)  *><  (sin  kx2)'*>  9 

<i^>2  <(sinkx1)2><  (sin  kx2)2>2  4 

(two-dimensional)  (4.3.7b) 


Obviously,  the  corresponding  quantity  for  three-dimensional 
systems  is  (3/2) s « 27/8. 

Eq.  (4.3.7)  gives  us  two  of  the  factors  in  Eq.  (4.2.13). 
To  obtain  the  third,  we  must  say  more  about  the  system.  How- 
ever, note  that  the  spatial  variance  alone  contributes  a 
factor  (3/2) 2 (9/4) 2 » 11,  so  that  we  must  have  of  the  order 
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of  10  inodes  interacting  or  sampled  to  bring  <j2/m2  down  to 
the  order  of  1 or  100  inodes  to  get  down  to  a2/m2  * 0.1.  In 
general,  the  higher  the  dimensionality  of  the  system,  the 
greater  is  the  contribution  to  variance  freon  the  spatial 
sampling  factors.  To  counteract  this,  however,  the  higher 
dimensional  systems  normally  have  greater  modal  density,  and 
this  reduces  variance. 

Suppose  that  the  beam  in  Fig.  4.4  is  made  of  aluminum 
and  is  1.3  m long,  .05  m wide,  and  .003  m thick.  The  average 
frequency  separation  between  modes  for  this  beam  is  according 
to  (4.1.23) 


<5fb  = cfb/Jl  * 130  Vf/1000  . (4.3.8) 

Thus,  at  1000  Hz,  the  average  separation  is  130  Hz  and  at 
250  Hz  it  is  65  Hz.  We  assume  that  the  plate  is  also  of 
aluminum,  is  .003  m thick,  and  has  an  area  of  1.7  m2.  The 
average  frequency  separation  for  plate  modes  is  then 


6f  * hc£//3X  * 5HZ  . (4.3.9) 

Also,  a reasonable  value  for  loss  factor  for  both  systems  is 
n1»n2=*0i*  Thus, 
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(4.3.10) 


Let  us  assume  that  our  data  is  taken  in  octave  bands,  so  that 
Af=f/  /2  . Thus,  combining  Eqs.  (4.3,  8,  9,  10),  we  get 


{nln2  2 
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(4.3.11) 
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Thus,  at  f*100  Hz,  this  factor  is  1 and  at  10,000  Hz 
it  is  10“  . The  product  of  this  modal  count  term  and  the 
spatial  factors  is  presented  in  Fig.  4.9  for  the  example. 

We  have  also  sketched  the  variance  expected  when  there  is 
no  uncertainty  in  modal  amplitude  on  the  beam,  reducing  the 
variance  by  a factor  of  (2/3) 2,  from  Eg.  (4.3.7a).  This  is 
appropriate  if  for  example  we  make  response  measurements  at 
the  free  end  of  the  beam  and  connect  the  plate  to  the  base 
of  the  beam. 

Now  we  can  use  this  estimate  of  variance  to  establish 
a mean  bracketing  confidence  interval  in  each  frequency  band. 

If  this  is  done  for  a confidence  coefficient  of  0.8,  we 
obtain  the  result  shown  in  Fig.  4.11.  We  have  used  the 
estimate  for  the  mean  from  Fig.  4.5  in  the  case  where 
nbp  < nb'  wh*cl1  is  not  quite  the  case  for  one  example  since 
n.  * w/4i  « .05/5  a .oi  » n$j*  The  confidence  interval  in 
nopway  depends  on  the  estimate  of  the  mean,  however. 

Altogether  then,  our  analysis  of  variance,  in  combination 
with  a mean  estimate  and  the  calculation  of  confidence  intervals 
based  on  normalized  variance,  provides  us  with  a way  of  making 
estimates  of  response  for  real  systems  that  includes  variations 
in  a realistic  manner.  The  weakest  link  of  this  chain  at 
present  is  the  estimation  of  variance.  There  is  good  reason 
to  believe,  for  example,  that  the  location  of  resonance 
frequencies  among  the  population  of  "similar"  systems  is  not 
well  described  by  a distribution  that  assumes  that  modes  occur 
with  equal  probability  along  the  frequency  axis.  Such  a dis- 
tribution predicts  a high  probability  for  very  closely 
separated  modes,  whereas  the  normal  coupling  among  modes  will 
tend  to  "split"  them,  as  demonstrated  by  Fig.  3.2.  Thus, 
interacting  groups  of  modes  should  have  a lower  probability  of 
very  small  differences  in  resonance  frequencies,  and  this 
should  result  in  smoother  distributions  of  response  and  lower 
variance. 


4.4  Coherence  Effects  - Pure  Tone  and  Narrow  Band  Response. 


In  this  final  section,  we  consider  the  effects  on  our 
analysis  of  the  bandwidth  of  excitation  Aw.  When  Aw  is  very 
broad,  then  we  have  shown  that  certain  frequency  integrations 
such  as  those  of  paragraphs  2.1  and  3.1  become  quite  simple, 
and  that  the  expected  variance  in  system  response  calculated 
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in  paragraph  4.2  is  reduced  as  the  noise  bandwidth  increases. 
When  Aw  gets  smaller  than  the  average  spacing  between  modes, 
we  may  expect  that  variance  will  increase  because  of  the 
smaller  number  of  modes  included  in  the  averages.  When  Aw 
becomes  of  the  order  of  a modal  bandwidth,  other  complications 
(or  possibly  simplifications)  may  develop. 

Narrow  band  excitation  of  resonator.  Let  us  suppose 
that  the  excitation  I (t)  of  the  single  ciof  system  of  Fig.  2.1 

is  the  pure  tone  signal  L e”lwt.  The  response,  according  to 
Eq.  (2.1.20)  is 


v (t) 


. LY 
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(4.4.1) 


In  Chapter  2,  we  averaged  j Y ) 2 over  an  interval  Aw  (large 
compared  to  w0n)  to  determine  the  noise  response  of  the 
resonator.  Let  us  now  suppose  that  w is  fixed  and  we  are 
selecting  resonators  out  of  a population  having  resonance 
frequencies  Wg  uniformly  distributed  over  the  interval  Aw 
(this  is  one  feature  of  the  SEA  model  introduced  in 
paragraph  3.2).  Then  the  average  m.s.  response  will  be 
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where  iu/2) wn,  the  effective  bandwidth.  This  is  quite  a 
slowly  varying  function  of  w,  so  that  if  we  now  say  that  Jt(t) 
is  a narrow  band  signal  of  bandwidth  Sw  so  that  1/2  Lz-*-2Si  (w)  6w , 
we  have  (remember  that  spectral  densities  in  w are  defined  from 

“ 00  < w < + ® ) . 


2S0  (w)  (6w/Aw) 

<v2>  » ^ (4.4.3) 
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which  should  be  compared  to  Eg.  (2.1.39).  We  may  say  that 
the  effect  of  restricting  the  bandwidth  to  Aid,  is  to  reduce 
the  effective  spectral  density  of  the  excitation  in  the  ratio 
<5w/Aw.  When  6 w+Aw,  we  obtain  the  original  response  of  the 
resonator  to  the  noise  of  bandwidth  Aid.  The  response  of  a 
resonator  to  a signal  of  any  bandwidth  is  of  the  same  form  as 
a noise  response/  if  we  average  over  the  resonance  frequency 
of  the  resonator. 

Interacting  Modes.  The  result  of  Eq.  (4.4.3)  allows  us 
to  interpret  the  result  of  Eq.  (3.1.32)  in  an  interesting 
way.  The  interaction  of  mode  "2"  with  a single  mode  "1"  was 
found  to  be  equivalent  to  excitation  of  mode  by  a white  noise 
source.  Let  us  re-examine  this  result  in  the  light  of  Eq. 
(4.4.3).  To  do  this,  imagine  for  simplicity  that  the 
resonators  in  Fig.  3.1  have  gyroscopic  coupling  only.  Then, 
noise  excitation  with  spectral  density  Sj^  of  resonator  1 with 
y2=0,  results  in  a m.s.  velocity  of  this  resonator  given  by 
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This  m.s.  velocity  has  an  equivalent  bandwidth  Ag  . * i\A^/2, 

since  its  spectral  form  is  that  of  the  admittance  6r  resonator  1. 

From  Eq.  (3,1.3b)  a m.s.  velocity  <v£>  will  produce  a m.s. 
force  on  resonator  2 given  by  <v£>  G2  = <v|>  y2M^M2.  This  m.s. 
force  is  "spread"  over  a bandwidth  Ae/i,  so  that  its  spectral 
density  is 


Thus 


2 (oi)  - <v*>  Y*  «iM2/4le> 


ttS.  (w) 

*1 


AiMi 


y2m1m2 

V/i"~ 


(4.4.5) 
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is  the  spectral  density.  Since , according  to  Eq.  (4.4.3),  such 
a signal  is  equivalent  to  a broad  band  noise  over  a band  Aw, 
with  spectrum  Aei  /Aw,  the  equivalent  broad  band  excitation 
produced  by  the  interaction  is 


S 


M2  ttA2/2 
S*1  Ab) 


(4.4.7) 


where  A*y2/A1A2  for  gyroscopic  coupling  only,  which  is  con- 
sistent with  Eg.  (3.1.32)  for  weak  coupling. 

We  have  shown  two  important  features  of  the  interaction 
with  this  calculation.  First,  we  can  properly  account  for 
the  interaction  when  the  coupling  is  weak  (X<<1)  as  excitation 
of  the  indirectly  excited  system  by  noise,  through  a filter 
of  bandwidth  Ale  and  transfer  magnitude  determined  by  the 
coupling  strength,  if  the  resonance  frequency  of  the  receiving 
system  is  allowed  to  take  random  values  over  the  interval  Aw. 
Secondly,  if  the  effective  bandwidth  of  the  excitation  of 
system  2 is  Ae]_,  then  if  more  than  one  mode  may  be  excited 
simultaneously,  we  may  expect  important  coherence  effects  to 
arise  from  such  narrow  band  excitation.  We  now  discuss  just 
what  those  effects  may  be. 

Coherence  Effects  in  Multi-Modal  Response  to  Band-Limited 
NoiseT  We' suppose  that  the  system  described  by  Eq.  (2.2.1)  is 
excited  by  noise  source  of  bandwidth  Aw.  This  source  may  be  a 
vibration  exciter  or  a damped  resonator  driven  by  noise  as  dis- 
cussed in  the  preceding  paragraphs.  The  displacement  response 
is  of  the  form  given  in  Eq.  (2.2.3).  The  density  average  of  m.s. 
displacement  is,  of  course 


<y 


</*(  t) 
m 


(4.4.8) 


We  now  assume  that  all  modes  in  the  interval  ow  have  the  same 
damping,  with  effective  bandwidth  Ae.  If  we  excite  this  system 
at  a location  where  all  modes  have  nearly  the  same  amplitude 
(on  a free  plate,  this  would  correspond  to  a corner),  then  we 
can  assume  that  the  m.s.  response  of  each  mode  is  identical. 
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Accordingly, 


<y2>P,t  - N ^2>t 


(4.4.9) 


where  N is  the  number  of  resonantly  excited  modes  and  <y2>t 
is  the  m.s.  response  of  each. 

Eq.  (4.4.9)  gives  the  average  response  of  the  system  in 
the  spirit  of  paragraph  4.1,  assuming  completely  incoherent 
response.  Let  us  suppose,  however,  that  the  response  of 
every  mode  were  perfectly  coherent,  as  it  would  be  under 
sinusoidal  excitation  (Aw<<Ae) . In  this  event,  the  m.s. 
response  is  at  any  point  is 


<y2>t  = 21  < Ym(t)  Yn(t)>t  ^m(x)^n(x).  (4.4.10) 

m,n 

If  there  is  a location  where  every  mode  has  an  antinode 
(if/=Vmax'  and  every  modal  vibration  is  in  phase,  then  the  m.s. 
response  at  that  location  is 


max  <y2>t  = N2.  <y2>t  i|;2ax  (4.4.11) 

The  ratio  of  the  rms  response  at  this  location  to  the  incoherent 
rms  response  is 


y 

•'max 

yrms 


= *, 


max 


✓ft 


(4.4.12) 


Thus,  if  N is  large,  this  coherent  "hot  spot"  of  response 
may  be  substantially  greater  than  the  incoherent  estimate. 
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Let  us  now  assume  that  the  noise  bandwidth  Aw  is  large 
compared  to  Ae.  The  modal  density  of  the  system  is  n(w). 

Thus,  the  size  of  each  group  of  coherent  excited  modes  in  nAe, 
which  by  presumption  must  be  large  for  coherent  effects  to  be 
of  any  importance*  The  coherent  response  maximum,  therefore, 
according  to  Eq.  (4.4.11),  is  (nAe) 2 <Y2>  t^max*  / ***3  a m*s* 
response  of  nAe<Y2>t,  according  to  Eq.  (4.4.9).  Each  possible 
coherent  peak,  therefore,  exists  in  an  incoherent  "background" 
given  by  n(Aw-Ae),  so  that  the  total  response  at  a hot  spot  is 


max  <y2>t  * (nAft) ^max*****  t+n(Aw“Ae*  <y2>t 

(4.4.13) 


The  ratio  of  maximum  to  average  response  then  is 


max<y2> 

^2>x,t 


* nV max  ; Au<Ae  (4.4.14b) 

As  an  example,  consider  the  plate  of  the  example  in 
section  4.3.  In  that  example  nAe*(Tr/2)  fnp/6fp*10”  Trf/10. 

At  1000  Hz,  the  number  of  coherent  modes  would  be  ^3.  At 
10,000  Hz,  we  would  have  31  modes  in  each  coherent  group. 

The  statistical  response  concentrations  for  this  plate  under 
1/3  OB  and  pure  tone  response  (V;  «=2  for  2-dimensional  mode) 

at  1000  Hz  are 


*+  (1-  -jS)  ,•  Aw>A< 


(4.4.14a) 


<y2> 


10”2Tt*103  N , 

r 4+  (1 TF75 ) =1 
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max<y2>. 

(pure  tone)  = tr • 4 = 13 


<y2> 


At  1000  Hz,  there  is  no  discernible  response  concentration 
due  to  coherence  for  1/3  OB  excitation,  but  a concentration 
factor  of  about  3.5  will  occur  for  pure  tone  excitation. 

At  10,000  Hz,  we  have 


„/3  maX<Y  _ 7T2  *102  *4  ^ ..  tt*10“2  *103» 
(1/3  OB)  - wTT-rr + (1 ) 


<y2> 


- 1.6  + 1 * 2.6 


max<y2>. 

(pure  tone)  = it  *10*4  = 126 

<y:>  * 


The  concentration  for  the  band  of  noise  is  still  negligible, 
but  for  pure  tone  excitation,  the  response  concentration 
is  greater  than  11. 

It  is  clear  that  statistical  response  concentrations 
are  sizeable.  We  have  not  discussed  the  probability  that 
such  a coherent  peak,  as  represented  by  Eq.  (4.4.12),  will 
occur.  Although  the  point  is  not  completely  resolved,  it 
appears  that  such  concentrations  have  a good  likelihood  of 
occurring.  As  a practical  matter,  such  coherence  effects 
will  be  important  when  large  systems  with  high  modal  density 
and  moderate  to  high  damping  are  excited  by  pure  tones.  When 
these  "statistical"  concentrations  occur,  they  may  be  con- 
siderably larger  than  the  "stress  concentrations"  that  occur 
near  boundaries  as  exemplified  by  Fig.  4.3. 
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(0)  MODE  SHAPE  %(*2)  FOR  CLAMPED -FREE  BEAM 


(i>)  M.  S.  DISPLACEMENT  OF  CLAMPED  - FREE  BEAMS&*  (x*) 


(c)  M.S.  SURFACE  STRAIN  OF  CLAMPED- FREE  BEAM { ^*(x2) 


F!'<3  4.1 

SPATIAL  DISTRIBUTION  OF  RESPONSE  FOR  CLAMPED- FREE  BEAM 


k sir.  4>  k, 

FIG.  4.2 


DISTRIBUTION  OF  2-DIMENSIONAL  MODES  EXCITED  BY  BAND  OF  NOISE  Aw 

! 

1 


1 


FIG.  4.3 

MULTI-MODAL  TEMPORAL  MEAN  SQUARE  RESPONSE  OF  RECTANGULAR  PLATE 
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FIG  4.4 

BEAM-PLATE  SYSTEM  .BEAM  INDIRECTLY  EXCITED 


FIG  4.5 

COMPARISON  OF  MEASURED  RESPONSE  RATIOS  WITH  THOSE 

PREDICTED  BY  EO.  (4.1.25) 
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NEUTRAL  AXIS 


FIG  4.6 

BENDING  STRAIN  IN  A BEAM  OF  THICKNESS  h 


a! 


Eor^oco- 

yy  yvyv 


A 


U»- 


FIG  4.7 

THE  SUM  IN  EQ.(4.2.5)  AS  A FUNCTION  OFo>o* 
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FIG  4.8 


UPPER  BOUND  OF  SIMPLE  EXCEEDAf 
AS  A FUNCTION  OF  NORM* 


FIG  4.9 

HALF-WIDTH  OF  MEAN -BRACKETING  ESTIMATION  INTERVALS 
AS  A FUNCTION  OF  NORMALIZED  VARIANCE 
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K)  log 


FREQUENCY  (Hz) 

FIG  4.11 


MEAN  BRACKETING  80%  CONFIDENCE  INTERVAL  BASE  ON  FIG  4.9 
FOR  SYSTEM  OF  FIG 4.4,  MEAN  VALUE  ESTIMATE  OF  FIG 4.5 
ANO  VARIANCE  OF  FIG  4. 10  (DASHED  LINE ) 
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power  flow  between  coupled  resonators  is  discussed  in 

(bw)  D.  E.  Newland,  "Calculation  of  Power  Flow  Between  Coupled 
Oscillators",  J.  Sound  Vib. , Vol.  3,  No.  3,  pp  262-276 
(1966). 

A brief  summary  of  the  problems  of  measuring  various  SEA 
parameters  is  given  in 

(bx)  R.  H.  Lyon,  "Analysis  of  Sound-Structural  Interaction  by 
Theory  and  Experiment",  Contribution  to  Proceedings  of 
Purdue  Noise  Control  Conference,  July  14-16,  1971, 

Purdue  University. 

Modal  density  measurements  of  cylinders  are  presented 
in  (bb,  bs) . Measurements  for  flat  plates  are  given  in 
(ar,  at) , and  for  mass  loaded  plates  in 

(by)  R.  W.  Sevy  and  D.  A.  Earls,  "The  Prediction  of  Internal 
Vibration  Levels  of  Flight  Vehicle  Equipments  Using 
Statistical  Energy  Methods,"  U.S.  Air  Force  Technical 
Report  AFFDL-TR-  69-54,  January  1970. 

Calculations  of  modal  densities  are  given  in  (ba  through 
bh)  for  a variety  of  structural  elements  or  subsystems. 
Input  impedances  to  beams  and  plates  are  found  in  (ao, 
ay,  av) . Input  impedances  to  sound  fields  may  be  found 
in  (am,  aw,  ax,  bd,  bg) , and  in 

(bz)  R.  H.  Lyon,  "Statistical  Analysis  of  Power  Injection  and 
Response  in  Structures  and  Rooms"  J.  Acoust.  Soc.  Am., 

Vol.  45,  No.  3,  pp  545-565  (1969). 

The  derivation  of  transmission  loss  and  its  relation  to 
coupling  loss  factor  may  be  found  in 

(ca)  I.  L.  Ver  and  C.  I.  Holmes,  "Interaction  of  Sound  Waves 
with  Solid  Structures",  Chapter  11  of  Noise  and 
Vibration  Control,  Ed.  by  L.  L.  Beranek  (McGraw-Hill 
Book  Co. , Inc.  New  York  1971). 

(cb)  C.  I.  Holmes,  "Sound  Transmission  Through  Structures: 

A Review"  M.Sc.  Thesis,  John  Carroll  University, 

Cleveland,  Ohio  1969. 
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A useful  "dictionary"  of  transmission  loss  for  various 
practical  wall  structures  may  be  found  in 

(cc)  R.  b.  Berendt,  G.  E.  Winzer  and  C.  B.  Burrough,  "A 
Guide  to  Airborne,  Impact  and  Structure  Borne  Noise 
Control  in  Multifamily  Dwellings."  U.  S.  Dept,  of 
Housing  and  Urban  Development,  Washington,  D.  C., 
September  1967. 

The  relation  of  vibrational  velocity  to  stress  in  a 
variety  of  structural  elements  is  presented  in 

(cd)  f.  V.  Hunt,  "Stress  and  Strain  Limits  on  the  Attainable 
Velocity  in  Mechanical  Vibration",  J.  Acoust.  Soc.  Am., 
Vol.  32,  No.  9,  pp  1123-1128  (1960) " 

(ce)  E.  E.  Ungar,  "Maximum  Stresses  in  Beams  and  Plates 

Vibrating  at  Resonance",  J . Eng . Ind . Vol.  84,  No.  1, 
pp  149-155  (February  19627^  ' 

Analysis  of  response  variance  is  developed  in  (aq>  ar, 
as,  bz) . A basic  discussion  of  confidence  intervals 
and  coefficients  may  be  found  in 

(cf)  A.  M.  Mood,  Introduction  to  the  Theory  of  Statistics 
(McGraw-Hill  Book  Co.,  Inc.,  New  York,  1950).  chapter  11. 


1 . 3 Future  Developments  of  SEA 

The  range  of  professional  workers  making  use  of  SEA  is 
discussed  in  (aa) . That  reference  prompted  additional 
discussion. 

(eg)  J.  L.  Zeman  and  J.  L.  Bogdanoff,  Letter  to  the  Editor, 

The  Shock  and  Vibration  Digest,  Vol.  3,  No.  1,  p 2 (1971). 

(ch)  C.  T.  Morrow,  "Can  Statistical  Energy  Analysis.be 
Applied  to  Design?"  The  Shock  and  Vibration  Digest, 

Vol.  3,  No.  5,  Pt.  1 TT57T). 

and  a reply 

(ci)  R.  H.  Lyon,  Letter  to  the  Editor.  The  Shock  and  Vibration 
Digest,  Vol.  3,  No.  10,  pp  3-4  (1971) . 
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A computer  solution  of  simultaneous  equations  in  SEA 
has  been  carried  out  and  is  reported  in 

(cj)  D.  M.  Wong,  "Transmission  of  Noise  and  Vibration  Between 
Coupled  Cylinders,"  Lockheed  Sunnyvale  Report,  1969. 

The  development  of  a range  of  models  more  complex  than 
the  SEA  models  presented  here  that  have  greater 
"information  capacity"  has  been  discussed  in 

(ck)  R.  H.  Lyon,  "Application  of  a Disorder  Measure  to  Acoustical 
and  Structural  Models"  J.  Eng.  Ind.,  Trans.  ASME,  Vol.  93, 
Ser.  B,  No.  3,  pp  814-8TE  (August  T97TT7 

The  allowance  for  modal  interference  effects  the  modal 
resonance  frequency  spacing  has  been  discussed  in  (al, 
bz) . 


1 .4  Organization  of  the  Report 


Various  surveys  of  SEA  have  been  published  for  various 
audiences.  They  include  (aa,  ao,  at,  bq,  br,  bx)  and 
also 

(cl)  R.  H.  Lyon  and  G.  Maidanik,  "Statistical  Methods  in 
Vibration  Analysis",  AIAA  Journal,  Vol.  2,  No.  6, 
pp  1015-1024  (1964). 

(cm)  P.  A.  Pranken  and  R.  H.  Lyon,  "Estimation  of  Sound- 
Induced  Vibrations  by  Energy  Methods,  with  Applications 
to  the  Titan  Missile",  Shock  and  Vibration  Bulletin, 

No.  31,  Part  III,  pp  12-16  (1963). 

(cn)  R.  H.  Lyon,  "Basic  Notions  in  Statistical  Energy  Analysis," 
Contribution  to  Synthesis  of  Vibrating  Systems,  Ed.  by 

V.  H.  Neubert  and  J.  P.  Raney,  (ASME,  New  York,  1971) . 

(co'  P.  W.  Smith,  Jr,  and  R.  H.  Lyon,  "Sound  and  Structural 
Vibration",  NASA  Contractor  Report  CR-160,  March  1965. 


Chapter  2.  Energy  Description  of  Vibrating  Systems 
2.0  Introduction 

Additional  discussion  of  energy  variables  in  vibration 
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may  be  found  in  (cb) . The  relation  of  impedance 
functions  to  energy  variables  is  discussed  in 

(cp)  E.  A.  Guillemin,  Introductory  Circuit  Theory , (Chapman 
and  Hall,  Ltd,  New  York  1953)  Chapter  1. 

Discussion  of  the  relation  between  modal  and  wave 
descriptions  may  be  found  also  in  (ag,  am,  co) . 


2.1  Modal  Resonators 

Solutions  for  the  single  dof  vibrator  may  be  found  in 
(ab,  ad,  co,  cp) . Relations  among  various  measures  of 
damping  may  be  found  in  (ay) . Reverberation  time  as 
a measure  of  the  damping  of  sound  waves  is  discussed 
in  (ah,  co) . The  idea  of  viewing  variuos  frequency 
regions  of  response  is  "mass-controlled" , etc.  is 
discussed  in  (ca,  co) . Random  excitation  of  simple 
resonators  is  extensively  discussed  in  (ab,  ac) . The 
notion  of  a white  noise  generator  as  a constant  power 
source  is  in  (am; . The  removal  of  the  divergent  mass 
law  response  of  a mode  to  white  noise  may  be  found  in 

(cq)  R.  H.  Lyon,  "Shock  Spectra  for  Statistically  Modeled 
Structures",  Shock  and  Vibration  Bulletin,  No.  40, 

Part  4,  pp  17-23,  December  1969. 


2.2  Modal  Analysis  of  Distributed  Systems 


The  theory  of  the  vibration  of  distributed  systems  by 
using  modal  expansions  may  be  found  in  (an)  and  also  in 

(cr)  P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical 
Physics  (McGraw-Hill  Book  Col,  Inc.,  New  York,  1953) . 

When  damping  is  not  proportional  to  the  mass  density 
an  irreducible  modal  coupling  results  which  is 
extensively  discussed  in  (bn).  Modal  lattices  are 
used  extensively  in  acoustics  (ag,  ah)  and  in  x-ray 
diffraction  theory  and  solid  state  physics: 

(cs)  L.  Brillouin,  Wave  Propagation  in  Periodic  Structures 
(McGraw-Hill  Book  Co.,  Inc.,  New  York,  1946). 
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An  elementary  discussion  of  energy  and  group  velocity 
may  be  found  in  (ai).  The  expansion  of  the  structure 
admittance  function  as  a sum  of  modal  admittances  is 
given  in  (bz).  The  point  input  resistance  of  an  infinite 
plate  is  derived  in  (ay,  co) . The  use  of  conductance 
measurement  as  a means  of  evaluating  modal  density  was 
first  proposed  in  (bx) . 


2.3  Dynamics  of  Infinite  Systems 


The  development  of  dispersion  relations  and  their  inter- 
pretation in  terms  of  phase  and  group  velocities  may 
be  found  in 

(ct)  L.  M.  Brekhovskikh , Waves  in  Layered  Media  (Academic 
Press,  New  York  1960j"I  Derivations  for  the  "stiffness 
operator"  for  strings,  beams,  and  sound  fields  may  be 
found  in 

(cu)  P.  M.  Morse,  Vibration  and  Sound,  Second  Edition, 
(McGraw-Hill  Book  Co . , Inc . 1 94 8 ) . The  theory  of 
Fourier  integrals  and  the  use  of  contour  integrals 
to  evaluate  the  inverse  transform  may  be  found  in 
(cr) . The  derivation  of  mean  free  path  for  a two 
dimensional  system  is  to  be  found  in  (co) . The  dis- 
cussion of  resonator  interaction  with  a finite  plate 
is  adapted  from  (bo) . 

2.4  Modal -Wave  Duality 


The  decomposition  of  a wave  field  into  "direct"  and 
"reverberant"  components  is  discussed  in  (ag,  ah,  ay, 
bz,  co) . The  formation  of  the  direct  field  as  a 
result  of  modal  coherence  may  be  found  in  (ag) . 

Various  integral  expressions  for  and  relations  between 
bessel  functions  may  be  found  in  fcr) . 


Chapter  3.  Energy  Sharing  by  Coupled  Systems 


A discussion  of  the  environmental  excitation  of 
structures  may  be  found  in  (bd,  bu)  and  in 
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(cv)  R.  H.  Lyon/  "An  Energy  Method  for  Prediction  of  Noise 
and  Vibration  Transmission",  Shock  and  Vibration 
Bulletin,  No.  33,  Part  II  pp  13-25  (February  1964). 

The  definition  of  uncoupled  systems  has  been  an  important 
element  in  several  studies  (be,  bk,  br,  bs) . An 
interesting  application  of  reciprocity  to  cevelop 
engineering  formulas  for  nedse  transmission  is 

(cw)  M.  A.  Heckl  and  E.  J.  Rathe,  "Relationship  Between 
the  Transmission  Loss  and  the  Impact-Noise  Isolation 

of  Floor  Structures",  J.  Acoust.  Soc.  Am.,  Vol  35,  No.  11 
pp  1825-1830  (1963). 


3.1  Energy  Shari ng  Among  R^ionators 


The  discussion  in  this  section  is  taken  from  (ap) , 
which  in  turn  is  based  on 


(cx)  T.  p.  Scharton,  "Random  Vibration  of  Coupled  Oscillators 
ahd  Coupled  Structures,"  ;•>«-. D.  Thesis,  MIT  Dept,  of 
Mechanical  Engineering,  15)6  5\. 

Additional  discussion  of  the  effects  of  blocking 
vibrating  systems  may  be  found  in  (bq,  bn)  and  in 

(cy)  D.  E.  Newland,  ’’Power  Flow  Between  A Class  of  Coupled 
Oscillators,"  J.  Acoust.  Soc.  Am.,  Vol.  43,  No.  3, 

pp  553-559  (1966')  7 ' 

The  example  of  two  resonators  coupled  by  stiffness  only 
is  taken  from  (cn) . The  behavior  of  coupled  resonators 
in  free  vibration  may  be  found  in  (ad) . The  analysis 
of  power  flow  in  this  section  and  in  (ao)  is  carried 
out  in  the  "frequency  domain".  For  a time  domain 
analysis,  see  (am).  Integrations  over  the  frequency 
response  functions  in  Eq.  (3.1.14)’  may  be  found  in 
(ab) . The  result  of  Eq.  (3.1.19)  was  found  in  (am) 
for  weak  coupling,  although  the  result  (for  weak 
coupling)  was  also  applied  to  relations  like  Eq.  (3.1.25) 
The  essential  correctness  of  Eq.  (3.1.25)  must  be 
found  from  (ao,  ap,  cx) . The  equivalence  of  average 
modal  interaction  and  white  noise  excitation  was  first 
put  forward  in  (am) . 
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3.2  Energy  Exchange  In  Multl-Degrce-of-Frcedom  Systems 

Additional  discussion  of  the  assumptions  involved  in 
forming  the  ’uncoupled"  modes  for  the  blocked  condition 
may  be  found  in  (bn) . A discussion  of  the  treatment 
of  damping  as  a random  parameter  is  to  be  found  in 
(eg) . Quite  a thorough  analysis  of  sound-structural 
coupling,  carried  out  on  ,a  completely  modal  basis  Is 
to  be  found  in 

(cz)  F.  Fahy,  "vibration  of  Containing  Structures  by  Sound 
in  the  Contained  Fluid",  J.  Sound  Vib.,  Vol.  10,  No.  3 
pp  490-512  (1969).  “ 

(da)  F.  Fahy,  "Response  of  a Cylinder  to  Random  Sound  in  the 
Contained  Fluid",  J.  Sound  Vib.,  Vol.  13,  No.  2, 
pp  171-194  (1970). 

Discussions  of  the  generation  of  multi-modal  inter- 
actions from  the  two-mode  interaction  are  to  be  found 
in  (am,  ao,  at).  The  first  published  use  of  a diagram 
like  that  of  Fig.  3.5  appears  to  be  (cv) . 


3.3  Reciprocity  and  Energy  Exchange  In  .the  Wave  Description 

For  a discussion  of  diffusion,  reverberation  and 
energy  density  in  sound  fields,  see  (ah,  cu) . The 
general  conditions  for  reciprocity  in  continuous 
systems  are  given  in  (cr,  ct)  and  for  "lumped" 
systems  in  (cp) . The  use  of  reciprocity  to  develop 
response  estimates  of  the  SEA  form  may  be  found  in 
(cv) . The  use  of  impedance  relations  to  determine 
coupling  loss  factors  may  be  found  in  (am,  as,  au, 
be) . A study  of  two  plates  connected  at  a point  is 
presented  in  (bq) . The  relation  between  coupling 
loss  factor  and  transmission  coefficients  (or 
transmission  loss)  is  developed  in  (bk,  cn)  and  in 

(db)  R,  H.  Lyon,  Aerospace  Noise  Vibration.  Notes  for  a 
-'ecture  series  offered  by  Bolt  Beranek  and  Newman, 
1966. 

Transmission  loss  data  are  to  be  found  in  (cc) . 


3 . 4  Some  Sample  Applications  of  SEA 

This  discussion  of  the  excitation  of  a resonator  by  a 
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sound  field  may  be  compared  to  that  in  (an) . The 
analysis  of  coupled  beams  may  be  compared  with 
similar  discussions  in  (ap;  br)  and 

(dc)  H.  6.  Davies,  "Exact  Solutions  for  the  Response  of  Some 

Coupled  Multimodal  Systems",  J.  Acoust.  Soc.  Am.,  Vol.  51, 
No.  1 (Pt.  2),  pp  387-392  (197TT: 

The  relation  between  "mass  law"  transmission  loss  and 
the  coupling  loss  factor  is  discussed  in  (bk,  db) . 


Chapter  4.  The  Estimation  of  Response  in  Statistical  Energy 
Analysi s . 


4.0  Introduction 


Additional  discussion  of  the  requirements  for  response 
estimation  may  be  found  in  (aa,  ao,  bu,  eg,  ch,  ci) 
and  in 

(dd)  J.  L.  Bogdanoff,  "Meansquare  Approximate  Systems  and 
Their  Application  in  Estimating  Response  in  Complex 
Disordered  Linear  Systems",  J.  Acoust.  Soc..  Am., 

Vol.  38,  No.  2,  pp  244-252  (l9$5) . 

General  source  books  for  statistical  estimation  pro- 
cedures are  (bf)  and 

(de)  H.  Cramer,  Mathematical  Methods  of  Statistics  (Princeton 
University  Press.  Princeton,  194(5).  J“""~ 


4. 1 Mean  Value  Estimates  of  Dynamical  Response 


The  mean  square  pressure  in  reverberant  sound  fields 
has  been  studied  extensively  in 

(df)  R.  V.  Waterhouse,  "Statistical  Properties  of  Reverberant 
Sound  Fields",  J.  Acoust.  Soc.  Am.,  Vol.  43,  No.  6, 
pp  1436-1444  (15^T  . 

Spatial  coherence  effects  in  multi-modal  response  in 
the  neighborhood  of  the  drive  point  are  discussed  in 
(bz)  and  in 
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(dg)  L.  Wittig,  "Random  Vibration  of  Point-Driven  Strings 
and  Plates",  Ph.D.  Thesis,  MIT  Dept,  of  Mech.  Eng., 

May  1971. 

The  problem  of  adequate  sampling  of  a reverberant 
field  for  reducing  uncertainty  in  the  energy  content 
of  the  field  has  been  discussed  by 

(dh)  R.  V.  Waterhouse  and  D.  Lubman,  "Discrete  Versus 
Continuous  Averaging  in* a Reverberant  Sound  Field", 

J.  Acoust.  Soc.  Am.  Vol.  48,  No.  1 (Pt.l)  pp  1-5 

jrmr. 

The  beam  plate  system  studied  in  this  section  is  similar 
to  that  in  (as) . There  has  been  much  confusion  regarding 
the  use  of  infinite  system  impedances  to  evaluate  coupling 
loss  factors  in  finite  systems,  as  was  done  in  (am,  as, 
an,  aw).  Hopefully  the  discussion  in  this  section  will 
clarify  the  procedure.  The  expression  of  system  strain 
in  terms  of  a particle  "Mach  number"  is  found  in  (cd, 
ce) . 


4.2  Calculation  of  Variance  In  the  Temporal  Mean  Square 
Response^ 


The  analysis  of  variance  is  an  extensive  and  important 
topic  in  statistical  theory  (cf,  de) . The  procedures 
used  in  this  section  were  first  developed  in  (ag,  ar, 
as)  and  are  based  on  the  theory  of  variance  of  Poisson 
pulse  processes  as  presented  in 

(di)  S.  6.  Rice,  "Mathematic  Analysis  of  Random  Noise". 
Contribution  to  Noise  and  Stochastic  Processes , Ed. 
by  N.  Wax  (Dover  Publications,  Inc.  New  York,  1954). 

The  effects  of  temporal  variation  in  multi  dof  system 
variance  are  discussed  by 

(dj)  H.  Andres,  "A  Law  Describing  the  Random  Spatial 
Variation  of  Sound  Fields  in  Rooms  and  Its  Application 
to  Sound  Power  Measurements."  Acustica,  Vol.  16, 

p 278-294  (1965).  Translation  by  I.  L.  Ver,  BBN-TIR-70, 
(1968). 

Eq.  (4.2.9)  is  found  from  a relation  in  (di) . Additional 
discussion  of  the  role  of  mode  shapes  and  sampling 
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strategies  for  locations  of  excitation  and  observation 
points  in  affecting  response  variance  may  be  found 
in  (as,  bz).  The  use  of  multiple  observation  points 
or  line  averages  to  reduce  variance  are  discussed  in 
(dh) . Modifications  to  the  variance  calculation  for 
non-Poisson  intervals  between  resonance  frequencies 
are  discussed  in  (bz).  Note  that  the  effective  number 
of  modal  interactions  given  in  Eq.  (4.2.13)  is  not  the 
total  number  of  possible  interactions  N.N,  as  suggested 
in  (cz,  da).  1 1 


4 . 3 Calculation  of  Confidence  Coefficients 


The  approach  and  terminology  in  this  section  are  taken 
from' (cf) . The  gamma  and  related  distributions  are 
discussed  in  (cf,  de)  and  in 

(dk)  E.  Parzen,  Stochastic  Processes  (Holden-Day,  Inc., 

San  Francisco,  1962) . 

The  two  particular  confidence  intervals  derived  in 
this  section  were  first  presented  in  (as,  cv) . The 
incomplete  gamma  function  is  discussed  in 

(dl)  A.  Erdelyi,  et  al.,  Higher  Transcendental  Functions, 

Vol.  2 (McGraw-Hill  Book  Co.,  Inc.  1953).  Chapter  IX. 

Confidence  coefficients  for  the  determination  of 
environmental  levels  based  oh  loads  and  response  data 
are  used  in 

(dm)  P.  T.  Mahaffey  and  K.  W.  Smith,  "A  Method  for  Predicting 
Environmental  Vibration  Levels  in  Jet-Powered  Vehicles", 
S.  & V.  Bulletin,  No.  28,  Pt.  4,  pp  1-14  (1960). 


4.4  Coherence  Effects  - Pure  Tone  and  Narrow  Band  Response 


The  discussion  in  this  section  is  based  on 

(dn)  R.  H.  Lyon,  "Spatial  Response  Concentrations  in  Extended 
Structures".  Trans.  ASME,  J,  Eng,  Ind.,  November  1967. 
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Additional  discussion  of  extreme  value  statistics 
for  multi-dimensional  sinusoids  may  be  found  in 


(do)  R.  H.  Lyon,  "Statistics  of  Combined  Sine  Waves" 
J.  Acoust.  Soc.  Am. , Vol.  48,  No.  1 (Pt.  2) 
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